Reproduced  From 
Best  Available  Copy 


r 


OPTIMAL  SILO  ATTACK  PLAN  FOR  THE  RED 
INTEGRATED  STRATEGIC  OFFENSIVE  PLAN 
(RISOP) 

THESIS 

Richard  Charles  Pace 
Captain,  USA 

_ AFIT/GQR/ENS/93M-15  _ 


DEPARTMENT  OF  THE  AIR  FORCE 
y  AIR  UNIVERSITY 

{  AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


s 


DTIC 

ELECTE 
mo  5  1933 

B 


Wright-PatLer$on  Air  Force  Base,  Ohio 

olD0009,l^  ID4 


AFIT/GOR/ENS/93M-15 


OPTIMAL  SILO  ATTACK  PLAN  FOR  THE  RED 
INTEGRATED  STRATEGIC  OFFENSIVE  PLAN 
(RISOP) 


THESIS 

Richard  Charles  Pace 
Captain,  USA 


AFIT/GOR/ENS/93M-15 


Approved  for  public  release;  distribution  unlimited 


4  02  162 


93-07014 


AFIT/GOR/ENS/93M-15 


OPTIMAL  SILO  ATTACK  PLAN  FOR  THE  RED  INTEGRATED 
STRATEGIC  OFFENSIVE  PLAN  (RISOP) 

THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 
In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Meister  of  Science  in  Operations  Research 


j 

Richard  Charles  Pace,  B.A. 
Captain,  USA 


March,  1993 

jfnC  QUA1IT7  INSPECTED 


HTIS 
DTIC  T,!B 
UnwruiouDced 
Ju'jt :  e.jt  lon_ 


sr 


By . _ 


}_Dl3tMbutlon/_ 
ivallability  Codes 
Avail  and/or 
Special 


Dial 

J/I'l 


Approved  for  public  releaise;  distribution  unlimited 


□  □ 


THESIS  APPROVAL 


STUDENT:  CPT  Richard  C.  Pace  CLASS:  GOR-93M 

THESIS  TITLE:  OPTIMAL  SILO  ATTACK  PLAN  FOR  THE  RED 

INTEGRATED  STRATEGIC  OFFENSIVE  PLAN  (RISOP) 

DEFENSE  DATE:  March  4,  1993 

COMMITTEE: 

Advisor 

Reader 

Reader 


NAME/DEPARTMENT  SIGNATURE 


Lt  Col  James  T.  Moore/ENS 


Dr.  James  W.  Chrissis/ENS 


Maj  Bruce  W.  Morlan/ENC 


A  cknowledge  merits 


I  have  had  a  great  deal  of  help  in  writing  this  thesis.  I  would  especially  like  to 
thank  my  advisor,  I  t  Col  Moore,  for  always  being  available  when  I  held  a  question 
2uid  for  making  sure  I  stayed  on  track.  Thanks  is  also  due  to  my  readers.  Dr.  Chrissis 
and  Maj  Morlan,  for  their  timely  input.  I  am  deeply  indebted  to  my  point  of  contact 
at  J-8  (DISA/DSSO),  CPT  Pete  Davidson,  who  proposed  this  problem  to  begin  with 
and  who  provided  a  tremendous  amount  of  information  cuid  other  help  for  the  last  six 
months.  I  also  owe  thanks  to  Dr.  Roy  Marsten  of  the  Georgia  Institute  of  Technology 
for  making  ZOOM  available  to  me  for  this  research.  Finally,  and  most  critical  of  all, 
my  deepest  gratitude  goes  to  my  wife,  Beth,  who  sacrificed  more  than  I  did  over  the 
last  six  months  in  order  to  ensure  this  project  was  successfully  completed. 

Richard  Charles  Pace 


II 


Table  of  Contents 


Page 

Acknowledgements  . ii 

List  of  Figures  . . vi 

List  of  Tables . vii 

Abstract . viii 

1.  Introduction .  1 

1.1  Background . 1 

1.2  Statement  of  the  Problem .  3 

1.3  Research  Objectives  .  3 

1.4  Scope  . 4 

1.5  Description  of  the  Unclassified  Sample  Problem .  5 

1.6  Approach . .  .  5 

1.7  Format . .  .  6 

IL  Literature  Review  .  .  .  ....  ...  ...  .  . . .  .  .  .  7 

2.1  Introduction . 7 

2.2  Model  Terminology  . . .  ^  7 

2.3  Mode!  Elements .  7 

2.3.1  The  Weapon  System .  8 

2.3.2  The  Target  Complex .  9 

2.3.3  The  Engagement  Submodel .  10 

2.3.4  The  Damage  Submodel .  11 

2.3.5  The  Solution  Algorithm  . 11 

iii 


Page 

2.4  The  J-8  Missile  Allocation  Problem .  12 

2.5  Solution  Techniques .  13 

2.5.1  Goal  Programming  . . 13 

2.5.2  Disjunctive  Constraints  .  .  14 

2.6  Conclusion . . . .  ,  15 

III.  Model  Formulation  and  Methodology .  16 

3.1  Introduction . : .  16 

3.2  Variable  Definitions .  16 

3.3  Model  Formulation . . .  .  18 

3.3.1  Objective  Function .  18 

3.3.2  Constraints . 18 

3.3.3  The  Complete  Model . 30 

3.4  Solution  Methodology  .  33 

3.4.1  Procedure .  33 

IV.  Results  and  Findings . . .  37 

4.1  Introduction . .  37 

4.2  Solution  to  the  Unclassified  Problem .  38 

4.2.1  Allocation  Summary  and  Explanation .  39 

4.3  Model  Verification  and  Validation .  42 

4.3.1  Comparison  to  GAMS/ZOOM .  42 

4.3.2  Face  Validation . 44 

4.4  Solution  Using  Number  of  Warheads .  44 

4.4.1  Allocation  Summary  and  Explanation .  45 

4.5  Parametric  Analysis  .  47 


IV 


Page 

V.  Conclusions  and  Recommendations . 49 

5.1  Conclusions . 49 

5.2  Model  Limitations  . 50 

5.3  Recommendations . 51 

Appendix  A.  Sample  Problem-Unclassified  Data  .  . . 53 

A.l  Launch  Fields  and  Missile  Data . 53 

A. 2  Target  Complexes  and  Target  Data  .  . . •  >  ■  •  53 

A.  3  Flight  Time  Data .  54 

Appendix  B.  FORTRAN  Code .  55 

B. l  Introduction  . . 55 

B.2  Limitations . .  .  55 

B. 3  FORTRAN  Code .  56 

Appendix  C.  Output  Files .  99 

C. l  Solution  Summary  for  J-8  Review .  99 

C.2  Solution!  immary  for  SINBAC . .  .  99 

Appendix  D.  GAMS  Input  File . 100 

Appendix  E.  Complete  Solutions  and  Allocation  Summaries .  103 

Appendix  F.  User’s  Guide . 119 

F.l  Introductio.i .  119 

F.2  How  to  Use  SILOATT .  119 

Bibliography  . 124 

Vita . 126 


V 


\ 


List  of  Figures 


Figure  Page 

1.  Flow  Diagram  .  35 

2.  Flow  Diagram  (Continued)  . . . .  .  36 


VI 


List  of  Tables 

Table  Page 

1.  Sample  Problem  Allocation  Summary . 40 

2.  Allocation  Summary  for  the  GAMS/ZOOM  Solution .  44 

3.  Allocation  Summary  Using  Warheads  Instead  of  Missiles .  46 

4.  Parametric  Analysis  Summary  .  . .  47 

5.  Launch  Fields  and  Missile  Data  . . 53 

6.  Target  Complexes  and  Target  Data .  53 

7.  Flight  Time  Data . .  54 

8.  Case  2  Allocation  Summary . 104 

9.  Case  3  Allocation  Summary . 106 

10.  Case  4  Allocation  Summary .  108 

11.  Case  5  Allocation  Summary . 110 

12.  Case  6  Allocation  Summary . 112 

13.  Case  7  Allocation  Summary . 114 

14.  Case  8  Allocation  Summary  . .  116 

15.  Case  9  Allocation..Summary . 118 


AFIT/GOR/ENS/93M-15 


Abstract 

The  Joint  Staff  Directorate  for  Force  Structure,  Resources,  and  Assessment 
(J-8)  sought  a  procedure  which  could  be  used  to  generate  an  optimal  missile  al¬ 
location  for  the  silo  attack  portion  of  the  Red  Integrated  Strategic  Offensive  Plan 
(RISOP).  Their  current  solution  procedure  is  a  manual  heuristic  which  is  time- 
consuming  and  is  not  guaranteed  to  lead  to  an  optimal  solution.  J-8  defines  an 
optimal  solution  as  a  feasible  solution  which  minimizes  both  the  flight  time  of  the 
missile  that  impacts  first  and  the  duration  of  the  attack.  J-8  defined  several  input 
rules  which  limit  how  missiles  may  be  allocated,  k  mathematical  model  of  the  J-8 
missile  allocation  problem  was  developed  that  uses  a  goal  programming  approach 
to  solve  the  problem.  The  input  rules  defined  the  constraints  for  the  proble.  J-8 
provided  unclassified  sample  data  tc  use  as  a  test  case.  The  model  was  used  to  solve 
the  sample  problem  with  nine  different  variations  of  the  data.  The  model  developed 
is  a  flexible  tool  designed  to  solve  missile  allocation  problems  whose  objective  is  to 
minimize  the  first  impact  time  or  minimize  the  duration  of  the  attack.  The  model 
uses  binary  variables,  so  it  may  be  impractical  to  use  if  the  number  of  binary  vari¬ 
ables  gets  large.  However,  in  the  foreseeable  future,  the  size  of  the  problem  should 
decretise,  thus  making  the  model  usable  for  at  least  several  years. 
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OPTIMAL  SILO  ATTACK  PLAN  FOR  THE  RED  INTEGRATED 
STRATEGIC  OFFENSIVE  PLAN  (RISOP) 

1.  Introduction 

1.1  Background 

The  problem  of  allocating  available  warheads  to  targets  has  existed  since  the 
development  of  long-range  missiles.  This  problem  has  received  much  greater  atten¬ 
tion  with  the  development  of  long-range,  nuclcar-armed  missiles.  As  the  build-up 
of  nuclear  •  lissi’e  arsenals  in  the  United  States  and  the  Soviet  Union  progressed,  so 
did  the  urgency  placed  on  finding  optimal  solutions  to  missile  allocation  problems 
(MAPs). 

Today,  with  the  dissolution  of  th6  Soviet  Union,  the  prevailing  sentiment  may 
be  that  this  problem  is  irrelevant.  However,  as  long  as  nations  jiossess  nuclear 
r’issiles,  the  threat  of  nuclear  attack,  although  perhaps  greatly  diminished,  still 
exists. 

Currently,  the  Joint  Staff  Directorate  for  Force  Structure,  Resources,  and  As¬ 
sessment  (J-8)  at  the  Pentagon  generates  the  Red  Integrated  Strategic  Offensive  Plan 
(RISOP),  which  is  used  to  validate  the  Single  Integrated  Operations  Plan  (SIOP) 
(Davidson,  1992a:2).  The  SIOP,  developed  by  U.  S.  Strategic  Command,  is  the  con¬ 
tingency  plan  the  United  States  has  for  the  case  of  all-out  nuclear  war.  The  RISOP 
is  the  scenario  which  is  believed  to  be  the  woist-case  attack  an  enemy  could  conduct 
against  the  United  States. 

In  order  to  generate  the  RISOP,  J-8  uses  a  computer  model  called  SINBAC 
(System  for  Integrated  Nuclear  Battle  Analysis  Calculus).  SINBAC  determines  the 


optimal  allocation  of  nuclear  warheads  by  maximizing  the  expected  amount  of  dam¬ 
age  that  attack  could  inflict  on  the  United  States.  Due  to  the  target  values 
assigned,  SINBAC  does  not  allocate  any  weapons  for  attacking  U.  S.  missile  silos. 
Since  J-8  realizes  that  any  adversary  wishing  to  prevent  a  U.  S.  counterstrike  will 
attack  U.  S.  missile  silos,  they  have  added  a  routine  lo  SINBAC  that  can  be  used  to 
analyze  a  user-specified  silo  attack  plan. 

•In  order  to  generate  the  user-specified  missile  silo  attack  portion  of  the  RISOP, 
J-8  must  solve  a  missile  allocation  problem.  For  this  MAP,  J-8  designates  in  advance 
which  groups  of  missiles  from  which  launch  fields  will  be  used  for  the  attack.  These 
missiles  are  chosen  because  J-8  knows  that  they  achieve  the  desired  level  of  damage  on 
the  targets  (missile  silos).  Thus,  for  this  portion  of  their  MAP,  they  are  not  concerned 
with  maximizing  the  amount  of  damage  done.  Instead,  J-8  has  two  objectives  for 
their  MAP;  1)  minimize  the  flight  time  of  the  missile  that  impacts  first;  and  2) 
minimize  the  duration  of  the  attack.  For  this  MAP,  “duration  of  the  attack"  is 
defined  as  the  difference  between  the  time  of  the  first  warhead  impact  and  the  time 
of  the  last  warhead  impact.  It  is  assumed  that  the  attacker  would  want  to  minimize 
the  flight  time  of  the  first  missile  in  order  to  attack  the  targets  as  soon  as  possible. 
It  is  assumed  that  the  attacker  would  want  to  minimize  the  duration  of  the  attack  in 
order  to  limit  the  reaction  lime  of  the  target  complexes  after  the  first  complex  has 
been  attacked. 

Currently,  J-8  uses  a  manual  heuristic  to  allocate  enemy  missiles  against  U.  S. 
missile  silos  (Davidson,  1992a:2).  A  heuristic  is  a  method  for  finding  a  feasible 
solution  that  is  not  guaranteed  to  be  optimal.  The  heuristic  used  by  J-S  basically 
consists  of  visually  examining  the  flight  times  by  missile  type  from  each  launch  field 
to  each  target  complex  and  assigning  warheads  to  targets  in  such  a  way  as  to  attempt 
to  strike  each  target  as  early  as  possible.  This  procedure  is  done  manually  and  is 
time-consuming. 
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l.S  Statement  of  the  Problem 

J-8  seeks  a  procedure  which  can  be  used  to  generate  an  optimal  missile  cillo- 
cation  plan;  however,  their  current  heuristic  is  not  guaranteed  to  lead  to  an  optimal 
solution.  J-8  defines  an  optimal  solution  <is  a  feasible  solution  which  minimizes  both 
the  flight  time  of  the  missile  that  impacts  first  and  the  duration  of  the  attack  (David¬ 
son,  1992b:2).  Simply  put,  the  problem  is  ..o  develop  a  procedure  which  determines 
an  optimal  missile  allocation  plan.  The  current  heuristic,  which  is  time-consuming 
and  docs  not  necessarily  generate  an  optimal  solution,  needs  to  be  replaced. 

l.S  Resear^'h  Objectives 

The  major  objectives  of  the  research  are  two-fold.  One  objective  is  to  develop 
a  procedure  which  determines  the  optimal  missile  allocation  plan  and  the  other  is  to 
meet  the  needs  of  J-8.  As  a  minimum,  J-8  needs  a  description  of  how  the  problem 
may  be  solved  using  existing  commercially  available  software  as  long  as  it  can  run 
on  VAX  or  SUN/UNIX  machines.  Given  this  description,  they  want  a  structured 
algorithm  describing  how  the  problem  may  be  solved.  Finally,  J-8’s  greatest  need 
is  for  a  FORTRAN  routine  that  can  be  incorporated  into  SINBAC  that  determines 
the  optimal  silo  attack  plan  (Davidson,  l992a:2). 

Unfortunately,  finding  the  “optimal”  solution  is  not  as  simple  as  it  may  sound. 
This  is  due  to  there  being  two  objectives  for  the  solution  of  J-8’s  MAP  that  conflict 
somewhat:  the  shorter  the  flight  time  of  the  first  missile,  the  longer  the  duration  of 
the  missile  attack;  and  the  shorter  the  duration,  the  longer  the  flight  time  of  the  first 
missile.  The  “optimal”  solution  will  therefore  have  a  flight  time  of  the  first  missile 
that  is  perhaps  greater  than  the  minimum  possible  flight  time,  and  a  duration  that 
is  perhaps  greater  than  the  minimum  that  could  possibly  be  achieved.  To  resolve 
this  conflict,  the  decision  maker  must  specify  how  short  the  duration  should  be  in 
relation  to  the  minimum  flight  time.  To  do  this,  the  decision  maker  must  define 
a  percentage  of  the  minimum  flight  time  that  the  duration  is  allowed  to  be.  For 
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instance,  if  the  decision  maker  wants  the  duration  to  be  no  longer  than  10  percent  of 
the  minimum  flight  time,  then  if  the  flight  time  of  the  first  missile  was  40  minutes, 
the  duration  should  be  no  longer  than  four  minutes.  In  addition,  the  decision  maker 
must  define  weighting  factors  associated  with  both  the  flight  time  of  the  missile  which 
impacts  first  and  the  duration  of  the  attack  which  indicate  the  relative  importance 
of  minimizing  each  of  them. 

1.4  Scope 

In  order  to  keep  this  document  unclassified,  J-8  has  modified  the  data.  The 
launch  fields  and  target  complexes  are  fictitious.  The  data  for  number  of  missiles, 
number  of  targets,  distances,  and  flight  times  are  also  notional.  In  addition,  missile 
types  are  identified  only  as  “good”  and  “fair,”  and  targets  are  not  identified  but  only 
classified  as  “good,”  “fair,”  or  “poor.” 

The  specific  missile  allocation  problem  that  J-8  needs  to  solve  has  15  allocation 
rules  and  restrictions  which  more  precisely  define  the  problem.  A  complete  discus¬ 
sion  of  these  rules  and  restrictions  and  their  mathematical  modeling  is  presented  in 
Chapter  III.  Below,  the  most  restrictive  rules  are  listed: 

1.  A  wave  (set  of  missiles  launched  at  the  same  time)  coming  from  one  launch 
field  is  better  than  a  wave  launched  from  two  different  launch  fields. 

2.  For  fair  missiles,  it  is  better  to  target  a  given  complex  with  two  waves  from  the 
same  launch  field  than  with  one  wave  from  each  of  two  different  launch  fields. 

3.  In  the  initial  wave,  good  targets  must  be  hit  before  fair  targets,  and  fair  targets 
must  be  hit  before  poor  ones. 

4.  The  attackers  will  keep  at  least  10  percent  of  their  missiles  as  backup  (strategic 
reserve).  A  “backup”  missile  is  one  that  is  not  allocated  for  launch  in  either 
the  initial  or  follow-on  wave.  It  is  best  to  have  10  percent  of  the  missiles  in 
each  launch  field  kept  as  backup. 
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5.  No  “strays”  are  allowed.  This  means  that  the  cillocation  should  be  somewhat 
balanced.  In  other  words,  the  attackers  will  not  allocate  90  percent  of  a  given 
target  complex’s  targets  to  missiles  from  one  launch  field  and  10  percent  to 
missiles  from  another.  The  allocation  will  be  closer  to  50  percent  allocated  to 
each  launch  field. 

6.  During  each  wave,  every  target  in  a  given  complex  must  be  hit  at  the  same 
time.  This  may  result  in  some  missiles  not  being  launched  at  the  same  time 
as  others.  In  effect,  all  missiles  allocated  to  a  given  target  complex  have  the 
same  time  of  impact. 

1.5  Description  of  the  Unclassified  Sample  Problem 

The  data  for  the  sample  problem  is  given  in  Appendix  A.  Some  of  the  launch 
fields  have  only  one  type  of  missile  and  the  others  have  both  types.  Each  target 
complex  has  only  one  type  of  target.  However,  J-8  has  staled  that  Tybee  and 
Atlanta  are  close  enough  together  to  be  treated  as  a  single  complex.  If  so  combined, 
the  resulting  combined  complex  has  both  good  and  fsur  targets. 

The  launch  fields  and  target  complexes  are  somewhat  aggregated  in  that  each 
missile  of  a  given  type  from  a  given  launch  field  has  the  same  flight  time  to  each 
target  complex.  There  is  no  data  for  the  flight  times  of  each  individual  missile  to 
each  individual  target.  All  flight  time  data  is  from  the  center  of  mass  (centroid)  of 
the  launch  field  to  the  center  of  mziss  (centroid)  of  the  target  complex. 

1.6  Approach  \ 

This  problem  was  formulated  as  a  mixed-integer  program  (MIP).  The  MIP’s 
objective  function  is  to  minimize  the  flight  time  of  the  missile  that  impacts  first  and  to 
minimize  the  duration  of  the  attack.  The  other  input  rules  establish  the  parameters 
and  constraints  of  the  problem.  In  formulating  the  problem,  a  constraint  which 


ensures  that  the  duration  is  no  longer  than  a  given  percentage  of  the  minimum  flight 
time  is  needed. 

Once  the  problem  was  formulated,  a  procedure  was  developed  that  solves  it. 
This  procedure  was  implemented  in  a  FORTRAN  program  which  required  four  major 
parts:  one  that  reads  in  all  the  appropriate  data,  a  second  that  processes  the  data 
and  formulates  the  mixed-integer  program  (see  Chapter  III),  a  third  that  solves  the 
problem,  and  a  fourth  that  produces  the  solution  in  a  format  that  J-8  can  use. 
Ideally,  this  approach  allows  for  any  number  of  launch  fields,  missiles,  missile  types, 
target  complexes,  targets,  and  target  types.  In  the  solution  phase,  the  mixed  integer 
solver  called  ZOOM  (Zero/One  Optimization  Methods)  was  used.  This  solver  was 
chosen  because  it  is  written  in  FORTRAN  and  may  be  included  in  the  program  as 
a  set  of  subroutines. 

1.7  Format 

In  Chapter  II,  much  of  the  literature  pertaining  to  missile  allocation  problems 
(MAPs)  is  reviewed,  terminology  and  components  common  to  missile  allocation  prob¬ 
lems  are  described,  and  solution  techniques  that  are  used  to  solve  the  problem  are 
discussed.  In  Chapter  III,  the  problem  formulation  and  the  method  used  to  solve 
this  problem  are  discussed.  In  Chapter  IV,  the  results  and  findings  are  summarized. 
In  Chapter  V,  conclusions  and  recommendations  are  listed.  In  Appendix  A,  the 
sample  problem  unclassified  data  that  J-8  provided  are  presented.  In  Appendix  B, 
the  FORTRAN  code  which  reads  in  the  input  problem  data  and  solves  the  problem 
IS  listed.  In  Appendix  C,  the  output  files  that  the  FORTRAN  program  produced  in 
solving  the  sample  problem  are  listed.  In  Appendix  D,  the  GAMS  (General  Algebraic 
Modeling  System)  input  file  is  presented.  In  Appendix  E,  the  complete  solutions  and 
allocation  summaries  for  each  of  the  parametric  analysis  cases  discussed  in  Chap¬ 
ter  IV  are  listed.  In  Appendix  F,  the  User’s  Guide  which  explains  how  to  use  the 
FORTRAN  program  to  solve  a  missile  allocation  problem  is  presented. 
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II.  Literature  Review 

2.1  Introduction 

The  purpose  of  this  literature  review  is  to  examine  the  currently  available 
MAP  models  and  solution  methodologies  to  see  if  one  of  them  could  be  used  to  solve 
the  J-8  missile  allocation  problem.  This  review  first  discusses  the  terminology  and 
components  that  are  common  to  MAP  models,  then  discusses  the  specifics  of  the 
problem  proposed  by  J-8  in  terms  of  the  components  of  a  MAP,  and  finally  discusses 
solution  techniques  which  may  be  used  to  solve  this  problem. 

2.2  Model  Terminology 

The  following  terms  are  common  throughout  the  literature. 

1.  Damage  Expectancy  (DE):  DE  is  the  expected  amount  of  target  daimage  that 
a  given  weapon  produces  against  a  specific  target.  It  is  based  on  the  probabil¬ 
ities  that  the  weapon  will  launch,  arrive,  detonate,  and  cause  various  levels  of 
damage  (Seiler,  1983:11). 

2.  Hedge:  A  hedge  is  a  constraint  (input  rule)  that  specifies  side  goals  or  limita¬ 
tions  for  allocating  missiles  (Bunnell  and  Takacs,  1984:14). 

3.  Target  Optimization:  Target  optimization  consists  of  allocating  the  available 
missiles  to  the  targets  in  order  to  best  meet  the  main  objective.  The  main 
objective  is  usually  to  maximize  the  expected  amount  of  damage  done  to  the 
set  of  targets  (Bunnell  and  Takacs,  1984:14). 

2.3  Model  Elements 

All  missile  allocation  models  have  five  submodels  in  common:  the  weapon  sys¬ 
tem,  the  target  complex,  the  engagement  model,  the  damage  model,  and  the  solution 
algorithm  (Matlin,  1970:337).  Each  submodel  has  a  certain  level  of  complexity.  In 
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genertJ,  the  more  complex  each  submodel  is,  the  harder  it  is  to  find  the  optimal 
solution.  Thus,  no  model  is  robust  enough  to  handle  each  submodel  at  its  highest 
level  of  complexity  (Matlin,  1970:337). 

2.3.1  The  Weapon  System.  The  weapon  system  consists  of  three  dis¬ 
tinct  categories:  the  scope  (number  of  weapon  types  and  penetration  aids),  weapon 
“reach”  (which  targets  each  weapon  may  strike),  and  weapon  commitment  policy 
(the  number  of  waves  launched,  quality  of  bomb-damage  assessment,  and  weapon 
availability  uncertainties)  (Matlin,  1970:337-339). 

2.3. 1.1  Scope.  The  scope  of  the  weapon  system  is  defined  by  the 
number  of  weapon  types  available  and  whether  or  not  penetration  aids  are  available 
(Matlin,  1970:339).  Penetration  aids  are  those  items  which  help  a  weapon  system 
defeat  the  target  defensive  system.  These  may  include  chaff,  terminal  decoys,  and 
area  decoys  (Bunnell  and  Takacs,  1984:10).  The  simplest  models  would  assume  a 
single  weapon  type  with  no  penetration  aids.  The  most  complex  models  allow  for 
multiple  weapon  types  and  the  use  of  penetration  aids  which  are  usually  modeled  as 
an  overall  probability  of  breaching  tl,  defenses  and  detonating  (Matlin,  1970:339). 
Grotte  (1982:430)  allowed  for  up  to  four  types  of  weapons,  while  the  arsenal  exchange 
model  (AEM)  allows  for  up  to  25  (Bozovich  and  others,  1973:7).  Some  models 
express  weapons  in  equivalent  units  which  permits  easier  solution  of  the  problem 
(Lemus  and  David,  1963:789;  Wambsganss,  1982:10). 

2.3. 1.2  Weapon  Reach.  Weapon  reach  is  usually  given  as  an  inci¬ 
dence  matrix  where  a  “1”  entry  indicates  the  weapon  can  reach  a  given  target  while 
a  “0”  entry  indicates  it  cannot.  The  simplest  models  assume  all  missiles  can  reach  all 
targets.  More  complex  models  consider  range  restrictions.  The  most  complex  mod¬ 
els  consider  various  payloads  which  missiles  may  carry  to  different  ranges  (Matlin, 
1970:339-340,  358-360). 


2.S.1.S  Weapon  Commitment  Policy.  Weapon  commitment  strate¬ 
gies  may  be  divided  into  two  categories:  how  many  waves  are  allowed  in  the  model 
and  whether  the  commitment  is  deterministic  or  probabilistic.  Deterministic  com¬ 
mitment  means  all  missiles  are  available  and  launch  reliably,  and  bonib-damage 
assessment  is  perfect.  Probabilistic  commitment  means  missile  boosters  may  fail  or 
enemy  action  may  eliminate  missiles  before  launch.  The  AEM  allows  for  probabilis¬ 
tic  commitment  (Bozovich  and  others,  1973:IV-B-11).  The  simplest  models,  such 
as  Day’s  targeting  complex  problem  (Day,  1966:992-1013),  assume  a  single  attack 
wave  with  all  missiles  available.  The  most  complex  models  allow  for  more  than 
one  attack  wave  with  each  missile  having  a  probability  of  not  launching  (Matlin, 
1970:340-341).  Grotte’s  model  (Grotte,  1982:433)  went  so  far  as  to  consider  the 
possibility  of  a  catastrophic  failure  for  an  entire  segment  of  the  weapon  inventory. 

2.3.2  The  Target  Complex.  The  target  complex  also  consists  of  three 
distinct  characteristics:  the  types  of  targets  Considered,  the  target  values,  and  the 
defenses  available  to  the  target  complex  (Matlin,  1970:341). 

2.3.2. 1  Types  of  Targets.  Targets  may  be  classified  as  either  point  or 
area,  and  either  dependent  or  independent.  A  point  target  is  one  that  is  small  enough 
that  a  single  weapon  may  destroy  it.  A  missile  silo  is  an  example  of  a  point  target. 
If  the  target  is  large  enough  that  it  requires  more  than  one  weapon  to  destroy  it, 
it  is  an  area  target.  Military  bases  and  cities  are  axea  targets.  A  dependent  target 
is  a  collection  of  point  targets  or  an  area  target  that  has  been  divided  into  several 
individual  aim-points  such  that  a  single  weapon  could  destroy  more  than  one  of 
them.  Otherwise,  point  and  area  targets  are  considered  independent.  The  simplest 
models  consider  only  independent  targets.  The  most  complex  models,  such  as  Day’s 
model  (Day,  1966:993-999),  consider  dependent  area  and  point  targets  (Bunnell  and 
Takacs,  1984:12-13). 


2. 3. 2. 2  Target  Values.  The  usual  measure  of  effectiveness  chosen  is 
the  expected  target  value  killed.  The  simplest  models  consider  (sometimes  implicitly) 
that  all  targets  have  the  same  value.  More  complex  models  rank  targets  in  order 
of  priority  (ordinal  value  scale),  but  do  not  assign  numerical  values.  In  the  next 
level  of  complexity,  models  assign  numericaJ  values  to  each  target  (cardinal  value 
scale).  At  the  highest  level  of  complexity,  models  cillow  for  both  targets  that  have 
indirect  or  intrinsic  value  and  targets  that  have  direct  or  extrinsic  value.  Targets 
with  indirect  (intrinsic)  value  are  those  with  no  value  cissigned  for  killing  them,  but 
if  all  such  targets  are  killed,  then  targets  with  assigned  value  may  be  eliminated 
either  automatically  or  more  easily  (Matlin,  1970:341-342).  Examples  of  targets 
with  indirect  value  would  be  surface-to-air  missile  sites  and  command-and-conlrol 
facilities. 

2.3. 2.3  Available  Target  Defenses.  The  simplest  models  assume  there 
are  no  target  defenses.  More  complex  models  treat  defenses  by  eliminating  some 
of  the  attacking  missiles  before  they  arrive.  This  may  be  done  by  either  having  an 
“admission  price”  (the  defenses  will  shoot  down  a  given  number  of  attacking  missiles) 
or  with  a  probability  that  an  attacking  missile  gets  through  (Bunnell  and  Takacs, 
1984:14). 


2.3.3  The  Engager-ent  Submodel.  The  engagement  submodel  determines 
the  probability  that  a  weapon  brd^ks  through  the  target  defenses  (if  present).  This 
is  based  on  whether  the  offense  is  deterministic  or  probabilistic  and  whether  the 
defense  is  deterministic  or  prcbabil^ic.  Deterministic  offenses  have  missiles  that  hit 
exactly  where  they  are  aimed  and  penetration  aids  that  work  perfectly.  Probabilistic 
offenses  have  missiles  that  may  not  Hit  exactly  where  they  are  aimed  and  penetration 
aids  that  work  imperfectly.  Deterniinistic  defenses  have  radars  that  always  detect 
incoming  missiles  and  anti-missile  missiles  (AMMs)  that  always  work.  Probabilistic 
defenses  have  imperfect  radars  and  AMMs.  The  simplest  models  assume  both  the 
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offense  and  the  defense  are  deterministic.  The  most  complex  models  assume  both 
are  probabilistic  (Matlin,  1970:343-344). 

5.3.4  The  Damage  Submodel.  The  damage  submodel  determines  how 
the  target  damage  or  value  accumulates  as  a  function  of  the  number  and  types  of 
attacking  missiles.  Damage  may  be  either  deterministic  or  probabilistic  and  either 
partial  or  total.  The  simplest  models  assume  a  given  number  of  weapons  will  always 
destroy  a  given  target.  More  complex  models  use  a  probability  that  a  given  target 
will  be  destroyed  by  a  given  number  of  missiles.  A  little  higher  complexity  level 
allows  for  partial  damage  to  accumulate  before  a  target  is  considered  destroyed. 
The  most  complex  models,  like  the  CODE  50  Nuclear  Exchange  Model  (Hillermaa, 
1971:67-71),  allow  for  different  weapon  types,  each  with  different  probabilities  for 
different  levels  of  damage,  and  a  weapon  commitment  strategy  which  permits  attack 
of  targets  by  several  types  of  weapons  (Matlin,  1970:344-345).  Bracken  and  McGill 
use  both  total  and  partial  damage  (Bracken  and  McGill,  1973:31-32). 

2.3.5  The  Solution  Algorithm.  (Matlin,  1970:345-346).  Analysts  have 
used  many  different  algorithms  or  computational  procedures  to  solve  MAPs.  These 
include  analysis  (differentiating  the  expected  damage  equations,  setting  ,ach  to  zero, 
and  solving  the  resulting  system  of  equations),  game  theory,  graphical  techniques, 
graph  theory,  linear  programming,  dynamic  programming,  nonlinear  programming, 
exhaustive  searches,  Monte  Carlo  techniques,  and  combinations  of  these.  Some  key 
properties  to  consider  when  comparing  solution  techniques  are: 

1.  whether  the  solution  is  optimal  or  only  near-optimal,  and  how  optimality  is 
established; 

2.  whether  the  solution  is  integer  or  continuous; 

3.  whether  the  algorithm  also  determines  the  optimal  defensive  allocation;  and 


II 


\ 


4.  the  capability  of  the  algorithm  to  generate  optimal  solutions  and  the  compu¬ 
tational  complexity  of  the  algorithm. 

The  idecil  algorithm  should  provide  integer  solutions,  yield  a  proven 
optimal  solution,  derive  the  optimum  defensive  allocation  as  well  as  the 
optimum  offensive  allocation,  be  capable  of  handling  large  weapon  and 
target  complexes,  run  rapidly,  be  insensitive  to  small  variations  in  weapon 
and  target  numbers  and  associated  parameters,  and  provide  a  global 
rather  than  a  local  or  restricted  solution.  (Matlin,  1970:346) 

8.4  The  J~8  Missile  Allocation  Problem 

J-8  has  defined  certain  input  rules  for  their  specific  missile  allocation  prob¬ 
lem.  There  are  a  given  number  of  fixed  launch  fields  and  a  given  number  of  target 
complexes.  Each  launch  field  contains  a  known  number  of  missiles,  and  each  target 
complex  has  a  known  number  of  targets.  In  terms  of  the  components  of  a  MAP  listed 
above,  the  problem  proposed  by  J-8  may  be  summarized  as  the  following  (Davidson, 
1992b:2-3): 

1.  Weapon  System: 

(a)  Scope:  There  are  two  types  of  missiles:  good  and  fair.  No  penetration 
aids  are  available. 

(b)  Reach:  All  missiles  may  range  all  targets  with  no  degradation. 

(c)  Commitment  Policy:  There  may  be  up  to  two  v.'aves:  1)  an  initial  wave 
which  will,  providing  enough  resources  are  available,  attack  each  target; 
and  2)  a  follow-on  wave  which  will  attack  the  targets  that  were  hit  with 
a  fair  missile  in  the  initial  wave.  The  follow-on  wave  will  be  launched  30 
minutes  after  the  initial  wave  is  launched.  All  missiles  are  available  and 
launch  reliably. 

2.  Target  Complex: 
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(?.)  Types  of  Targets:  All  targets  are  missile  silos  or  command-and-control 
stations  and  are  considered  independent  point  targets. 

(b)  Target  Values:  There  are  three  categories:  good,  fair,  and  poor  (ordinal 
value  scale). 

(c)  Available  Defenses:  None. 

3.  Engagement  Model:  All  missiles  hit  exactly  where  they  are  aimed. 

4.  Damage  Model:  It  takes  one  good  missile  or  two  fair  missiles  to  destroy  any 
target. 

5.  Solution  Algorithm:  The  current  algorithm  employed  by  J-8  is  a  manual  heuris¬ 
tic  which  does  not  necessarily  yield  the  optimal  solution. 

6.  Other  Input  Rules:  Within  this  framework,  J-8  has  imposed  15  allocation  rules 
and  restrictions.  These  rules  and  restrictions  are  addressed  in  Chapter  III. 

2.5  Solution  Techniques 

In  formulating  the  J-8  missile  allocation  problem  as  a  mixed-integer  program 
(see  Chapter  III),  two  solution  techniques,  goal  programming  and  disjunctive  con¬ 
straints,  were  used.  These  techniques  are  discussed  below. 

2.5.1  Goal  Programming.  (Winston,  1991:175-178).  For  some  problems, 
there  may  be  more  than  one  objective  or  goal  that  the  decision  maker  desires  to 
meet.  If  the  decision  maker  can  determine  the  relative  importance  of  meeting  each 
goal,  then  goal  programming  may  be  used  to  solve  the  problem. 

Assume  there  are  n  goals  that  the  decision  maker  seeks  to  achieve.  Each  goal 
may  be  expressed  as  a  function  of  the  decision  variables: 

/i(x)  =  <7.  for  i  =  l,2,...,n 
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where  5,  is  the  decision  maker’s  desired  value  for  /,(x).  In  a  solution  of  this  problem, 
each  gocil  will  either  be  met  exactly,  be  underachieved  by  some  amount  (say  S,),  or 
be  overachieved  by  some  amount  (say  O,).  Then,  in  modeling  the  problem,  each 
goal  becomes  a  constraint  of  this  form: 

/.(x)  +  5i  -  Oi  =  Qi. 

If  we  let  be  the  relative  penalty  of  exceeding  goal  i,  and  W~  be  the  relative 
penalty  for  being  below  goal  t,  then  the  objective  function  may  be  written  as 

min  +  W^-5,). 

1=1 

If  there  are  goals  that  the  decision  maker  desires  to  achieve  or  exceed  (it  does  not 
matter  by  how  much),  set  =  0.  For  those  the  decision  maker  desiies  to  achieve 
but  not  exceed  (it  does  not  matter  how  much  below),  set  Wf  =  0. 

8.5.2  Disjunctive  Constraints.  (Garfirkel  and  Nemhauser,  1972:11).  Con¬ 
sider  a  linear  program  which  has  a  set  of  constraints  (say  m  of  them)  of  which  at 
least  k  constraints  must  hold  (1  <  ^  <  m  —  1).  In  general,  this  type  of  problem 
would  have  the  form: 

mai  (or  min)  {/(x)  :  x  €  S} 

(where  S  is  an  n-dimensional  space  defined  by  the  constraints  of  the  problem)  where 
at  least  k  of  the  following  constraints  must  hold  (1  <  ^  <  m  —  1): 

5i(x)  >  0  for  i  =  1,2, 

To  model  this  situation,  replace  the  above  constraints  with 
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5i(x)  >  SiQi  for  :  = 

m 

Si  <  {or  =)m  —  k 

1=1 

Si  =  Oorl  for  t  =  1,2,  ...,m 

where is  a  known,  negative,  finite  lower  bound  on  gi{'x).  The  first  set  of  constraints 
(  le  for  each  of  the  original  set  of  m  constraints^  states  that  sr,(x)  will  be  greater 
than  or  equal  to  0  if  d,-  is  zero  and  greater  than  or  equal  to  gi  if  Si  is  one.  The  second 
constraint  states  that  the  sum  of  the  SiS  must  be  less  than  or  equal  to  (or  equal  to) 
the  total  number  of  constraints  in  this  set  minus  the  number  of  constraints  of  this  set 
that  must  hold.  Therefore,  at  least  k  of  the  SiS  will  be  set  to  "  ’•o,  which  means  that 
at  least  k  of  the  original  constraints  must  hold.  The  last  constraint  simply  declares 
each  Si  to  be  a  binary  variable. 

2.S  Conclusion 

Many  analysts  have  attempted  to  model  the  missile  allocation  problem.  Some 
considered  only  the  defender’s  objective,  some  only  the  attacker’s  objective,  and  some 
both.  A)  of  the  models  reviewed  consider  the  measure  of  effectiveness  to  be  some 
variation  of  the  expected  value  of  damage  done  to  the  targets.  Several  models  could 
handle  portions  of  the  J-8  problem;  however,  none  can  fully  model  the  J-8  objective: 
to  minimize  the  flight  time  of  the  missile  that  impacts  first  and  to  minimize  the 
duration  of  the  attack. 
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III.  Model  Formulation  and  Methodology 

5.1  Introduction 

This  chapter  contains  the  mixed-integer  programming  formulation  of  the  J-8 
missile  allocation  problem.  First,  the  variables  and  parameters  of  the  model  are 
defined.  Then,  the  mathematical  model  of  the  problem  is  formulated.  Finally,  the 
solution  methodology  is  presented. 

5.2  Variable  Definitions 

In  the  following  discussion,  the  variables  of  the  J-8  missile  allocation  problem 
are  defined.  Let 

LF  =  the  number  of  launch  fields, 

TC  =  the  number  of  target  complexes, 

MT  =  the  number  of  missile  types, 

TT  —  the  number  of  target  types. 

Further,  let 

jo  =  the  set  of  target  complexes  that  contain  at  least  one  good  target, 
jp  =  the  set  of  target  complexes  that  do  not  contain  any  good  targets  but  do 
contain  at  least  one  fair  target, 

jp  =  the  set  of  target  complexes  that  do  not  contain  any  good  or  fair  targets. 

Define 

Xijki  =  the  number  of  missiles  at  launch  field  i  allocated  to  target  complex  j 
which  are  of  type  k  and  attack  type  /  targets 
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for  i  =  1,2,  ...,LF;  j  =  1,2, 


TC-,k  =  1,2, 


MT;  and  /  =  1,2,... ,TT. 


Let 

Thus,  Yijk  is  a  variable  that  indicates  whether  or  not  there  are  missiles  allocated 
from  launch  field  »  to  target  complex  j  of  type  k. 

Let 

So  = 

so  that  Sij  is  a  variable  that  indicates  whether  or  not  there  are  missiles  alloc  ited 
from  launch  field  i  to  target  complex  7. 

Let 


1 


0  ifEf=Tv^o^  =  0 


for  every  i,j 


1  ifES;A:ow>o 


0  i{ZZXijki  =  0 


for  every  i,j,  k. 


Zj  =  the  time  that  target  complex  j  is  attacked.  Zj  is  greater  than  ^  v 
equal  to  the  largest  flight  time  of  missiles  allocated  to  target 
complex  7. 

Tfint  —  the  flight  time  of  the  missile  that  impacts  first 

=  the  earliest  time  of  attack  of  the  target  complexes. 

Tiaat  =  the  flight  time  of  the  missile  that  impacts  last 

=  the  latest  time  of  attack  of  the  target  complexes. 

Dur  =  the  duration  of  the  attack  =  T/ajt  —  27, >,(. 

Oi  =  the  amount  T/.r,*  is  over  its  goal  (which  is  Tmin,  defined  below). 
O2  =  the  amount  Dur  is  over  its  goal  (which  is  0). 


The  parameters  of  the  problem  are: 

Tijk  =  the  flight  time  from  launch  field  t  to  target  complex  7  of  a  type  k  missile. 
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Pet  =  the  fraction  of  that  Dur  is  allowed  to  be. 

rriik  =  the  number  of  missiles  available  at  launch  field  i  of  type  k. 

Tiji  =  the  number  of  targets  at  target  complex  j  of  type  /. 

Tmax  =  the  longest  flight  time  possible  (max{7^jjt})- 
Tmin  =  the  shortest  flight  time  possible  {min{I7ijt}). 

Wi  =  the  weighting  factor  associated  with  Tjint- 

W2  =  the  weighting  factor  associated  with  Dur. 

3.3  Model  Formulation 

3.3.1  Objective  Function.  Since  there  are  two  objectives — mil  imize  the 
flight  time  of  the  missile  that  impacts  first  and  minimize  the  duration  of  the  attack — 
a  goal  programming  approach  is  used.  Since  J-8  desires  the  flight  time  of  the  missile 
that  impacts  first  to  be  as  short  as  possible  and  the  duration  of  the  attack  to  be  as 
short  as  possible,  the  goal  for  Tfint  is  Tmin  and  the  goal  for  Dur  is  0.  There  is  no 
penalty  for  underachieving  (which  is  impossible),  but  there  is  a  penalty  associated 
with  not  achieving  the  goals.  Thus,  the  objective  function  may  be  written  as: 

min  WiOi  +  W2O2 

where  Wi  and  W2  are  user-specified  weighting  factors  associated  with  and 

Dur,  respectively,  which  indicate  the  relative  importance  of  minimizing  each  of  these 
variables  (i.e.,  if  the  decision  maker  wants  to  make  sure  that  T/irat  gets  as  small  2is 
possible  and  feels  that  it  is  not  nearly  as  important  to  minimize  Dur,  ^lssign  a 
value  much  greater  than  W2). 

3.3.2  Constraints.  There  are  several  constraints  on  the  problem.  Each 
type  of  constraint  is  discussed  below. 


1.  Every  target  must  be  hit  by  a  missile  in  the  first  wave.  For  each  target  com¬ 
plex/target  type  combination,  this  can  be  enforced  by  requiring  the  total  num¬ 
ber  of  missiles  allocated  from  each  launch  field  to  a  target  complex  to  be  equal 
to  the  number  of  targets  of  the  given  type  at  the  target  complex.  Mathemati¬ 
cally,  this  becomes 

IF  MT 

«=1  k=l 

2.  The  allocation  of  more  missiles  of  a  given  type  from  a  particular  launch  field 
than  that  launch  field  has  available  is  prohibited.  For  each  launch  field/missile 
type  combination,  this  may  be  enforced  by  requiring  the  total  number  of  mis¬ 
siles  allocated  from  a  launch  field  to  each  target  complex  to  be  less  than  or 
equal  to  the  number  of  missiles  of  the  given  type  available  at  the  launch  field. 
The  appropriate  constraints  are: 

TC  TT 

H  H  ^•Jki  <  m.jfe  for  every  i,  k. 

>=i  (=1 

3.  In  order  to  determine  the  flight  time  of  the  missile  that  impacts  first,  the  small¬ 
est  flight  time  of  the  missiles  that  are  allocated  must  be  identified.  Variable 
Yijk  is  therefore  needed  to  indicate  whether  or  not  missiles  of  type  k  have  been 
allocated  from  launch  field  i  to  target  complex  j.  Two  constraints  are  needed 
for  each  Yijk  variable  to  ensure  that  it  is  assigned  1  when  YTJi  ^ijki  >  0  and 
0  otherwise.  The  first  constraint  should  state  that  for  each  launch  field/target 
complex/missile  type  combination,  the  associated  variable  Yijk  must  be  less 
than  or  equal  to  the  total  number  of  missiles  of  type  k  from  launch  field  i 
allocated  to  target  complex  j.  This  constraint  forces  Yijk  to  be  0  if  there  are 
no  missiles  of  type  k  allocated  from  launch  field  i  to  target  complex  j.  The 
second  constraint  should  state  that  the  associated  variable  Yijk  must  be  greater 
than  or  equal  to  the  total  number  of  missiles  of  type  k  from  launch  field  i  to 
target  complex  j  divided  by  the  total  number  of  targets  at  target  complex  j. 
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This  constraint  forces  Yijk  to  be  1  if  there  are  missiles  of  type  k  allocated  from 
launch  field  i  to  target  complex  j  because  the  number  of  missiles  of  type  k 
alloca  ted  from  launch  field  i  to  target  complex  j  can  never  be  greater  than  the 
total  number  of  missiles  at  t^irget  complex  j.  These  two  sets  of  constraints 
may  be  written  as: 


Yiik  <  EiJiXiSki 


Kifc  > 


1 

I  for  every  t,  j,  k. 

4 


4.  The  desire  that  each  target  complex  be  attacked  by  missiles  from  only  one 
launch  field  generates  a  need  to  determine  which  launch  fields  are  attacking 

each  target  complex.  This  can  be  done  by  setting  S{j  to  1  if  Yijk  >  0 

1 

and  to  0  otherwise.  To  do  this,  two  sets  of  constraints  are  needed.  For  every 
launch  field/target  complex  combination,  the  first  constraunt  should  state  that 
the  associated  variable  5,j  should  be  less  than  or  equal  to  the  sum  of  each 
associated  variable  Yijk-  This  forces  to  be  0  if  there  are  no  missiles  allocated 
from  launch  field  i  to  target  complex  j.  The  second  constraint  should  state  that 
the  associated  5^  must  be  greater  than  or  equal  to  the  sum  of  each  associated 
Yijk  divided  by  the  number  of  missile  types.  This  forces  5,j  to  1  if  there  are 
missiles  allocated  from  launch  field  i  to  target  complex  j  because  the  sum  of 
each  cissociated  Yijk  can  never  be  greater  than  the  number  of  missile  types. 
These  two  constraints  may  be  written  as: 


Sii  <  HSZYtit 


’  for  every  i,j. 

4 


5.  Each  target  complex  should  be  attacked  by  missiles  from  only  one  launch  field. 
For  each  target  complex,  this  may  be  enforced  by  having  a  constraint  that 
states  the  sum  of  the  associated  5yS  must  be  equal  to  1.  This  means  that 


exactly  one  launch  field  attacks  each  target  complex.  This  constrednt  may  be 
expressed  as: 

LF 

Sij  =  1  for  each  j. 

1=1 

If  this  is  infeasible,  then  it  would  be  desirable  for  as  few  target  complexes  as 
possible  to  be  attacked  by  two  launch  fields.  If  the  right  hand  side  of  the  above 
constraint  were  changed  to  be  <  2,  it  would  allow  each  target  complex  to  be 
attacked  by  two  launch  fields.  Thus,  a  constraint  would  be  needed  to  limit  the 
number  of  target  complexes  that  could  be  attacked  by  two  launch  fields.  If 
each  target  complex  were  attacked  by  one  launch  field,  the  sum  of  all  the  Sij 
variables  would  be  equal  to  the  number  of  target  complexes  (TC).  Therefore, 
adding  a  constraint  that  forced  the  sum  of  all  the  5,jS  to  be  less  than  or  equal 
to  the  number  of  target  complexes  plus  one  would  allow  at  most  one  target 
complex  to  be  attacked  by  two  launch  fields.  Thus,  if  it  is  not  feasible  to  attack 
each  target  complex  with  only  one  launch  field,  change  the  right  hand  side  of 
the  above  constraint  to  <  2,  add  the  following  constraint  to  the  formulation: 

LF  TC 

EE^osrc+i, 

i=ij=i 

and  re-solve  the  problem.  This  h€is  the  effect  of  allowing  exactly  one  target 
complex  to  be  attacked  by  two  launch  fields.  If  this  is  also  infeasible,  increase 
the  right  hand  side  of  the  new  constraint  by  one.  Continue  this  until  a  feasible 
solution  is  obtained. 

NOTE:  The  constraint  to  allow  each  target  complex  to  be  attacked  by  missiles 
from  only  one  launch  field  may  also  be  expressed  in  terms  of  only  the 
variables,  which  eliminates  the  need  for  the  Sij  variables.  In  terms  of  just  the 
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Yijk  variables,  this  constraint  becomes 

Yi.j>k')  <  0  for  every  z*,  j*,  k* 

At*  k 

where  z*  is  a  particular  launch  field,  j*  is  a  particular  target  complex,  k*  is 
a  particular  missile  type,  and  M  is  a  large  number.  One  of  these  constraints 
is  needed  for  each  possible  combination  of  launch  field/target  complex/missile 
type.  What  these  constr^unts  do  is  this:  if  is  1,  then  no  other  launch 

field  may  shoot  missiles  of  any  type  at  target  complex  j*.  If  K*;*fc*  is  0,  and  if 
M  is  big  enough,  then  that  particular  constraint  does  not  limit  what  the  other 
YijkS  may  be.  For  M  to  be  “big  enough,”  it  needs  to  be  at  least  as  big  as  the 
total  number  of  possible  combinations  of  launch  fields  and  missile  types,  not 
including  launch  field  z*. 

There  are  advantages  and  disadvantages  to  each  formulation.  If  the  user  de¬ 
cides  to  use  SijS,  the  formulation  is  easier  to  understand,  but  it  requires  more 
binary  variables  which  may  cause  the  solution  time  to  get  very  large.  If  the 
user  decides  to  use  only  VIjfcS,  fewer  binary  variables  are  needed,  so  solution 
times  should  be  fcister.  However,  if  a  particular  set  of  data  happens  to  be  in¬ 
feasible,  this  constraint  is  not  as  easy  to  alter  to  allow  up  to  two  launch  fields 
to  attack  each  target  complex. 

The  number  of  constraints  must  also  be  considered.  If  5,;S  are  used,  the  pro¬ 
gram  requires  an  additional  2  X  LF  X  TC  -\-TC  constraints.  Without  the 
5.is,  the  program  requires  an  additional  LF  x  TC  x  MT  constraints.  For  the 
sample  data,  considering  only  the  four  launch  fields  that  may  be  used  in  the 
allocation,  and  considering  only  feasible  launch  field/missile  type  and  target 
complex/target  type  combinations,  there  is  not  much  of  a  difference  in  the 
total  number  of  constraints  for  either  formulation  (203  with  the  S,;S,  185  with- 
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out).  However,  if  LF,  TC,  or  MT  is  different,  there  could  be  an  appreciable 
difference  in  the  number  of  constraints  required  for  ^tich  formulation. 


6.  During  the  initial  wave,  every  target  in  a  given  complex  must  be  hit  at  the 
same  time.  Thus,  for  a  given  allocation,  the  time  a  given  target  complex  is 
attacked  is  greater  than  or  equal  to  the  flight  time  of  the  missile  that  takes 
the  longest  to  get  there.  All  other  missiles  attacking  that  target  complex  have 
their  launch  times  delayed  so  that  all  missiles  impact  simultaneously.  Since 
Zj  is  the  time  target  complex  j  is  attacked,  Zj  must  be  greater  than  or  equal 
to  the  flight  time  of  each  missile  allocated  to  target  complex  j.  Zj  may  be 
determined  with  this  set  of  constraints: 

Zj  >  TijkYijk  for  every  i,i,  k. 

This  ensures  that  Zj  is  greater  than  or  equal  to  the  flight  time  of  missiles 
allocated  to  target  complex  j  because  if  there  are  no  missiles  of  type  k  allocated 
from  launch  field  t  to  target  complex  j,  then  the  right  hand  side  of  the  above 
constraint  is  zero  because  the  associated  Yijk  is  zero.  Similarly,  if  there  are 
missiles  of  type  k  allocated  from  launch  field  i  to  target  complex  j,  then  Zj  is 
greater  than  or  equal  to  Tijk  because  Yijk  is  set  to  one. 

7.  Since  7/,r,<  is  the  flight  time  of  the  missile  that  impacts  first,  it  should  be  set 
equal  to  the  smallest  Zj.  This  can  be  done  with  a  set  of  disjunctive  constraints 
which  guarantee  that  Tjirst  is  greater  than  or  equal  to  at  least  one  of  the  Zjs: 

Tfir.t  >  Zj  +  Sj{-T^ax)  for  every  j 
TC 

E^;  =  TC- 1 

i=i 

where  6j  is  a  binary  variable  for  each  target  complex^  and  T^ax  is  the  maximum 
flight  time  possible  (the  largest  Tijk).  The  first  set  of  constraints  (one  for  each 
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target  complex)  states  that  T first  is  greater  than  or  equal  to  Zj  if  Sj  is  zero. 
This  is  true  because  T^ax  is  greater  than  or  equal  to  all  other  I^jits,  and  i', 
therefore  greater  them  or  equal  to  each  Zj.  The  last  constraint  states  that 
the  sum  of  the  SjS  must  be  equal  to  the  number  of  tcirget  complexes  minus 
one,  which  means  that  only  one  of  the  first  set  of  constraints  holds  is 

greater  than  or  equal  to  at  least  one  of  the  Zjs).  Since  Zj  is  defined  as  the 
time  that  target  complex  j  is  attacked,  Tfirst  should  be  equal  to  the  smallest 
Zj.  Therefore,  in  order  to  set  T/irst  equal  to  the  smallest  Zj,  another  set  of 
constraints  is  needed  to  ensure  that  T first  is  less  than  or  equcil  to  each  Zj .  The 
appropriate  constraints  are 

Tfirst  <Zj  for  every  j. 

8.  Since  Tust  is  the  flight  time  of  the  missile  that  impacts  last,  it  is  equal  to  the 
time  when  the  last  target  complex  is  attacked,  xhis  means  it  must  be  set 
greater  than  or  equal  to  each  Zj.  This  yields  the  following  constraints: 

Tiast  >  Zj  for  each  j. 

Since  part  of  the  objective  function  is  to  minimize  the  duration,  which  is  the 
difference  between  Tiast  and  Tfirst,  Tiast  is  set  as  small  as  possible.  Because  of 
this  and  the  above  constiaints,  Tiast  is  set  equal  to  the  greatest  Zj. 

9.  Good  targets  must  be  attacked  before  fair  targets,  and  fair  targets  must  be 
attacked  before  poor  ones.  Therefore,  the  time  of  attack  {Zj)  for  each  target 
complex  with  good  targets  must  be  less  than  or  equal  to  the  time  of  attack 
for  each  target  complex  that  does  not  contain  good  targets  but  does  contain 
fair  ones,  and  the  time  of  attack  for  target  complexes  with  fair  targets  and  no 
good  targets  must  be  less  than  or  equal  to  the  time  of  attack  for  each  target 
complex  that  does  not  contain  any  good  or  fair  targets.  Thus,  the  following 
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set  of  constraints  is  needed: 


< 

for  every  j*  €  jo,j' 

e  Jf 

z,. 

< 

Zp 

for  every  j'  6  jF,j^ 

where  jo  is  the  set  of  target  complexes  that  contain  at  least  one  good  target, 
jp  is  the  set  of  target  complexes  that  do  not  contain  any  good  tcirgets  but  do 
contain  fair  ones,  and  jp  is  the  set  of  target  complexes  that  do  not  contain  any 
good  or  fair  targets. 

10.  Constraints  are  also  needed  to  define  the  variables  in  the  objective  function  in 
terms  of  and  Dur.  Since  Oi  and  O2  are  the  amounts  by  which  T/,r,i  and 
Dur  exceed  their  respective  goals,  the  following  constraints  are  needed: 

Dur  =  Tlast  —  Tjirit 
Tjirtt  ~  0\  =  7m<n 

Dur  —  C?2  =  0 

where  Tmin  is  the  minimum  flight  time  possible  (smallest  Tijk)-  The  first  con¬ 
straint  sets  the  duration  [Dur)  to  be  equal  to  the  difference  between  when 
the  last  target  complex  is  attacked  {Tiaat)  and  when  the  first  target  complex  is 
attacked  (T/irat),  the  second  constraint  sets  Oi  to  be  equal  to  how  much 
exceeds  Tmin,  and  the  third  constraint  sets  O2  to  be  equal  to  how  much  Dur 
is  above  zero. 

11.  Since  J-8  also  wants  the  duration  of  the  attack  to  be  no  bigger  than  a  user- 
specified  percentage  of  the  minimum  flight  time,  the  following  constraint  is 
needed: 

Dur  <  Pet  X  Tjirsf 
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This  constraint  states  that  Dur  must  be  less  than  or  equal  to  a  user-specified 

percentage  (Pet)  of 

12.  All  the  input  rules  must  be  addressed  in  the  formulation.  Each  input  rule  is 

presented  here  with  a  discussion  of  how  it  is  addressed  by  the  model. 

(a)  Input  rule  says  there  are  six  launch  fields  and  seven  target  complexes. 
This  is  problem  dependent.  The  model  is  able  to  accommodate  different 
numbers  of  launch  fields  (LF)  and  target  complexes  (TC). 

(b)  Input  rule  says  there  are  two  types  of  missiles.  Agcun,  this  is  prob¬ 
lem  dependent,  and  the  model  handles  different  numbers  of  missile  types 
(MT). 

(c)  Input  rule  #3  says  there  are  three  types  of  targets.  The  model  allows  for 
different  numbers  of  target  types  (TT). 

(d)  Input  rule  #4  says  that  a  wave  coming  from  one  launch  field  is  better 
them  a  wave  launched  from  two  different  launch  fields.  This  is  enforced 
by  constraint  set  5. 

(e)  Input  rule  ^5  says  that  it  takes  one  good  missile  or  two  fair  missiles  to 
destroy  any  target,  and  input  rule  #6  says  that  there  may  be  up  to  two 

- waves:  an  initial  wave  and  a  follow-on  wave  which  must  impact  30  minutes 

after  the  initial  wave  impacts.  This  is  modeled  by  scheduling  the  second 
wave  after  the  initial  solution  is  found.  If  a  target  is  hit  by  a  good  missile 
in  the  initial  wave,  it  is  not  targeted  in  the  follow-on  wave.  If  a  target  is 
hit  by  a  fair  missile  in  the  initial  wave,  it  is  targeted  with  a  fair  missile  in 
the  follow-on  wave. 

(f)  Input  rule  #7  says  that  all  missiles  are  available  for  launch  and  launch 
reliably.  This  is  dealt  with  by  not  decreasing  the  supply  of  available 
missiles. 


26 


(g)  Input  rule  #8  says  that  for  fair  missiles,  it  is  better  to  target  a  given  target 
complex  with  two  waves  from  the  same  launch  field  than  with  one  wave 
from  each  of  two  different  launch  fields.  This  is  handled  by  considering 
only  50  percent  of  the  fair  missiles  as  being  available  for  the  first  wave. 
Thus,  prior  to  solution,  for  m,*  representing  fair  missiles,  m,*  is  set  to  half 
the  actual  number  of  fair  missiles  in  each  launch  field  i  (rounding  down 
for  any  fraction).  This  results  in  the  second  wave  having  an  allocation  of 
fair  missiles  identical  to  that  of  the  first  wave. 

(h)  Input  rule  #9  s^ys  that  each  missile  carries  10  warheads  and  so  may  be 

targeted  against  up  to  10  targets  within  the  same  target  complex.  This  is 

modeled  by  using  the  number  of  missiles  (not  warheads)  at  each  launch 

field  and  by  dividing  the  number  of  targets  (nji)  at  each  target  complex 

1  _ 

by  10  and  rounding  up  for  any  fraction  before  solving  the  problem. 

NOTE:  It  is  possible  that  some  time  in  the  future  there  will  be  different 

i 

missile  types  that  carry  different  numbers  of  warheads.  If  each  type  of 
missile  carries  ai  different  number  of  warheads,  this  may  be  handled  by 
defining  a  new  jparameter,  Numwavk,  which  would  be  the  number  of 

I 

warheads  that  each  type  k  missile  carries.  Then  each  Xijki  would  be  in 
terms  of  warheads  instead  of  missiles,  m,*  would  be  the  total  number  of 
type  k  warheads  at  launch  field  i,  and  Uji  would  be  the  total  number  of 
type  /  targets  at  target  complex  y.  The  problem  may  then  be  solved  in  the 
same  way.  The  resulting  solution  would  be  in  terms  of  warheads  instead 
of  missiles. 

One  difficulty  this  could  cause  is  that  each  Xijki  niay  not  be  an  exact 
multiple  of  the  appropriate  Numivark.  If  this  would  occur,  the  allocation 
would  indicate  that  Xijki  divided  by  Numxcark  (rounded  up  for  fractions) 
missiles  should  be  launched  from  launch  field  i  to  type  I  targets  at  target 
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complex  j.  Any  excess  warheads  should  be  indicated  so  that  the  user  may 
allocate  them  to  other  targets  in  the  vicinity  of  target  complex  j. 

(i)  Input  rule  #10  says  that  good  targets  must  be  attacked  before  fair  targets, 
and  fair  targets  must  be  attacked  before  poor  ones.  This  rule  is  enforced 
by  constraint  set  9. 

(j)  Input  rule  #11  says  that  the  attackers  keep  at  least  10  percent  of  their 

missiles  cis  backup  (strategic  reserve),  and  it  is  best  to  have  10  percent 
of  the  missiles  in  each  launch  field  kept  as  backup.  This  percentage  may 
change  in  the  future.  Thus,  this  number  is  treated  as  a  parameter  input 
by  the  user  (Pctbu).  The  requirement  for  at  least  Pctbu  of  each  launch 
field’s  missiles  to  be  reserved  as  backup  is  handled  by  multiplying  the 
number  of  available  missiles  of  type  k  at  launch  field  i  for  every  i 

and  fc,  by  (1  —  Pctbu)  (rounding  down  for  any  fraction)  before  solving  the 
problem. 

(k)  Input  rule  #12  says  that  two  specific  launch  fields  are  assigned  to  the 
strategic  reserve.  Before  solving  the  problem,  the  user  must  input  whether 
or  not  any  launch  fields  are  for  backup  only.  If  so,  the  user  must  also  input 
how  many  and  which  ones.  When  solving  the  problem,  any  launch  field 
designated  as  backup  only  are  ignored. 


(I)  Input  rule  ,',113  says  that  two  target  complexes  are  sufficiently  close  so  as 
to  be  treated  as  a  single  complex.  Before  solving  the  problem,  the  user 
must  specify  whether  or  not  any  target  complexes  should  be  combined.  If 
so,  the  usei  must  indicate  which  ones.  The  program  then  combines  them 
by  calculating  the  weighted  average  of  the  missile  flight  times: 


Tij-k  = 


TT  , 


for  every  i,  Ji*,  A: 


and  by  summing  the  number  of  targets  at  each  complexr 
Tiju  =  Mji;  +  Tijji  for  every  / 

where  ji  is  one  of  the  target  complexes  to  be  combined,  is  the  other 
one,  and  j*  represents  the  combination  of  the  two  target  complexes.  The 
weighted  average  is  used  because  the  flight  times  are  from  the  centroid  of 
the  launch  field  to  the  centroid  of  the  target  complex. 

(m)  Inpi  c  rule  ^14  says  that  no  “strays”  are  allowed.  This  means  that  the 
allocation  should  be  somewhat  balanced.  This  is  enforced  initially  by 
the  constraints  for  input  rule  ^4  (constraint  set  5),  which  allow  each 
target  complex  to  be  attacked  by  only  one  launch  field.  If  the  problem  is 
infeaisible,  and  those  constraints  are  relaxed  so  that  missiles  from  up  to 
two  launch  fields  may  attack  a  single  target  complex  (see  constraint  set 
5),  then  a  new  solution  is  obtained.  If  this  one  is  not  balanced,  constraints 
may  be  added  to  require  the  solution  to  oe  balanced. 

One  way  to  do  this  is  to  first  identify  which  target  complex  required  two 
laurch  fields  to  attack  all  of  its  targets.  Call  this  complex  ji*.  Then 
re  quire  that  any  launch  field  that  attacks  j*  to  attack  no  more  than  60 
percent  of  the  total  number  of  targets  at  j*.  This  may  be  enforced  by 
adding  a  constraint  for  each  launch  field  that  says  the  total  number  of 
missiles  of  any  type  from  this  launch  field  to  j*  attacking  any  type  of 
target  must  be  less  than  or  equal  to  0.6  times  the  total  number  of  targets 
at  j*.  Mathematically,  this  becomes: 

MT  TT  TT 

EE  ^ii‘ki  ^  0.6  X  for  every  i. 

Jt=i  /=i  /=i 

After  this  set  of  constraints  hcis  been  added,  the  problem  may  be  rc  solved. 
If  there  is  a  feasible  solution  to  this  formulation,  it  will  be  balanced  since 
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only  one  target  complex  may  be  attacked  by  two  launch  fields  and  this 
constraint  forces  to  be  that  one.  However,  the  resulting  solution  may 
not  be  as  good  in  terms  of  minimum  flight  time  or  minimum  duration  as 
the  previous  solution.  This  is  true  because  there  may  be  a  better  feasible, 
balanced  solution  with  a  different  target  complex  being  attacked  by  two 
launch  fields. 

The  only  way  to  know  if  such  a  solution  exists  is  to  force  a  solution  with 
a  different  target  complex  being  attacked  by  two  launch  fields.  One  way 
to  force  this  is  to  remove  the  above  added  constraints  and  add  one  that 
requires  j*  to  be  attacked  by  only  one  launch  field.  This  constraint  may 
be  expressed  in  terms  of  the  Sij  variables: 

LF 

t=l 

This  constraint  forces  target  complex  j*  to  be  attacked  by  only  one  launch 
field,  but  it  does  not  guarantee  anything  else  about  the  solution.  The 
solution  obtained  after  adding  this  constraint  may  be  both  feasible  and 
balanced  and  better  than  the  balanced  solution  where  j*  is  attacked  by 
two  launch  fields,  or  it  may  turn  out  to  be  infeasible,  or  the  solution 
obtained  may  be  worse  (longer  flight  time  or  duration),  or  this  solution 
may  be  unbalanced  as  well.  The  process  of  identifying  the  best  feasible, 
balanced  solution  may  be  time  consuming. 

(n)  Input  rule  #15  says  that  during  each  wave,  every  target  in  a  given  complex 
must  be  hit  at  the  same  time.  This  is  ensured  by  constraint  set  6. 

3.5.5  The  Complete  Model.  The  entire  mixed  integer  formulation  of  the 
problem  is 
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min  WiOi  +  W2O2 


Subject  To: 

Eift  EfclT  ^ijki  =  Tiji  for  every  j,  I 
Ej^i  EUi  Xijki  <  rn,k  for  every  t,  k 

Yijk  <  ElhXiiki  1 

y;,,  >  for  every  i,j,t 

E(,i } 


Sii 

Sii 

<  E^J,Yi:k 

>  for  every  z.  j 
^  1  ^MT  y  1 

—  MT  ^k=l  ^*jk  J 

X^LF  c 

= 

1 

for  every  j 

> 

TiJkYijk 

for  every  i,j,  k 

T/irat 

> 

Zj  +  Sj{—Tmax) 

for  every  j 

z= 

TC-  1 

Tjirat 

< 

for  every  j 

Tu,t 

> 

for  every  j 

Zi> 

< 

2:,. 

for  every  ;*  e  joj'  e  ;> 

Zy 

< 

2,, 

for  every  j'  G  j>, e  ;> 

Variable  Types: 


Dur 

Tjir$t  —  0\ 

Dur  —  O2 
Dur 


—  Tla,t  —  Tfirst 

—  Tmin 
=  0 

<  Pet  X  Tjirst 


^ijki  >  0  for  every  A:,/ 
Yijk  =  0  or  1  for  every  i,j,k 
Sit  =  0  or  1  for  every  i,j 


> 

0 

for  every  j 

Si 

= 

0 

or  1  for  every  j 

Tfirat 

> 

0 

Tlast 

> 

0 

Dur 

> 

0 

Ox 

> 

0 

O2 

> 

0 

Since  Xijki  is  the  number  of  missiles  at  launch  field  i  allocated  to  target  complex 
j  which  are  of  type  k  and  attack  type  /  targets,  each  Xijki  should  be  an  integer 
variable.  Although  Xijki  is  not  required  to  be  integer-valued  in  this  formulation  of 
the  model,  each  assume  integer  value  when  every  Uji  is  an  integer  and  when  each 
target  complex  is  attacked  by  only  one  launch  field  and  only  one  missile  type.  If  this 
is  not  the  case,  Xijki  could  have  non-integer  value  in  the  optimal  solution. 

Since  each  riji  used  in  the  model  is  am  integer,  the  only  way  this  could  happen 
is  if  more  than  one  launch  field  is  attacking  the  target  complex  or  if  a  launch  field 
is  using  more  than  one  missile  type  to  attack  a  target  complex.  If  the  latter  case 
results  in  two  Xijkis  having  non-integer  value,  simply  round  one  of  them  down  cind 
round  the  other  one  up.  This  results  in  a  feasible  solution  since  there  must  be  an 
integer  number  of  missiles  of  each  type  available  at  each  launch  field.  If  there  are 
two  launch  fields  striking  the  same  target  complex  ^lnd  the  associated  Xijkis  are 
non-integer,  then  round  the  larger  Xijki  down  and  the  smaller  one  up.  This  is  again 
feasible  because  there  must  be  an  integer  number  of  missiles  of  each  type  available  at 
each  launch  field.  This  also  helps  ensure  the  allocation  is  balanced  [see  the  handling 
of  input  rule  #14,  section  3.3.2,  item  12(m),  p.  29-30]. 
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If  each  Xijtci  is  required  to  have  integer  value,  the  model  has  an  additional 
LF  X  TC  X  MT  X  TT  integer  variables.  This  could  cause  a  significant  increase  in 
the  solution  time. 

S.4  Solution  Methodology 

The  mathematical  formulation  presented  in  section  3.3.3  is  a  mixed-integer 
program  (MIP)  because  it  has  both  integer  and  continuous  variables.  The  integer 
variables  in  the  above  formulation  are  all  binary.  Thus,  a  MIP  solver  is  needed  to 
solve  this  problem.  However,  the  solver  needs  the  problem  data  entered  into  the 
above  formulation  before  it  can  solve  the  problem.  Since  J-8  desired  a  FORTRAN 
routine  that  they  con'd  incorporate  into  SINBAC,  a  FORTRAN  computer  program 
was  written  to  read  in  the  problem  data,  transform  this  data  for  use  in  the  above 
formulation,  and  then  transmit  the  formulated  problem  to  the  solver.  The  program 
also  takes  the  solution  produced  by  the  solver  and  provides  it  in  an  easily  understood 
format. 

3.4 -i  Procedure.  The  following  is  an  outline  of  a  procedure  which  may  be 
used  to  solve  J-8’s  problem.  Figures  1  and  2  present  a  flow  diagram  of  the  procedure. 

STEP  1:  Read  in  all  the  data. 

a.  Read  in  the  number  of  launch  fields  (LF),  number  of  target  complexes 
(TC),  number  of  missile  types  {MT),  number  of  target  types  {TT),  the  percent  of 
the  minimum  flight  time  that  the  duration  can  be  {Pet),  the  percent  backup  re¬ 
quired  at  each  launch  field  {Petbu),  and  the  relative  importance  of  minimizing  the 
flight  time  and  the  duration  of  the  attack  {W^  and  W2  respectively) 

b.  Read  in  the  flight  time  data  {Tijk),  number  of  missiles  by  type  at  each 
launch  field  {ruik),  and  number  of  targets  by  type  at  each  target  complex  (uj/). 


STEP  2:  Formulate  the  model  presented  in  section  3.3.3  using  the  data  collected 
in  STEP  1. 

STEP  3:  Solve  the  problem  using  the  Zero/One  Optimization  Methods  (ZOOM). 

STEP  4:  Write  the  solution  to  a  file.  If  the  user  is  satisfied  with  the  solution, 
then  go  to  STEP  5.  If  the  user  desires  to  modify  the  solution  directly,  then  perform 
the  modifications  and  go  to  STEP  5.  If  the  user  wants  to  modify  the  problem  and 
re-solve  it,  then  make  the  modifications  and  go  to  STEP  3. 

STEP  5:  Process  and  print  the  final  solution. 


Figure  1.  Flow  Diagram 


IV.  Results  and  Findings 

4.1  Introduction 

Chapter  III  presented  a  mathematical  formulation  of  the  J-8  missile  allocation 
problem  as  a  mixed-integer  program.  In  this  chapter,  the  results  and  findings  of  the 
procedure  presented  in  Chapter  III  are  discussed.  The  presentation  begins  with  a 
more  detailed  description  of  the  procedure. 

To  accomplish  Steps  1  and  2  of  the  procedure  (see  Section  3.4.1),  a  FORTRAN 
computer  program  (see  Appendix  B)  was  written  that  receives  as  input  the  data  for 
the  problem.  This  data  is  then  put  into  the  mixed-integer  formulation  as  specified 
in  Chapter  III  and  the  problem  is  written  to  a  Mathematical  Programming  System 
(MPS)  formatted  file.  At  Step  3,  the  Zero/One  Optimization  Methods  (ZOOM) 
routines  read  the  MPS  file  and  solve  the  missile  allocation  problem.  At  Step  4,  the 
solution  produced  by  ZOOM  is  sent  to  two  output  files:  one  that  is  in  the  format 
J-8  wants  so  that  they  can  review  the  solution,  and  one  that  is  used  to  store  the 
allocation  in  the  format  that  SINBAC  needs. 

If  the  problem  is  infeasible,  the  user  may  modify  either  the  input  data  or  the 
FORTRAN  code.  If  the  user  chooses  to  modify  the  code,  one  or  more  of  the  following 
sets  of  constraints  could  be  relaxed: 

1.  Allow  a  target  complex  to  be  attacked  by  two  launch  fields,  so  change 


to 


add  the  constraint 


IF 

Sij  =  1  for  each  j 

1=1 


LF 

<  2  for  each  j, 
1=1 

LF  TC 
1=1  j=i 
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and  re-solve; 

2.  Do  not  require  that  target?  attacked  by  fair  missiles  in  the  initial  wave  be 
attacked  by  missiles  from  the  same  launch  field  in  the  second  wave.  This 
constraint  may  be  relaxed  by  not  dividing  m,*  for  fair  missiles  in  haJf  prior  to 
solving  (see  the  handling  of  input  rule  #8,  section  3.3.2,  item  12(g),  p.  27); 

3.  Do  not  require  each  launch  field  to  have  a  Pctbu  backup,  so  do  not  multiply 
each  rriik  by  (1  —  Pctbu),  but  instead  add  the  following  constraint: 

LF  TC  MT  TT  LF  MT 

<  (1  -  /’ct6u) 

t=l  i=:l  t=l  /=t  «=1  A:=l 

which  ensures  that  no  more  than  (1  —  Pctbu)  percent  of  the  total  number  of 
available  missiles  are  used  in  the  allocation. 

If  the  solution  is  unacceptable  to  the  user,  either  modify  the  solution  fi!<*s  :• 
rectly,  modify  the  input  data  file  and  re-solve  the  problem,  or  modify  the  FORTRAN 
routines  directly  ^uld  re-solve  the  problem. 

4-2  Solution  to  the  Unclassified  Problem 

The  FORTRAN  program  (with  ZOOM)  was  used  to  solve  the  unclassified 
problem.  Using  10  percent  of  the  earliest  flight  time  for  the  limit  as  to  how  long 
the  duration  could  be  {Pet  =  0.1),  using  weighting  factors  Wi  =  10  and  W2  =  \ 
(tissuming  it  is  much  more  important  to  minimize  the  flight  time  of  the  first  missile 
than  it  is  to  minimize  the  duration  of  the  attack),  and  requiring  each  launch  field  to 
maintain  a  10  percent  backup  {Pctbu  =  0.1),  the  following  solution  was  obtained: 

Objective  Function  Value  =  22.575 

Tfir,t  =  34.027  Tiast  =  37.43  Dur  =  3.403 

O,  =  1.917  O2  =  3.403 
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^1211 

= 

3 

-i^l212 

= 

18 

■^^1221 

= 

15 

A^2321 

= 

21 

■^3422 

= 

12 

■^3623 

=z 

2 

-^4111 

= 

5 

-^4121 

= 

14 

•^4511 

= 

24 

5^121 

— 

1 

1^122 

= 

1 

I232 

1 

5^2 

= 

1 

^362 

z= 

1 

= 

1 

^412 

= 

1 

5^451 

= 

1 

5i2 

= 

1 

S23 

= 

1 

■^34 

=: 

1 

S36 

= 

1 

S4I 

= 

1 

•S’45 

= 

1 

Zi 

= 

35.15 

Z2 

= 

35.98 

^3 

= 

34, 

Z4 

= 

37.43 

Zf, 

= 

37.43 

Ze 

= 

37, 

Si 

= 

1 

S2 

:= 

1 

S4 

= 

1 

Ss 

1 

Se 

=: 

1 

all  other  variables  =  0. 


4-2.1  Allocation  Summary  and  Explanation.  Table  1  summarizes  the  allo¬ 
cation  specified  by  the  solution  to  the  unclassified  problem. 

This  solution  means  that: 

1.  Since  Xuii  =  3  and  A’1212  =  18,  New  Bern  launches  three  good  missiles  at  good 
targets  at  Atlanta/Tybee  and  18  good  missiles  at  fair  targets  at  Atlahta/Tybee. 
Also,  Z2  =  35.98,  so  these  missiles  impact  at  time  35.98,  which  is  the  actual 
flight  time  of  good  missiles  from  New  Bern  to  Atlanta/Tybee.  New  Bern 
started  with  a  total  of  24  good  missiles,  so  it  is  left  with  three  good  missiles  as 
backup. 
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Table  1.  Sample  Problem  Allocation  Summary 


LAUNCH 

FIELD 

MISSILE 

TYPE 

TARGET  COMPLEX 

TARGET  TYPE 

Macon 

Atl/Tyb 

Savannah 

Columbus 

Brunswick 

Athens 

good 

good 

fair 

good 

fair 

good 

poor 

New  Bern 

good 

fair 

3 

15 

18 

Durham 

fair 

21 

Raleigh 

fair 

12 

2 

Wilmington 

good 

fair 

5 

14 

24 

Charlotte 

Lumberton 

2-  ^1221  =  15,  which  means  New  Bern  also  launches  15  fair  missiles  at  good  tar¬ 
gets  at  Atlanta/Tybee.  Since  Z2  —  35.98,  these  missiles  are  launched  in  order 
to  impact  at  35.98,  so  their  launch  time  is  delayed  by  2.74  minutes  (the  actual 
flight  time  of  fair  missiles  from  New  Bern  to  Atlanta/Tybee  is  33.24  minutes). 
This  delay  should  occur  so  that  all  missiles  launched  at  Atlanta/Tybee  impact 
at  the  same  time.  Also,  in  the  second  wave.  New  Bern  should  launch  15  fair 
missiles  at  Atlanta/Tybee.  These  missiles  should  be  launched  at  time  32.74 
(30  minutes  after  the  first  wave  is  launched)  so  that  they  detonate  30  minutes 
after  the  ones  in  the  first  wave.  Since  New  Bern  started  with  a  total  of  34  fair 
missiles,  this  leaves  it  with  four  fair  missiles  as  backup.  - 

3.  Since  A’2321  =  21,  Durham  launches  21  fair  missiles  at  the  21  good  targets  at 
Savannah.  Also,  Z3  =  34.027,  so  these  missiles  should  be  launched  in  order 
to  impact  at  time, 34.027.  This  means  that  their  launch  should  be  delayed  by 
1.577  minutes  (the  actual  flight  time  of  fair  missiles  from  Durham  to  Savannah 
is  32.45  minutes).  This  delay  should  occur  in  order  to  keep  the  duration  of 
the  attack  within  the  given  percentage  {Pci  =  0.1)  of  the  earliest  flight  time. 
During  the  second  wave,  Durham  also  launches  21  fair  missiles  at  Savannah. 
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These  second  wave  missiles  should  be  launched  at  time  31.577.  Durham  started 
with  51  fair  missiles,  so  this  leaves  it  with  nine  fair  missiles  as  backup. 

4.  ^3422  =  12,  so  Raleigh  launches  12  fair  missiles  at  the  12  feur  targets  at  Colum¬ 
bus.  Since  Z\  =  37.43,  these  missiles  should  be  launched  in  order  to  impact  at 
time  37.43,  so  their  launch  must  be  delayed  by  4.01  minutes  (the  actual  flight 
time  of  fciir  missiles  from  Raleigh  to  Columbus  is  33.42  minutes).  This  delay 
should  occur  so  that  all  good  targets  are  attacked  no  later  than  each  fair  tar¬ 
get.  Also,  .^3623  =  2  means  that  Raleigh  launches  two  fair  missiles  at  the  two 
poor  targets  at  Athens.  Since  Zq  =  37.43,  these  missiles  should  be  launched  in 
order  to  impact  at  time  37.43,  so  their  launch  must  be  delayed  by  2.09  minutes 
(the  actual  flight  time  of  fair  missiles  from  Raleigh  to  Athens  is  35.34  minutes). 
Again,  this  delay  should  occur  so  that  good  targets  are  attacked  at  lea^t  as 
early  as  fair  targets.  During  the  second  wave,  Raleigh  also  launches  12  fair 
missiles  at  Columbus  and  two  fair  missiles  at  Athens.  These  missiles  should 
launch  at  times  34.01  and  32.09,  respectively,  in  order  to  impact  30  minutes 
after  the  first  wave  missiles  impact.  Raleigh  started  with  34  fair  missiles,  so 
this  leaves  it  with  six  fair  missiles  as  backup. 

5.  Since  =  5,  Wilmington  should  launch  five  good  missiles  at  Macon.  Zi  = 
35.15,  so  these  missiles  should  be  launched  in  order  to  impact  at  time  35.15, 
which  is  the  actual  flight  time  of  good  missiles  from  Wilmington  to  Macon. 
Also,  A^45h  =  24,  so  Wilmington  launches  24  good  missiles  at  the  24  good' 
targets  at  Brunswick.  Since  Zs  =  37.43,  these  missiles  should  be  launched  in 
order  to  impact  at  time  37.43,  which  is  the  actual  flight  time  of  good  missiles 
from  Wilmington  to  Brunswick.  Wilmington  started  with  33  good  missiles,  so 
this  leaves  it  with  four  good  missiles  ^LS  backup. 

6.  A^412i  =  14  means  that  Wilmington  launches  14  fair  missiles  at  good  targets 
at  Macon.  Since  Zi  =  35.15,  these  missiles  should  also  be  launched  in  order 
to  impact  at  time  35.15,  so  their  launch  must  be  delayed  by  2.85  minutes  (the 
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actual  flight  time  of  fair  missiles  from  Wilmington  to  Macon  is  32.30  minutes). 

.  This  delay  should  occur  so  that  all  missiles  attacking  Macon  impact  at  the  same 
time.  During  the  second  wave,  Wilmington  should  launch  14  additional  fair 
missiles  at  the  Scune  targets  at  Macon.  These  missiles  should  be  launched  at 
time  32.85  in  order  to  impact  30  minutes  after  the  first  wave  does.  Wilmington 
started  with  38  fair  missiles,  so  this  leaves  it  with  ten  fair  missiles  as  backup. 

7.  Since  Charlotte  and  Lumberton  were  designated  as  backup  launch  fields,  they 
do  not  launch  any  missiles,  so  they  have  all  of  their  missiles  as  backup. 

This  solution  is  summarized  by  the  output  files  in  the  formats  that  J-8  re¬ 
quested  (see  Appendix  C).  This  solution  was  obtained  in  57.15  minutes  of  CPU 
(Central  Processing  Unit)  time  on  a  VAX/VMS  6420  computer  system.  Obviously, 
a  larger  problem  may  take  longer  to  solve.  Fortunately,  the  actual  problem  will 
probably  decrease  in  size  in  the  future; 

4-3  Model  Verification  and  Validation 

4.3.1  Comparison  to  GAMS/ZOOM.  In  order  to  verify  that  the  model 
presented  in  Chapter  III — and  implemented  by  the  FORTRAN  program — solved  the 
problem  properly,  this  problem  was  also  input  into  GAMS  (the  General  Algebraic 
Modeling  System).  GAMS  solves  mixed-integer  problems  with  a  modified  version  of 
ZOOM  which  they  call  GAMS/ZOOM.  The  solution  found  by  GAMS/ZOOM  is  as 
follows  (see  Appendix  D  for  the  GAMS  input  file): 

Objective  function  value  — •  22.575 


Tfirst  =  34.027  Tiast  =  37.43  Dur 

O,  =  1.917  Oo  =  3.403 


^  (  / 


■^1211 

= 

3 

Afi2l2 

= 

18 

.^1221 

15 

■^2121 

= 

19 

-^2623 

= 

2 

X3422 

12 

-^4311 

= 

6 

.^4321 

= 

15 

.^4511 

= 

23 

^4521 

= 

1 

5^121 

= 

1 

yi22 

= 

1 

V212 

1 

5^262 

= 

1 

5^342 

= 

1 

J^31 

= 

1 

^432 

.  = 

1 

i^451 

= 

1 

^52 

= 

1 

5i2 

= 

1 

521 

1 

526 

1 

•S'34 

= 

1 

S43 

= 

1 

Sas 

= 

1 

Zx 

34.027 

Z2 

37.43 

Z3 

37, 

z. 

=: 

37.43 

Zs 

= 

37.43 

Ze 

= 

37, 

6% 

= 

1 

S3 

= 

I 

Sa 

1 

s. 

1 

Se 

1 

all  other  variables  =  0 


4-3. 1.1  Allocation  Summary  and  Explanation.  Table  2  summarizes 
the  allocation  specified  by  the  GAMS/ZOOM  solution  to  the  unclassified  problem. 

The  solution  obtained  by  GAMS/ZOOM  is  clearly  different  from  the  one  ob¬ 
tained  by  ZOOM.  In  particular,  flight  times  to  various  complexes  differ.  This  is  the 
result  of  a  different  allocation.  For  example,  Macon  is  now  attacked  by  Durham 
rather  than  Wilmington.  However,  they  are  the  same  in  terms  of  the  objective  func¬ 
tion  value  (to  within  the  ZOOM  and  GAMS  specified  tolerance  of  O.i).  This  means 
that  there  could  be  alternate  optimal  solutions  for  any  given  problem.  There  is  a 
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Table  2.  Allocation  Summary  for  the  GAMS/ZOOM  Solution 


LAUNCH 

FIELD 

MISSILE 

TYPE 

TARGET  COMPLEX 

TARGET  TYPE 

Macon 

Atl/Tyb 

Savannah 

Columbus 

Brunswick 

Athens 

good 

good 

fair 

good 

fair 

good 

poor 

New  Bern 

good 

3 

18 

fair 

15 

Durham 

fair 

19 

2 

Raleigh 

fair 

12 

Wilmington 

good 

6 

23 

fair 

15 

1 

Charlotte 

Lumbertoi. 

slightly  different  missile  allocation  specified  by  each  approach,  but  the  flight  time  of 
the  first  missile  and  the  duration  of  the  attack  (Dur)  are  the  same. 


4-3.8  Face  Validation.  This  model  and  the  results  of  solving  the  unclas¬ 
sified  sample  problem  were  shown  to  J-8  representatives.  Each  of  them  believed 
that  the  solution  obtained  by  the  model  was  a  valid  optimal  solution  to  the  prob¬ 
lem.  Based  upon  expert  opinion,  the  solutions  obtained  using  this  model  are  valid 
solutions  to  the  J-8  missile  allocation  problem. 


4.4  Solution  Using  Number  of  Warheads 

It  is  possible  that  sometime  in  the  future  J-8  could  designate  missiles  for  this 
problem  that  do  not  have  10  warheads  each.  In  fact,  each  missile  type  could  have 
a  unique  number  of  warheads.  Therefore,  this  problem  was  re-solved  with  ZOOM 
using  the  same  unclcissified  data,  but  using  the  number  of  warheads  at  each  launch 
field  and  the  number  of  targets  (not  divided  by  10)  at  each  target  complex.  The 
results  are 


Objective  function  value  =  22.575 


Tjirst 

z= 

34.027 

Tlast 

= 

37.43 

Dur 

= 

3.403 

Ox 

= 

1.917 

O2 

= 

3.403 

-^^1211 

= 

41 

-^1212 

= 

175 

^1221 

= 

134 

■^2121 

= 

185 

-^2623 

= 

20 

■^3422 

= 

120 

= 

210  ■ 

^4511 

= 

87 

■^4521 

= 

153 

5^121 

= 

1 

yi22 

= 

1 

^212 

1 

^262 

1 

5^42 

= 

1 

^431 

= 

1 

^451 

= 

1 

^452 

= 

1 

5i2 

= 

1 

521 

= 

1 

526 

= 

1 

534 

= 

1 

543 

= 

1 

545 

= 

1 

34.027 

^2 

= 

35.98 

^3 

= 

35.60 

^4 

= 

37.43 

^5 

= 

37.43 

= 

37.43 

62 

= 

1 

^3 

= 

1 

^4 

— 

1 

6, 

= 

1 

^6 

== 

1 

all  other  variables  =  0 


4-4-i  Allocation  Summary  and  Explanation.  Taule  3  summarizes  the  al¬ 
location  specified  by  the  solution  to  the  unclassified  problem  when  the  number  of 
warheads  at  each  launch  field  is  used  instead  of  the  number  of  missiles. 

The  measure  of  performance  is  the  same  for  this  problem  as  it  was  for  the 
original  problem.  In  other  words,  this  formulation  yields  the  same  solution  in  terms 
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Table  3.  Allocation  Summary  Using  Warheads  Instead  of  Missiles 


TARGET  COMPLEX 
TARGET  TYPE 


LAUNCH 

MISFi^LE 

Macon 

1  Atl/Tyb  j 

Savannah 

Columbus 

Brunswick 

Athens 

FIELD 

TYPE 

good 

good 

fair 

good 

fair 

good 

poor 

New  Bern 

good 

fair 

41 

134 

175 

Durham 

fair 

185 

20 

Raleigh 

fair 

120 

Wilmngton 

good 

fair 

1 

210 

87 

153 

ChMlotte 

Lumberton 

of  the  flight  time  of  the  missile  that  impacts  first  and  the  duration  of  the  attack. 
However,  the  allocation  is  slightly  different.  | 

Since  Xnn  =  41,  this  would  indicate  that  New  Bern  should  launch  41  good 
warheads  at  Atlanta/Tybee.  However,  since  each  good  mipsile  carries  10  warheads, 
this  would  mean  that  New  Bern  would  be  launching  4.1  good  missiles.  This  is 
impossible,  so  New  Bern  would  really  launch  either  four!  or  five  good  missiles  at 
Atlanta/Tybee.  A  closer  examination  of  the  solution  reveals  that  New  Bern  should 
also  launch  134  fair  warheads  at  Atlanta/Tybee  (since  X1221  =  134),  which  would 
mean  a  launch  of  13  or  14  fair  missiles.  Since  the  number  of  warheads  carried  by  both 
good  and  fair  missiles  is  10,  each  launch  field  started  with  a  multiple  of  10  warheads 
available  for  the  attack  (not  reserved  for  backup).  Therefore,  New  Bern  started  with 
at  least  50  good  warheads  (five  good  missiles)  and  at  least  140  fair  warheads  (14 
fair  missiles)  available  for  the  attack.  This  means  that  in  the  actual  allocation,  New 
Bern  would  either  be  scheduled  to  launch  five  good  and  13  fair  missiles  or  four  good 
and  14  fair  missiles  at  Atlanta/Tybee.  Each  of  these  options  is  feasible,  so  it  would 
be  the  decision  maker’s  choice  as  to  which  allocation  would  be  used. 

A  similar  result  was  obtained  in  the  allocation  of  warheads  from  Wilmington 
to  Brunswick  (X^sn  =  87  and  X4521  =  153).  Here  the  decision  maker  would  have  to 
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choose  between  launching  eight  good  and  16  fair  missiles,  or  nine  good  and  15  fair 
missiles  from  Wilmington  to  Brunswick. 

Sample  output  files  generated  by  the  FORTRAN  program  cire  available  in 
Appendix  C. 


4.5  Parametric  Analysis 

In  order  to  test  the  sensitivity  of  the  unclassified  sample  problem  to  changes 
in  the  input  parameters,  the  parameters  Pet,  Wi,  and  W2  were  alternately  varied 
and  the  problem  re-solved.  The  parameter  Petbu  was  not  varied  since  it  is  currently 
fixed  at  10  percent. 

Nine  different  allocations  were  performed  using  varying  values  for  Pet,  Wi,  and 
Wj.  Table  4  summarizes  the  pertinent  results  of  the  nine  different  solutions. 

Table  4.  Parametric  Analysis  Summary 


CASE 

Pet 

Wi 

W2 

Tfirst 

DUR 

CPU  TIME 

1 

0.1 

10 

1 

34.027 

3.403 

0:57.15 

2 

0.05 

10 

1 

35.648 

1.782 

1:39.13 

3 

0.15 

10 

1 

32.548 

4.882 

0:53.71 

4 

0.1 

2 

1 

34.027 

3.403 

1:39.45 

5 

0.1 

1 

2 

37.43 

0 

1:42.58 

6 

0.05 

2 

1 

35.648 

1.782 

1:37.45 

7 

0.05 

1 

2 

37.43 

0 

1:39.94 

8 

0.1 

1 

1 

34.027 

3.403 

1:46.52 

9 

0.05 

1 

1 

X 

35.648 

1.782 

1:39.08 

Case  1  is  the  original  choices  of  the  input  parameters  with  the  full  results  sum¬ 
marized  in  section  4.2.  All  other  cases  have  their  complete  solutions  and  allocation 
summaries  listed  in  Appendix  E.  The  CPU  times  for  each  of  these  cases  were  ob¬ 
tained  on  a  VAX/VMS  6420  and  are  in  hoursrminutes  (to  the  nearest  one-hundredth 
of  a  minute). 
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Case  2  is  notable  for  its  interesting  results.  Pet  was  set  at  0.05  {Dur  was 
limited  to  five  percent  of  and  it  resulted  in  the  same  value  for  Tiaat  with  T/frjt 

increased  so  that  the  percentage  restriction  on  the  duration  of  the  attack  was  still 
met. 

The  result  from  case  3  seems  to  indicate  that  for  a  given  set  of  problem  data, 
as  the  duration  is  allowed  to  increase  (Pet  gets  larger),  the  flight  time  of  the  first 
missile  decreases  while  the  flight  time  of  the  last  missile  (Tiast)  stays  the 

same. 

The  result  of  case  1  compared  with  the  result  of  case  4,  and  the  result  of  Ccise 
2  compared  with  6,  seems  to  indicate  that  the  optimal  solution  is  the  same  when 
is  greater  than  (i.e.,  minimizing  the  flight  time  of  the  missile  that  impacts  first 
is  more  important  than  minimizing  the  duration  of  the  attack)  for  a  fixed  duration. 
Similarly,  the  results  of  cases  5  and  7  seem  to  indicate  that  if  W2  is  greater  than  Wx 
(it  is  more  important  to  minimize  the  duration),  then  the  value  of  T/,v,t  is  set  to 
equal  Tia,t  so  that  the  duration  is  zero. 
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V.  Conclusions  and  Recommendations 


5. 1  Conclusions 

The  solution  to  the  unclassified  sample  problem  obtained  by  the  model  pre¬ 
sented  in  Chapter  III  was  compared  to  the  solution  obtained  by  GAMS/ZOOM.  Since 
the  measure  of  performance  was  the  same,  the  solution  was  verified.  Also,  expert 
opinion  was  consulted  and  the  solution  obtained  was  considered  valid.  Therefore, 
the  model  presented  in  Chapter  III  may  be  solved  to  derive  optimal  solutions  to 
missile  allocation  problems. 

This  model  is  flexible  because  it  allows  for  various  numbers  of  launch  fields, 
target  complexes,  missile  types,  and  target  types.  This  model  also  allows  the  decision 
maker  to  specify  the  relative  importance  of  each  of  the  two  objectives:  minimize  the 
first  impact  time  and  minimize  the  duration  of  the  attack.  This  flexibility  makes  the 
model  usable  for  a  variety  of  scenarios. 

However,  this  model  would  not  be  practical  to  use  if  the  number  of  binary 
variables  becomes  significantly  large.  The  number  of  binary  variables  increases  by 
{MT  y.  TC)  TC  \{  the  number  of  launch  fields  increases  by  one.  Similarly,  if 
the  number  of  target  complexes  or  missile  types  increases  by  one,  the  number  of 
binary  variables  increases  by  {LF  x  MT)  -|-  LF  -|- 1  or  LF  xTC  respectively.  Since 
the  solution  times  could  potentially  increase  exponentially  as  the  number  of  binary 
variables  increases,  the  model  may  become  impractical  to  use  once  the  number  of 
binary  variables  exceeds  around  100  or  so.  However,  in  the  foreseeable  future,  the 
numbers  of  launch  fields  and  target  complexes  should  decrease,  and  the  number  of 
missile  types  could  drop  to  one,  so  this  model  should  be  usable  for  at  least  several 
years. 

Of  course,  the  future  is  unknown.  Input  rules  could  change;  new  ones  could  be 
added  or  current  ones  deleted.  These  changes  to  the  input  rules  could  require  the 
model  to  be  modified.  It  seems  fairly  simple  to  accommodate  foreseeable  changes  by 

49 


I 


changing,  adding,  or  deleting  constraints  as  appropriate.  Therefore,  even  if  changes 
in  the  world  situation  dictate  significant  changes  to  the  model  as  presented  here,  it 
should  still  serve  as  a  solid  basis  for  building  a  model  which  would  represent  the  new 
situation. 

5.2  Model  Limitations 

This  model  is  a  good  representation  of  the  J-8  missile  allocation  problem,  but 
it  is  not  perfect.  There  is  an  input  rule  that  is  not  modeled  explicitly.  Input  rule 

says  that  no  “strays”  are  allowed.  This  means  that  the  solution  should  be 
somewhat  balanced  (see  the  handling  of  input  rule  #14,  section  3.3.2,  item  12(m), 
p.  29).  If  there  is  a  feasible  solution  that  has  each  target  complex  attacked  by 
only  one  launch  field,  this  requirement  will  be  met.  However,  if  there  is  no  such 
solution,  then  any  target  complex  attacked  by  more  than  one  launch  field  may  have 
an  unbalanced  allocation  (i.e.,  it  may  have  one  launch  field  attacking  90  percent 
of  its  targets  and  another  attacking  10  percent).  There  is  nothing  in  the  model  to 
prevent  this  from  happening.  In  Chapter  III,  one  approach  that  may  be  tried  was 
discussed  (see  section  3.3.2,  item  12(m),  p.  29-30).  However,  this  approach  is  not 
guaranteed  to  yield  the  optimal  solution.  In  fact,  it  is  not  guaranteed  to  yield  even 
a  feasible  solution.  This  is  because  it  is  theoretically  possible  to  have  a  set  of  data 
that  has  no  balanced  solutions  that  are  feasible. 

When  initially  modeling  the  J-8  missile  allocation  problem,  the  constraints  that 
require  the  flight  time  of  the  missile  that  impacts  first  (7/,rj«)  to  be  less  than  or  equal 
to  the  time  of  attack  of  each  target  complex  (Zj)  were  not  included.  It  was  observed 
that  some  solutions  that  were  generated  had  a  first  impact  time  that  was  greater 
than  one  or  more  of  the  target  complex  attack  times  {Tf,rst  >  Zj  for  at  least  one  j). 
This  solution  is  not  consistent  with  the  definition  of  Zj,  which  is  the  time  of  attack 
of  tcirget  complex  j.  Therefore,  the  constraints  mentioned  above  were  added.  This 
had  the  effect  of  producing  solutions  that  are  consistent  with  the  variable  definitions. 
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but  the  solution  times  were  observed  to  increase.  The  solution  times  obtained  for 
the  imclassified  sample  problem  before  these  constraints  were  added  were  in  the 
neighborhood  of  25  to  30  minutes.  After  these  constraints  were  added,  solution  times 
jumped  to  between  55  and  60  minutes.  This  is  «in  increase  by  a  factor  of  two.  For 
this  problem,  this  does  not  cause  the  solution  times  to  become  prohibitive.  However, 
for  a  larger  problem,  including  this  set  of  constraints  may  result  in  the  solution  time 
becoming  too  large  for  the  model  to  be  used  practically.  If  this  is  the  case,  dropping 
these  constraints  may  result  in  the  solution  time  becoming  acceptable.  If  these 
constraints  are  dropped,  the  launch  times  of  the  missiles  allocated  to  impact  prior 
to  T/irat  would  have  to  be  manually  delayed  in  order  to  meet  all  of  the  constraints. 

5.3  Recommendations 

Since  the  input  rule  for  requiring  a  balanced  solution  is  not  explicitly  modeled, 
further  research  should  concentrate  on  finding  a  way  to  model  this  requirement  so 
as  to  guarantee  an  optimal  solution  that  is  both  balanced  and  feasible. 

This  model  generates  optimal  allocations  based  on  the  decision  maker’s  inputs 
for  the  weighting  factors  for  the  two  objectives  {Wi  and  1^2)  and  the  fraction  of 
the  first  impact  time  that  the  duration  is  allowed  to  be  {Pet).  Although  parametric 
analysis  was  conducted  and  the  results  reported  in  Chapter  IV,  further  research  may 
reveal  insights  as  to  potentially  “good”  values  for  these  parameters  based  on  the 
decision  maker’s  desires. 

Since  the  solution  time  for  solving  a  given  problem  is  heavily  dependent  on 
the  number  of  binary  variables,  every  effort  should  be  made  to  eliminate  ones  that 
are  unnecessary.  For  instance,  the  model  calls  for  a  binary  variable  {Yijk)  for  each 
launch  field/target  complex/missile  type  combination.  However,  since  most  of  the 
launch  fields  do  not  have  all  of  the  missile  types,  not  all  the  YijkS  are  necessary.  For 
the  unclassified  sample  problem,  the  basic  model  calls  for  6  x  6  x  2  =  72  different 
Yijk  variables.  However,  only  three  of  the  launch  fields  (New  Bern,  Wilmington, 


and  Chcirlotte)  have  both  missile  types,  so  the  number  of  Yijk  variables  can  be 
reduced  to  (3  x  6  x  2)  +  (3  x  6  x  1)  =  54.  In  addition,  two  of  the  launch  fields  are 
reserved  for  backup  only,  so  the  number  of  necessary  Yijk  variables  may  be  reduced 
to  (2  X  6  X  2)  +  (2  X  6  X  1)  =  36.  If  all  such  reductions  are  made,  the  numbe-  of 
necessary  binary  variables  in  the  unclassified  sample  problem  may  be  reduced  from 


Appendix  A.  Sample  Problem-Unclassified  Data 


A.l  Lz^unch  Fields  and  Missile  Data 

There  are  six  launch  fields:  New  Bern,  Durham,  Raleigh,  Wilmington,  Char¬ 
lotte,  and  Lumberton.  Each  launch  field  has  one  or  both  missile  types  as  follows: 


Table  5.  Launch  Fields  and  Missile  Data 


LAUNCH  FIELD 

TYPE 

MISSILES 

WARHEADS 

New  Bern 

good 

24 

i 

fair 

34 

mESm 

;  Durham 

fair 

51 

510 

Raleigh 

fair 

34 

340 

Wilmington 

good 

33 

fair 

38 

1  Charlotte 

good 

24 

1 

fair 

25 

Lumberton 

fair 

58 

580 

A. 2  Target  (^omplexes  and  Target  Data 

There  Are  seven  target  complexes:  Macon,  Atlanta,  Savannah,  Columbus, 
Brunswick,  Tybee,  and  Athens.  Each  target  complex  has  either  good,  fair,  or  poor 
targets.  Table  6  contains  this  data. 


Table  6.  Target  Complexes  and  Target  Data 


TARGET  COMPLEX 

TYPE 

NUMBER 

Macon 

good 

.  185 

Atlanta 

good 

175 

Savannah 

good 

210 

Columbus 

fair 

120 

Brunswick 

|||||^^]| 

240 

Tybee 

fair 

175 

Athens 

poor 

20 

A.S  Flight  Time  Data 


Table  7  lists  the  flight  time  from  each  launch  field  to  each  target  complex  by 
missile  type. 

Table  7.  Flight  Time  Data 


LAUNCH 

MISSILE 

TARGET  COMPLEX 

FIELD 

TYPE 

Macon 

Atlanta 

Savannah 

Columbus 

Brunswick 

Tybee 

Athens 

New  Bern 

good 

34.16 

34.81 

33.89 

35.33 

36.49 

37.15 

38.22 

fair 

32.61 

32.76 

32.66 

33.10 

33.61 

33.72 

34.82 

Durham 

fair 

32.11 

32.03 

32.45 

32.51 

32.96 

32.78 

34.71 

Raleigh 

fair 

32.95 

33.02 

33.09 

33.42 

33.92 

33.93 

35.34 

Wilmington 

good 

35.15 

35.26 

35.60 

36.26 

37.43 

34.42 

40.78 

fair 

32.30 

32.20 

32.67 

32.70 

33.14 

32.93 

34.94 

Charlotte 

good 

35.45 

35.66 

35.76 

36.58 

37.76 

37.87 

37.81 

fair 

32.58 

32.52 

32.90 

33.00 

33.46 

33.30 

35.19 

Lumberton 

fair 

33.41 

33.44 

33.60 

33.88 

34.38 

34.33 

35.91 
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Appendix  B.  FORTRAN  Code 

B.l  Introduction 

The  program  given  in  Appendix  B  was  written  in  its  entirety  by  the  author  of 
this  thesis.  This  program  calls  ZOMAIN,  which  is  the  top  level  routine  for  ZOOM,  as 
a  subroutine.  ZOMAIN  and  all  other  subroutines  contained  in  ZOOM  were  written 
by  Dr.  Roy  E.  Marsten,  currently  of  the  Georgia  Institute  of  Technology. 

In  order  to  get  ZOOM  to  interact  with  the  FORTRAN  program  that  this  author 
wrote,  and  to  get  it  to  run  on  the  VAX/VMS  system  at  the  Air  Force  Institute  of 
Technology,  the  following  changes  were  made  to  the  ZOOM  code: 

1.  ZOMAIN  was  changed  from  a  PROGRAM  to  a  SUBROUTINE. 

2.  All  references  to  the  subroutine  TIMER  were  commented  out. 

3.  In  the  subroutine  ZDRIVE,  a  variable  was  initialized  to  one  instead  of  zero 
because  that  value  is  used  as  an  array  dimension. 

4.  The  parameter  MAXVAR  and  the  array  VARVAL(MAXVAR)  are  passed  as 
arguments  into  ZOMAIN,  then  ZDRIVE,  and  then  MPSOUT  in  order  to  collect 
the  values  of  the  solution  variables  to  psiss  back  to  the  main  program. 

5.  Statements  were  added  to  various  subroutines  to  print  error  messages  to  the 
file  SILOATT.ERR  whenever  the  program  encounters  difficulties. 

B.2  Limitations 

The  FORTRAN  program  written  for  this  thesis  implements  this  model  and 
solves  missile  allocation  problems.  However,  it  has  some  limitations: 

1.  The  subroutine  that  reads  in  the  data  is  “fragile”  in  that  if  the  user  inputs  the 
wrong  type  of  data  (i.e.,  entering  integer  data  where  real  data  is  required),  the 
program  terminates.  Future  research  should  include  modifying  this  subroutine 
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so  that  if  the  user  accidentally  enters  the  wrong  type  of  data,  the  program  does 
not  terminate  but  instead  gives  an  error  message  as  to  what  the  problem  is. 

2.  If  the  user  Wcints  to  modify  a  constraint,  the  FORTRAN  code  must  be  modified 
in  each  part  of  the  program  that  applies  to  that  constraint,  particularly  in  the 
formulation  subroutine  (FORMUL)  which  formulates  the  data  into  the  appro¬ 
priate  mixed-integer  program  and  writes  it  to  an  MPS  file.  Future  research 
should  concentrate  on  making  it  easier  to  modify  constraints.  This  could  in¬ 
clude  having  a  main  menu  that  would  allow  the  user  to  choose  inputting  a  new 
problem,  modifying  an  existing  input  data  file,  modifying  the  current  solution 
files  directly,  or  modifying  constraints  on  the  problem. 


B.3  FORTRAN  Code 

PROGRAM  SILOATT 

C***it^******’¥**********t**********mt************ii*****************t*******^**m* 

RAME;  SILOATT 
C 

C  REV  DATE  OF 

C  HO.  CHARGE  DESIGNER  CCR  DESCRIPTION 

C  -  -  - -  - 

C  000  9FEB93  R.  PACE  ORIGINAL  RELEASE 

C 

C  DESCRIPTIOH:  DETERMINES  THE  OPTIMAL  MISSILE  ALLOCATION  FOR  THE  SILO  ATTACK 

C  _ 

C  CALLED  BY:  RONE 
C 

C  CALLS: 

C 

C  IRPUT  -  COLLECTS  THE  DATA  FOR  THE  SILO  ATTACK  PLAH 
C  PROCESS  -  PROCESSES  THE  DATA  FOR  THE  SILO  ATTACK  PLAH 
C  FORMUL  -  WRITES  THE  DATA  INTO  AN  MPS-FORMAT  FILE 

C  ZOMAIN  -  INVOKES  ZOOM  TO  SOLVE  THE  ALLOCATION 

C  PRHTSOL  -  PRINTS  THE  SOLUTION  INTO  A  FILE 
C 

C  CALLING  SEQUENCE:  NONE 
C 

C  FILES:  RONE 
C 

C  LOCAL  VARIABLES: 

C 

C  BU  -  ARRAY  STORING  ID  OF  BACKUP  LAUNCH  FIELDS 


/' 


C  LF  -  THE  HUMBER  OF  LAUKCH  FIELDS 

C  LFMT  -  THE  HUMBER  OF  LAUHCH  FIELD/MIS  .ILE  TYPE  COMBIHATIOHS  THAT  EXIST 

CM  -  ARRAY  STORIHG  THE  HUMBER  OF  MISSILES  OF  EACH  TYPE  AT  EACH  LAUHCH 

C  FIELD 

C  MADJ  -  ARRAY  STORIHG  THE  ADJUSTED  HUMBER  OF  MISSILES  OF  EACH  TYPE  AT 
C  EACH  LAUHCH  FIELD 

C  MT  -  THE  HUMBER  OF  MISSILE  TYPES 

C  H  -  ARRAY  STORIHG  THE  HUMBER  OF  TARGETS  OF  EACH  TYPE  AT  EACH  TARGET 

C  COMPLEX 

C  HADJ  -  ARRAY  STORIHG  THE  ADJUSTED  HUMBER  OF  TARGETS  OF  EACH  TYPE  AT  EACH 
C  TARGET  COMPLEX 

C  HUMBU  -  THE  HUMBER  OF  BACKUP  LAUHCH  FIELDS 

C  HUMFAIR  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  COHTAIH  HO  GOOD  TARGETS  BUT 

C  DO  COHTAIH  FAIR  TARGETS 

C  HUMGOOD  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  COHTAIH  GOOD  TARGETS 

C  HUMPOOR  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  COHTAIH  HEITHER  GOOD  HOR 

C  FAIR  TARGETS  BUT  DO  COHTAIH  POOR  TARGETS 

C  HUMTAR  -  ARRAY  STORIHG  THE  TOTAL  HUMBER  OF  TARGETS  AT  EACH  TARGET  COMPLEX 
C  PCT  -  PERCEHT  OF  THE  MIHIMUM  FLIGHT  TIME  THAT  THE  DURATIOH  IS  ALLOWED 

C  TO  BE 

C  PCTBU  -  THE  PERCEHT  OF  THE  AVAILABLE  MISSILES  AT  EACH  LAUHCH  FIELD 
C  THAT  MUST  BE  RESERVED  FOR  BACKUP 

C  T  -  ARRAY  STORIHG  THE  FLIGHT  TIME  OF  A  TYPE  K  MISSILE  FROM  LAUHCH 

C  FIELD  I  TO  TARGET  COMPLEX  J 

C  TC  -  THE  HUMBER  OF  TARGET  COMPLEXES 

C  TCTT  -  THE  HUMBER  OF  TARGET  COMPLEX/TARGET  TYPE  COMBIHATIOHS  THAT  EXIST 

C  TMAX  -  THE  LARGEST  POSSIBLE  FLIGHT  TIME 

C  TMIH  -  THE  SMALLEST  POSSIBLE  FLIGHT  TIME 

C  TT  -  THE  HUMBER  OF  TARGET  TYPES 

C  VARVAL  -  ARRAY  THAT  STORES  THE  SOLUTIOH  VALUES 

C  W1  -  THE  RELATIVE  WEIGHT  OF  MIHIMIZIHG  THE  EARLIEST  FLIGHT  TIME 

C  W2  -  THE  RELATIVE  WEIGHT  OF  MIHIMIZIHG  THE  DURATIOH  OF  THE  ATTACK 

C 

C  ERROR  MESSAGES:  HOME 
C 

C***0****m***********t*************<^*l*^:^^^* 

C 

C  DIMEHSIOH  PARAMETERS: 

C 

C  MAXLF  =  MAX  HUMBER  OF  LAUHCH  FIELDS 

C  MAXMT  =  MAX  HUMBER  OF  MISSILE  TYPES 

C  MAXTC  =  MAX  HUMBER  OF  TARGET  COMPLEXES 

C  MAXTT  =  MAX  HUMBER  OF  TARGET  TYPES 

C  MAXVAR  =  MAX  HUMBER  OF  VARIABLES 

C 

IHTEGER  MAXLF,  MAXMT,  MAXTC,  MAXTT,  MAXVAR 

PARAMETER  (  MAXLF  =  10,  MAXMT  =  S,  MAXTC  =  10,  MAXTT  =  6, 
ft  MAXVAR  =  600  ) 

C  LOCAL  VARIABLE  DECLARATIOHS: 
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IKTEGER  BU(MAXLF) ,  LF,  LFMT*  M(MAXLF,MAXMT) ,  MADJ(MAXLF,MAXMT) , 
ft  MT,  XMAXTC.MAXTT),  HADJ(MAXTC.MAXTT) ,  BUMBU,  NUMFAIR, 

ft  lURGOOD,  lUMPOOR,  >UMTAR(MAXTC) ,  TC,  TCTT,  TT 

REAL  PCT,  PCTBU,  T(MAXLF,MAXTC,MAXMT) .  TMAX,  TMI8,  Wl,  W2 

DOUBLE  PRECISION  VARVAL(NAXVAR) 

CHARACTER'S  CPLX(MAXTC).  FIELD (MAXL7) ,  MANE 

CHARACTERS  MTYPE(HAXMT) ,  TTYFE(MAXTT) 

C  INPUT  COLLECTS  THE  DATA  FOR  THE  SILO  ATTACK  PLAN: 

CALL  INPUT(CPLX,  FIELD,  LF.  MAXLF,  MAXHT,  MAXTC,  MAXTT,  MT, 
ft  HTYPE,  N.  PCT.  PCTBU,  T.  TC,  TT,  TTYPE,  HI,  H2) 

C  PROCESS  PROCESSES  THE  DATA  FOR  THE  SILO  ATTACK  PLAN: 

CALL  PROCESSCBU,  CPLX,  LF,  LFMT,  M,  MADJ,  MAXLF,  MAXMT,  MAXTC. 
ft  MAXTT,  MT,  H,  NADJ,  lUMBU,  NUMFAIR,  NUMGOOD,  NUMPOOR, 

ft  NUMTAR,  PCTBU,  T,  TC,  TCTT,  TMAX,  TMIN,  TT) 

C  FORMUL  WRITES  THE  DATA  INTO  AH  MPS-FORMAT  FILE: 

CALL  FORHULCBU,  LF,  LFMT,  MADJ.  MAXLF,  MAXMT,  MAXTC,  MAXTT,  MT, 
ft  NADJ,  NUMBU,  NUMFAIR,  NUMGOOD.  NUMPOOR,  NUMTAR,  PCT, 

ft  T,  TC,  TCTT,  TMAX,  TMIN,  TT,  HI,  H2) 

C  ZOMAIN  INVOKES  ZOOM  TO  SOLVE  THE  ALLOCATION; 

CALL  ZOHAIN(HAXVAR,  VARVAL) 

C  PRNTSOL  PRINTS  THE  SOLUTION  INTO  FILES: 

CALL  PRNTSOL(BU,  CPLX,  FIELD,  LF,  MADJ,  MAXLF,  MAXMT,  MAXTC, 
ft  MAXTT,  MAXVAR,  MT,  MTYPE,  NADJ,  NUMBU,  T,  TC,  TT, 

ft  TTYPE,  VARVAL) 

END 

C***********«******4>***********************************************4>*******«**** 

SUBROUTINE  INPUT(CPLX,  FIELD,  LF,  M,  MAXLF,  MAXMT,  MAXTC,  MAXTT, 
ft  MT,  MTYPE,  N,  PCT.  PCTBU.  T,  TC,  TT,  TTYPE,  HI, 

ft  H2) 

c 


NAME: 

INPUT 

REV 

DATE  OF 

NO. 

CHANGE 

DESIGNER 

CCR  DESCRIPTION 

000 

SFEBOS 

R.  PACE 

ORIGINAL  RELEASE 

c 
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C  DESCRIPTION:  COLLECTS  THE  DATA  FOR  THE  SILO  ATTACK  PLAN 
C 

C  CALLED  BY:  SILOATT 
C 

C  CALLS:  HONE 

C 

C  CALLING  SEQUENCE: 

C 

CIN  MAXLF  -  MAX  NUMBER  OF  LAUNCH  FIELDS 

C  MAXMT  -  MAX  HUMBER  OF  MISSILE  TYPES 

C  MAXTC  -  MAX  HUMBER  OF  TARGET  COMPLEXES 

C  MAXTT  -  MAX  NUMBER  OF  TARGET  TYPES 

COUT  CPLX  -  ARRAY  STORING  THE  NAMES  OF  THE  TARGET  COMPLEXES 

C  FIELD  -  ARRAY  STORING  THE  NAMES  OF  THE  LAUNCH  FIELDS 

C  LF  -  HUMBER  OF  LAUNCH  FIELDS 

CM-  ARRAY  STORING  NUMBER  OF  MISSILES  BY  TYPE  AT  EACH  LAUNCH  FIELD 
C  MT  -  HUMBER  OF  MISSILE  TYPES 

C  MTYPE  -  ARRAY  STORING  THE  NAMES  OF  THE  MISSILE  TYPES 

C  N  -  ARRAY  STORING  NUMBER  OF  TARGETS  BY  TYPE  AT  EACH  TARGET  COMPLEX 

C  PCT  -  PERCENT  OF  TFIRST  THAT  DUR  IS  ALLOWED  TO  BE 

C  PCTBU  -  THE  PERCENT  OF  THE  AVAILABLE  MISSILES  AT  EACH  LAUNCH  FIELD 

C  THAT  MUST  BE  RESERVED  FOR  BACKUP 

C  T  -  ARRAY  STORING  FLIGHT  TIMES  OF  MISSILES  (BY  TYPE)  FROM  EACH  LAUNCH 

C  FIELD  TO  EACH  TARGET  COMPLEX 

C  TC  -  NUMBER  OF  TARGET  COMPLEXES 

C  TT  -  NUMBER  OF  TARGET  TYPES 

C  TTYPE  -  ARRAY  STORING  THE  NAMES  OF  THE  TARGET  TYPES 
C  W1  -  WEIGHTING  FACTOR  FOR  TFIRST 

C  W2  -  WEIGHTING  FACTOR  FOR  DUR 

C 

C  FILES : 

C 

C  SILOATT.DAT 

C  SILOATT. ERR 

C 

C  LOCAL  VARIABLES: 

C 

C  I,  J,  K,  KK,  L,  LL,  q  -  COUNTERS 

C  TERROR  -  STORES  THE  VALUE  OF  lOSTAT 

C  LINE  -  STORES  THE  CONTENTS  OF  A  LINE  FRON  'SILOATT.DAT' 

C  MISS  -  TEMPORARILY  STORES  THE  NAME  OF  A  MISSILE  TYPE 

C  NAME  -  TEMPORARILY  STORES  THE  NAME  OF  A  LAUNCH  FIELD  OR  TARGET  COMPLEX 

C  NUM  -  TEMPORARILY  STORES  THE  NUMBER  OF  MISSILES  OR  TARGETS 

C  TARG  -  TEMPORARILY  STORES  THE  NAME  OF  A  TARGET  TYPE 

C  TEMP  -  ARRAY  THAT  TEMPORARILY  STORES  FLIGHT  TIME  DATA 
C 

C  ERROR  MESSAGES:  NONE 
C 

C******^********************'*>****1'***************4'***#*****4>**<^*4'*4>*>)'^********** 

c 

C  LOCAL  VARIABLE  DECLARATIONS: 
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INTEGER  I,  lERROR,  J,  K,  KK,  L,  LF,  LL,  M(MAXLF.MAXMT) ,  MT, 
ft  KMAXTC.MAXTT),  HUM.  Q,  TC,  TT 

REAL  PCT,  PCTBU,  T(MAXLF,HAXTC,MAXMT) ,  TEMP(IO),  Wl,  W2 

CHARACTER^SO  LINE 

CHARACTER*8  CPLX(MAXTC),  FIELD(MAXLF) .  NAME 

CHARACTER*4  MISS,  TARG,  MTYPE(MAXMT),  TTYPE(MAXTT) 

1  =  0 
J  =  0 
K  =  0 
L  =  0 
LF  =  0 
MT  =  0 
TC  =  0 
TT  =  0 

0PEN(UNIT=10,  FILE=' siloatt.dat',  ACCESS= 'SEQUENTIAL ' , 
ft  STATUS='OLD',  IOSTAT=IERROR,  ERR=9991) 

10  READ(10,1000,EKD=EOO)  LINE 

C  IF  THE  FIRST  CHARACTER  OF  A  LINE  IS  BLANK,  READ  THE  NEXT  LINE. 

IF  (LIHE(1:1)  .EQ.  '  ')  GO  TO  10 

C  IF  THE  FIRST  12  CHARACTERS  OF  THE  LINE  ARE  'MISSILE  DATA’,  THEN  THE 
C  FOLLOWING  LINES  GIVE  THE  NUMBER  OF  MISSILES  BY  TYPE  AT  EACH  LAUNCH  FIELD 
IF  (LIHE(1;12)  .EQ.  'MISSILE  DATA’)  THEN 

C  SKIP  THE  NEXT  SIX  LINES: 

DO  20  q  =  1,  6 
READdO,*) 

20  CONTINUE 

26  READ (10, 1010, END=SOO)  NAME,  MISS,  NUM 

C  IF  THE  FIRST  CHARACTER  OF  A  LINE  IS  BLANK,  THIS  SECTION  IS  DONE. 

IF  (NAME(1:1)  .EQ.  ’  ’)  GO  TO  10 

C  OTHERWISE,  READ  IN  THE  DATA.  IF  NO  LAUNCH  FIELDS  HAVE  BEEN  DETERMINED, 

C  MAKE  THIS  ONE  THE  FIRST  ONE. 

IF  (I  .EQ.  0)  THEN 
I  =  1 

FIELD (I)  =  NAME 
LF  =  LF  +  1 

C  ELSE  IF  THIS  IS  THE  SAME  LAUNCH  FIELD  AS  BEFORE.  CHECK  THE  MISSILE  TYPE. 
ELSE  IF  (NAME  .EQ.  FIELD(I))  THEN 
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GO  TO  30 


C  ELSE,  THIS  IS  A  HEW  LAUNCH  FIELD,  SO  ADD  ITS  HAME  TO  THE  LIST: 

ELSE 

1  =  1  +  1 
FIELD(I)  =  HAME 
LF  =  LF  +  1 
EHD  IF 

C  CHECK  THE  MISSILE  TYPE.  IF  THIS  IS  THE  FIRST  OHE,  ADD  IT  TO  THE  LIST: 

30  IF  (K  .EQ.  0)  THEH 
K  =  1 

MTYPE(K)  =  MISS 
MT  =  MT  +  1 

C  ELSE,  CHECK  THE  LIST.  IF  THIS  TYPE  IS  ALREADY  IK  THE  LIST,  ADD  THE  HUMBER 
C  OF  MISSILES  TO  THE  APPROPRIATE  LAUHCH  FIELD. 

ELSE 

DO  40  KK  =  1,  HT 

IF  (MISS  .EQ.  MTYPE(KK))  THEH 
K  =  KK 
GO  TO  60 
EHD  IF 

40  COHTIIUE 

K  =  MT  +  1 
MTYPE(K)  =  HISS 
MT  =  MT  +  1 
EHD  IF 

SO  M(I,K)  =  HUH 

GO  TO  25 

C  IF  THE  FIRST  11  CHARACTERS  OF  THE  LIKE  ARE  'TARGET  DATA’,  THIS  SECTIOH 
C  COHTAINS  THE  HUMBER  OF  TARGETS  BY  TYPE  AT  EACH  TARGET  COMPLEX: 

ELSE  IF  (LIHE( 1:11)  .EQ.  'TARGET  DATA')  THEH 

C  SKIP  THE  HEXT  FOUR  LINES; 

DO  60  Q  =  1,  4 
READClO,*) 

60  CONTINUE 

66  READClO, 1010, EHD=500)  NAME,  TARG,  HUM 

C  IF  THE  FIRST  CHARACTER  OF  A  LINE  IS  BLAHK,  THIS  SECTIOH  IS  DONE. 

IF  (HAME(1:1)  .EQ.  ’  ')  GO  TO  10 

C  OTHERWISE,  IF  THIS  IS  THE  FIRST  TARGET  COMPLEX,  ADD  ITS  NAME  TO  THE  LIST: 
IF  (J  .EQ.  0)  THEN 
J  =  1 

CPLX(J)  =  NAME 
TC  =  TC  +  1 
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C  ELSE  IF  THIS  COMPLEX  IS  ALREADY  IH  THE  LIST,  CHECK  THE  TARGET  TYPE. 
ELSE  IF  (HAME  .EQ.  CPLX(J))  THEH 
CO  TO  70 

C  ELSE,  THIS  IS  A  HEW  TARGET  COMPLEX,  SO  ADD  IT  TO  THE  LIST: 

ELSE 

J  =  J  +  1 
CPLX(J)  =  lAME 
TC  =  TC  +  1 
EHD  IF 

C  CHECK  THE  TARGET  TYPE.  IF  THIS  IS  THE  FIRST  OHE,  ADD  IT  TO  THE  LIST. 
70  IF  (L  .EQ.  0)  THEH 
L  =  1 

TTYPE(L)  =  TARG 
TT  =  TT  +  1 

C  OTHERWISE,  CHECK  THE  LIST  TO  SEE  IF  THIS  OHE  IS  ALREADY  OH  IT. 

ELSE 

DO  80  LL  =  1,  TT 

IF  (TARG  .EQ.  TTYPE(LL))  THEH 
L  =  LL 
GO  TO  90 
EHD  IF 

80  COHTIHUE 

L  =  TT  +  1 
TTYPE(L)  =  TARG 
TT  =  TT  +  1 
EHD  IF 

90  H(J,L)  =  HUM 
GO  TO  65 


C  IF  THE  FIRST  16  CHARACTERS  OF  THE  LIKE  ARE  'FLIGHT  TIME  DATA',  THE  HEXT 
C  SECTION  LISTS  THE  FLIGHT  TIMES  OF  MISSILES  (BY  TYPE)  FROM  EACH  LAUHCH 
C  FIELD  TO  EACH  TARGET  COMPLEX. 

ELSE  IF  (LIHE(1:16)  .EQ.  'FLIGHT  TIME  DATA')  THEH 
DO  lOO  Q  =  1,  6 
READdO,*) 

100  COHTIHUE 

110  RiEAD(lO,1020,EHD=500)  HAME,  MISS,  (TEMP(J),  J  =  1,  TC) 


C  IF  THE  FIRST  CHARACTER  OF  A  LIHE  IS  BLAHX,  THIS  SECTIOH  IS  DOME: 
IF  (IAME(1:1)  .EQ.  '  ')  GO  TO  10 


C  FIHD  THE  CURREHT  HAME  IH  THE  LIST  OF  LAUHCH  FIELDS: 
D0)l20  Q  =  1,  LF 

IF  (HAME  .EQ.  FIELD(Q))  I  =  Q 
120  COHTIHUE 


C  FIHD  THE  CURREHT  MISSILE  TYPE  IH  THE  LIST  OF  MISSILE  TYPES: 
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DO  130  Q  =  1 ,  MT 

IF  (MISS  .EQ.  MTyPE(Q))  K  =  Q 
130  COITIIUE 

C  ASSIGI  THE  FLIGHT  TIMES  FOR  TYPE  K  MISSILES  FROM  LAUJCH  FIELD  I  TO  EACH 
C  TARGET  COMPLEX  J: 

DO  140  J  =  1.  TC 

Td.J.K)  =  TEMP(J) 

140  COITIHUE 

GO  TO  110 

C  IF  THE  FIRST  HIE  CHARACTERS  OF  THE.  LIKE  ARE  'HEIGHTIHG',  THEN  THE  IFXT 
C  SECTION  CONTAINS  THE  WEIGHTING  FACTORS; 

ELSE  IF  (LINE(1:9)  .EQ.  'WEIGHTING')  THEN 
READdO,*) 

READ(10,1030)  W1 
READ(10,1030)  W2 

C  IF  THE  FIRST  10  CHARACTERS  OF  THE  LINE  ARE  'PERCENTAGE',  THE  NEXT  SECTION 
C  SPECIFIES  PCT: 

ELSE  IF  (LINE(1:10)  .EQ.  'PERCENTAGE')  THEN 
READClO,*) 

READ(10,1040)  PCT 

C  IF  THE  FIRST  14  CHARACTERS  OF  THE  LINE  ARE  'PERCENT  BACKUP',  THE  NEXT 
C  SECTION  SPECIFIES  PCTBO: 

ELSE  IF  (LIHE(l:14)  .EQ.  'PERCENT  BACKUP')  THEN 
READdO,*) 

READd0,1060)  PCTBU 
GO  TO  600 

C  OTHERWISE,  THE  FILE  'SILOATT.DAT'  HAS  AN  ERROR  OR  IS  THE  WRONG  FILE; 

ELSE 

0PEN(UNIT=16,  FILE='SILOATT.Er.R',  STATUS= ' UNKNOWN ' , 

A  ACCESS= 'SEQUENTIAL') 

WRITE(16,*)  'THIS  IS  THE  WRONG  INPUT  FILE' 

CLOSEdS) 

STOP 
END  IF 

GO  TO  10 

600  CLOSEClO) 

RETURN 

9991  0PEN(UNIT=16,  FILE= 'SILO ATT. ERR ' ,  STATUS= ' UNKNOWN ' , 

*  ACCESS=' SEQUENTIAL') 

WRITEdS,*)  'ERROR  IN  OPENING  FILE’ 

WRITEClS,*)  'ERROR  CODE  =  lERROR 
CLOSEdS) 
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STOP 


1000  F0RMAT(A80) 

1010  F0RMAT(A8.  IX,  A4.  4X,  13) 

1020  F0mT(A8.  IX,  A4,  4X,  10(F6.2,  3X)) 

1030  FORMATCSX.  F6.2) 

1040  F0RMAT(6X,  F4.2) 

1060  F0RMAT(8X,  F4.2) 

END 

C***********t***********************t********^***********i¥*****^*******i^****^^^ 
C******************************************************t*********************** 
SUBROUTINE  PaOCESS(BU,  CPLX,  LF.  LFMT,  K,  MADJ,  MAXLF,  MAXMT, 
k  MAXTC,  MAXTT,  KT.  N,  NADJ,  NUMBU,  NUMFAIR, 

k  NUMGOOD,  NUNPOOR,  NUKTAR,  PCTBU,  T,  TC, 

k  TCTT,  TMAX,  TMIN,  TT) 

c 

C  NAME;  PROCESS 
C 

C  REV  DATE  OF 

C  NO.  CHANGE  DESIGNER  CCR  DESCRIPTION 

C  -  -  - - -  - 

C  000  9FEB93  R.  PACE  ORIGINAL  RELEASE 

C 

C  DESCRIPTION:  PROCESSES  THE  DATA  FOR  THE  SILO  ATTACK  PLAN 
C 

C  CALLED  BY:  SILOATT 
C 

C  CALLS:  NONE 
C 

C  CALLING  SEQUENCE: 

C 

CIN 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 

CI/0 
C 
C 
C 

COUT 
C 


LF  -  NUMBER  OF  LAUNCH  FIELDS 

M  -  ARRAY  STORING  NUMBER  OF  MISSILES  BY  TYPE  AT  EACH  LAUNCH  FIELD 

MAXLF  -  MAX  NUMBER  OF  LAUNCH  FIELDS 

MAXMT  -  MAX  NUMBER  OF  MISSILE  TYPES 

MAXTC  -  MAX  NUMBER  OF  TARGET  COMPLEXES 

MAXTT  -  MAX  NUMBER  OF  TARGET  TYPES 

MT  -  THE  NUMBER  OF  MISSILE  TYPES 

N  -  ARRAY  STORING  NUMBER  OF  TARGETS  BY  TYPE  AT  EACH  TARGET  COMPLEX 

PCTBU  -  THE  PERCENT  OF  THE  AVAILABLE  MISSILES  AT  EACH  LAUNCH  FIELD 
THAT  MUST  BE  RESERVED  FOR  BACKUP 
TT  -  HUMBER  OF  TARGET  TYPES 

CPLX  -  ARRAY  STORING  THE  NAMES  OF  THE  TARGET  COMPLEXES 
TC  -  NUMBER  OF  TARGET  COMPLEXES 

T  -  ARRAY  STORING  FLIGHT  TIMES  OF  MISSILES  (BY  r/PE)  FROM  EACH  LAUNCH 

FIELD  TO  EACH  TARGET  COMPLEX 

BU  -  ARRAY  STORING  THE  INDICES  OF  THE  BACKUP  LAUNCH  FIELDS 

LFMT  -  THE  NUMBER  OF  LAUNCH  FIELD/MISSILE  TYPE  COMBINATIONS  THAT  ARE 
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C  POSSIBLE  (I.E.,  THE  HUMBER  OF  MISSILES  OF  THAT  TYPE  IS  HOT  ZERO) 

C  MADJ  -  ARRAY  STORIHG  THE  ADJUSTED  HUMBER  OF  MISSILES  OF  EACH  TYPE  AT 
C  EACH  LAUHCH  FIELD 

C  NADJ  -  ARRAY  STORING  THE  ADJUSTED  HUMBER  OF  TARGETS  OF  EACH  TYPE  AT 
C  EACH  TARGET  COMPLEX 

C  NUMBU  -  THE  HUMBER  OF  BACKUP  LAUNCH  FIELDS 

C  NUMFAIR  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  DO  HOT  CONTAIN  ANY  GOOD 

C  TARGETS  BUT  DO  COHTAIH  AT  LEAST  OHE  FAIR  TARGET 

C  MUMGOOD  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  COHTAIH  AT  LEAST  OHE  GOOD 
C  TARGET 

C  HUMPOOR  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  DO  HOT  COHTAIH  ANY  GOOD  OR 
C  FAIR  TARGETS 

C  NUMTAR  -  ARRAY  STORIHG  THE  TOTAL  HUMBER  OF  TARGETS  AT  EACH  TARGET  COMPLEX 

C  TCTT  -  THE  HUMBER  OF  TARGET  COHPLEX/TARGET  TYPE  COMBIHATIOHS  THAT  ARE 

C  POSSIBLE 

C  TMAX  -  THE  LARGEST  POSSIBLE  FLIGHT  TIME 

C  TMIH  -  THE  SMALLEST  POSSIBLE  FLIGHT  TIME 
C 

C  FILES: 

C 

C  SILOATT.DAT 
C  SILO ATT. ERR 
C 

C  LOCAL  VARIABLES: 

C 

C  AISl  -  STORES  RESPONSE  TO  WHETHER  OR  HOT  ANY  LAUHCH  FIELDS  ARE  FOR 
C  BACKUP  ONLY 

C  AIS2  -  STORES  RESPONSE  TO  WHETHER  OR  HOT  ANY  TARGET  COMPLEXES  SHOULD  BE 
C  COMBINED 

C  COMB  -  ARRAY  STORING  ID  OF  TARGET  COMPLEXES  TO  BE  COHBIHED 

C  LINE  -  STORES  THE  CONTENTS  OF  A  LINE  FROM  'SILOATT.DAT' 

C  NUHTGTl  -  THE  TOTAL  HUMBER  OF  TARGETS  AT  THE  FIRST  COMPLEX  TO  BE  COMBINED 

C  NUMTCT2  -  THE  TOTAL  NUMBER  OF  TARGETS  AT  THE  SECOND  COMPLEX  TO  BE  COMBINED 

C  TEMP  -  ARRAY  USED  TO  TEMPORARILY  STORE  THE  RESULT  OF  MULTIPLYING  EACH 
C  N(J,L)  BY  0.1 

C 

C  ERROR  MESSAGES:  NONE 
C 

C****'*'**'**'^******'*'************’''****'*'*********************' ***<•>******************* 

c 

C  LOCAL  VARIABLE  DECLARATIONS: 

INTEGER  BU(MAXLF),  C0MB(2),  LF,  LFMT,  M(MAXLF,MAXMT) , 
ft  MADJ(MAXLF,MAXMT),  MT,  M(MAXTC.MAXTT) ,  HADJ(HAXTC .MAXTT) , 

ft  HUMBU,  NUMFAIR.  NUMGOOD,  HUMPOOR,  HUMTAR(MAXTC) .  MUMTGTl, 

ft  NUMTGT2,  TC,  TCTT,  TT 

REAL  PCTBU,  T(MAXLF,MAXTC,MAXMT) .  TEMP (10, 5).  TMAX,  TMIN 

CHARACTER  ANSI ,  AHS2 
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CHARACTER*8  CPLX(MAXTC) 

,  CHARACTER*80  LIHE 

0PEI(UIIT=10,  FILE=* SILOATT.DAT*,  STATtIS=»OLD* . 
k  ACCESS- 'SEQUEITIAL'.  I0STAT=IER!10R.  ERR=991) 

10  READ (10. 1020}  LIRE 

C  IF  THE  FIRST  THREE  CHARACTERS  OF  A  LIHE  ARE  'ARE'.  THE  lEXT  SECTIOI  TELLS 
C  HBETHER  OR  HOT  AHY  LAUICH  FIELDS  ARE  FOR  BACKUP  OHLY. 

IF  (LIIE(1:3)  .EQ.  'ARE')  THEH .. 

REAO(IO.IOOO)  AISl 

C  IF  AHY  LAUHCH  FIELDS  ARE  FOR  BACKUP  OHLY.  READ  IN  HOV  NAHY  AHD  HHICH  OHES: 
IF  ((AHSl  .EQ.  'Y')  .OR.  (AHSl  .EQ.  'y'))  THEH 
READdO.*) 

READClO.*) 

READClO.*)  HUMBU 
READ(10.«) 

READClO.*) 

READdO.*)  (BU(I).  1  =  1.  HUMBU) 

C  OTHERWISE.  SKIP  THE  HEXT  SIX  LIHES: 

ELSE 

DO  20  I  =  1.  6 
READdO.*) 

20  COHTIHUE 

EHD  IF 

READdO,*) 

READClO,*) 

READdO,  1000)  AHS2 

C  IF  AHY  OF  THE  TARGET  COMPLEXES  SHOULD  BE  COMBINED,  READ  IH  WHICH  OHES  AHD 
C  WHAT  THE  NAME  OF  THE  COMBINED  COMPLEX  SHOULD  BE: 

IF  ((AHS2  .EQ.  'Y')  .OR.  (AHS2  .EQ.  'y'))  THEM 
READdO,  ♦) 

READClO,*) 

READClO,*)  CCOMB(I),  1=1,  2) 

READClO,*) 

READClO,*) 

READClO,  1010)  CPLXCCOMBCD) 

HUMTGTl  =  0 
HUMTGT2  =  0 

C  THIS  DO  LOOP  CALCULATES  THE  TOTAL  HUMBER  OF  TARGETS  AT  THE  TWO  TARGET 
C  COMPLEXES  THAT  WILL  BE  COMBINED: 

DO  30  L  =  1,  TT 

HUMTGTl  =  HUMTGTl  +  H(COMBd),L) 
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KUMTGT2  =  1IUMTGT2  +  N(C0MB(2),L) 

30  COHTIIUE 

C  THIS  IMBEDDED  DO  LOOP  CALCULATES  THE  WEIGHTED  AVERAGE  OF  THE  FLIGHT  TIMES 
C  FOR  EACH  LAUNCH  FIELD/MISSILE  TYPE  COMBIHATION  TO  THE  TWO  TARGET  COMPLEXES 
C  THAT  ARE  BEING  COMBINED: 

D060I=1,LF 

DO  40  K  =  1,  NT 

T(I,C0MB(1),X)  =  (HUMTGT1*T(I,C0MB(1),K) 
k  +  ■UMTGT2*T(I,C0MB(2),K))/(HUMTGT1  +  NDMTGT2) 

40  CONTINUE 

60  CONTINUE 

C  HERE  THE  PROGRAM  WILL  TREAT  THE  COMBINED  TARGET  COMPLEX  AS  HAVING  THE  LOWER 
C  INDEX  OF  THE  TWO  TARGET  COMPLEXES  IT  REPRESENTS.  AND  "SLIDES"  THE  OTHERS 
C  FORWARD  AS  APPROPRIATE: 

C  IF  THE  LARGER  OF  THE  TWO  INDICES  IS  NOT  THE  LAST  TARGET  COMPLEX,  THEN 
C  SLIDE  THE  FLIGHT  TIME  DATA  FORWARD  FOR  EACH  TARGET  COMPLEX  WITH  A  BIGGER 
C  INDEX: 

IF  (COMB (2)  .NE.  TC)  THEM 
DO  80  I  s  1,  LF 

DO  70  J  =  C0MB(2),  TC-1 
DO  60  K  =  1,  MT 

T(I.J.K)  =  T(I.J+1,K) 

60  CONTINUE 

CPLX(J)  =  CPLX(J+1) 

70  CONTINUE 

80  CONTINUE 

END  IF 

C  THIS  DO  LOOP  STORES  THE  NUMBER  OF  TARGETS  OF  EACH  TYPE  OF  THE  COMBINED 
C  COMPLEX  INTO  THE  ARRAY  PREVIOUSLY  USED  FOR  THE  COMPLEX  WITH  THE  LOWER  INDEX: 
DO  90  L  =  1 ,  TT 

N(C0HB(1).L)  =  H(C0MB(1),L)  +  N(C0MB(2),L) 

90  CONTINUE 

C  IF  THE  COMPLEX  WITH  THE  LARGER  INDEX  IS  NOT  THE  LAST  TARGET  COMPLEX,  THEN 
C  THE  DO  LOOP  SLIDES  THE  TARGET  DATA  FORWARD  FOR  ALL  COMPLEXES  WITH  A  URGER 
C  INDEX: 

IF  (C0MB(2)  .NE.  TC)  THEN 
DO  110  J  =  C0MB(2),  TC-1 
DO  100  L  =  1,  TT 

H(J,L)  =  N(J+1,L) 

100  CONTINUE 

110  CONTINUE 

END  IF 
TC  =  TC  -  1 
END  IF 
GO  TO  120 
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/ 

/ 


ELSE 

GO  TO  10 
EMD  IF 

C  HERE  THE  PROGRAM  CALCULATES  THE  HUMBER  OF  TARGET  FIELDS  THAT  HAVE 
C  GOOD.  FAIR,  AHD  POOR  TARGETS  RESPECTIVELY: 

120  HUMGOOD  =  0 
HUMFAIR  =0 
HUMPOOR  =  0 

C  THIS  DO  LOOP  ITERATES  THROUGH  EACH  TARGET  COMPLEX: 

DO  130  J  =  1,  TC 

C  IF  THERE  ARE  GOOD  TARGETS  AT  COMPLEX  J,  THEH  IICREASE  HUMGOOD  BY  1: 

IF  (I(J,1)  .GT.  0)  THEH 
HUMGOOD  =  HUMGOOD  +  1 

C  ELSE,  IF  THERE  ARE  FAIR  TARGETS  AT  COMPLEX  J,  THEH  HCr.EASE  HUMFAIR  BY  1: 
ELSEIF  (H(J,2)  .GT.  0)  THEH 
HUMFAIR  =  HUMFAIR  +  1 

C  ELSE,  THERE  ARE  OHLY  POOR  TARGETS  AT  COMPLEX  J.  SO  IHCREASE  HUMPOOR  BY  1: 
ELSE 

HUMPOOR  =  HUMPOLrt  +  1 
EHD  IF 
130  COKTIHUE 

C  HERE  THE  PROGRAM  FIHDS  TMIH  (THE  SMALLEST  FLIGHT  TIME)  AHD  TMAX 
C  (THE  LARGEST  FLIGHT  TIME) : 

TMIH  =  1000.0 
TMAX  =  0.0 

C  THIS  IMBEDDED  DO  LOOP  STEPS  THROUGH  ALL  COHBIHATIOHS  OF  LAUHCH  FIELD/ 

C  TARGET  COMPLEX/MISSILE  TYPE  LOOKTHG  FOR  THE  LARGEST  AHD  SMALLEST  T(I,J,K): 
DO  160  1  =  1,  LF 

DO  150  J  =  1,  TC 
DO  140  K  =  1,  MT 

C  IF  THERE  ARE  HO  MISSILES  OF  TYPE  K  AT  LAUHCH  FIELD  I,  T(I.J,K)  WILL  BE  0, 

C  TRY  THE  HEXT  MISSILE  TYPE: 

IF  (T(I,J,K)  .EQ.  0.0)  CO  TO  140 

C  IF  THE  CURREHT  T(I,J,K)  IS  GREATER  THAH  THE  CURREHT  TMAX,  SET  TMAX  EQUAL  TO 
C  THE  CURREHT  T(I,J,K): 

IF  (T(I,J,K)  .GT.  TMAX)  TMAX  =  T(I,J,K) 

C  IF  THE  CURREHT  T(I,J,K)  IS  LESS  THAH  THE  CURREHT  TMIH.  SET  TMIH  EQUAL  TO 
C  THE  CURREHT  T(I,J,K): 

IF  (T(I.J.K)  .LT.  TMIH)  TMIH  =  T(I,J,K) 

140  COHTIHUE 

ISO  COHTIHUE 
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160  COHTIHUE 


/ 


C  HERE  THE  PROGRAM  CALCULATES  THE  NUMBER  OF  LAUNCH  FIELD/MISSILE 
C  TYPE  COMBINATIONS  EXIST: 

LFMT  =  0 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  190  I  =  1,  LF 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  170  P  =  1.  NUMBU 

C  IF  THE  CURRENT  LAUNCH  FIELD  (I)  IS  ONE  OF  THE  BACKUP  ONES,  GO  TO  THE 
C  NEXT  LAUNCH  FIELD: 

IF  (I  .EQ.  BU(P))  GO  TO  190 
170  CONTINUE 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  180  K  =  1,  MT 

C  IF  THERE  ARE  NO  MISSILES  OF  TYPE  K  AT  LAUNCH  FIELD  I.  GO  TO  THE  NEXT 
C  MISSILE  TYPE.  OTHERWISE,  INCREASE  LFMT  BY  ONE: 

IF  (H(I,K)  .EQ.  0)  THEN 
GO  TO  180 

'  ELSE 

LFMT  =  LFMT  +  1 
END  IF 

180  CONTINUE 

190  CONTINUE 
TCTT  =  0 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  210  J  =  1,  TC 

C  STEP  THROUGH  EACH  TARGET  TYPE: 

DO  200  L  =  1,  TT 

C  IF  THERE  ARE  NO  TARGETS  OF  TYPE  L  AT  TARGET  COMPLEX  J,  GO  TO  THE  NEXT 
C  TARGET  TYPE.  OTHERWISE,  INCREASE  TCTT  BY  1: 

IF  (N(J,L)  .EQ.  0)  THEN 
GO  TO  200 

ELSE 

TCTT  =  TCTT  +  1 
END  IF 

200  CONTINUE 
210  CONTINUE 

C  HERE  THE  PROGRAM  ADJUSTS  THE  NUMBER  OF  MISSILES  AT  EACH  LAUNCH  FIELD. 

C  FOR  FAIR  MISSILES,  THE  PROGRAM  WILL  ONLY  CONSIDER  507,  (ROUNDED  DOWN) 

C  OF  THE  MISSILES  AVAILABLE  FOR  THE  ALLOCATION: 
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C  STEP  THROUGH  EACH  LAUHCH  FIELD: 
DO  230  I  =  1,  LF 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 
DO  220  K  =  1,  MT 

C  SET  THE  MADJd.K)  =  H(I,K): 

MADJd.K)  =  Md.K) 
220  COHTIHUE 
230  COHTIHUE 


C  IF  YOU  DO  HOT  WISH  TO  EHFORCE  THAT  EACH  TARGET  ATTACKED  BY  A  FAIR  TARGET  IH 
C  THE  IHITIAL  WAVE  BE  ATTACKED  BY  A  FAIR  MISSILE  FROM  THE  SAME  LAUHCH  FIELD  IH 
C  THE  SECOHD  WAVE,  THEH  COHMEHT  OUT  THIS  SECTIOH: 

C  STEP  THROUGH  EACH  LAUHCH  FIELD: 

DO  240  I  =  1.  LF 

C  FOR  FAIR  MISSILES  (K  =  2),  SET  MADJd,2)  EQUAL  TO  OHE-HALF  OF  Md,2), 

C  ROUHDIHG  DOWH  FOR  FRACTIOHS: 

MADJd,2)  =  IHT(0.S*MADJd,2)) 

240  COHTIHUE 

C**’***’**»*********************************************************************» 


C************************************t***************************************** 

C  IF  YOU  DO  HOT  WISH  TO  HAVE  A  PCTBU  BACKUP  OF  EACH  MISSILE  TYPE  AT  EACH 
C  LAUHCH  FIELD,  THEN  COMMENT  OUT  THIS  SECTIOH: 

C  SINCE  THERE  SHOULD  BE  A  PCTBU  BACKUP  OF  EACH  MISSILE  TYPE  AT  EACH  LAUNCH 
C  FIELD,  MULTIPLY  BY  (1  -  PCTBU)  AHD  ROUND  DOWN: 

C  STEP  THROUGH  EACH  LAUHCH  FIELD: 

DO  260  I  =  1,  LF 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  250  K  =  1,  MT 

C  SET  MADJd,K)  EQUAL  TO  0.9  TIMES  MADJd,K),  ROUHDIHG  DOWH  FOR  FRACTIOHS: 
MADJd,K)  =  IHTCd.O  -  PCTBU)*HADJd,K)) 

250  CONTINUE 
260  COHTIHUE 


C  HERE  THE  PROGRAM  DIVIDES  THE  HUMBER  OF  TARGETS  AT  EACH  COMPLEX  BY  10, 
C  ROUHDIHG  UP  FOR  ANY  FRACTION: 
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C  IF  EACH  TYPE  OF  MISSILE  CAH  HAVE  A  DIFFEREHT  HUMBER  OF  WARHEADS,  THEN  COMMENT 
C  OUT  THIS  SECTION  AND  INCLUDE  THE  NEXT  ONE: 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  280  J  =  1,  TC 

C  STEP  THROUGH  EACH  TARGET  TYPE: 

DO  270  L  =  1.  TT 

TEMP(J,L)  =  0.1*N(J.L) 

C  IF  0.1*N(J,L)  IS  NOT  AN  INTEGER,  ROUND  UP  FOR  THE  FRACTION: 

IF  (AM0D(TEMP(J,L),1.0)  .GT.  0)  THEN 
NADJ(J,L)  =  INT(TEMP(J,L))  +  1 

C  OTHERWISE,  SET  NADJ(J,L)  =  0.1*N(J,L): 

ELSE 

NADJ(J,L)  =  INT(TEMP(J,L)) 

END  IF 

270  CONTINUE 
280  CONTINUE 

C^^ttft***************************************************^'********************* 

C  DO  280  J  =  1,  TC 

C  DO  270  L  =  1,  TT 

C  NADJ(J,L)  =  N(J,L) 

C  270  CONTINUE 
C  280  CONTINUE 

C****************)t****r¥**************«******-tt-t*******0******<l‘*********^im****** 


C  CALCULATE  THE  NUMBER  OF  TARGETS  AT  EACH  COMPLEX: 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  300  J  =  1,  TC 
NUMTAR(J)  =  0 

C  STEP  THROUGH  EACH  TARGET  TYPE,  SETTING  THE  NUMBER  OF  TARGETS  AT  COMPLEX  J 
C  EQUAL  TO  THE  CURRENT  NUMBER  PLUS  THE  ADJUSTED  NUMBER  OF  TARGETS  OF  TYPE  L 
C  AT  COMPLEX  J: 

DO  290  L  =  1,  TT 

NUMTAR(J)  =  NUMTAR(J)  +  NADJ(J,L) 

290  CONTINUE 
300  CONTINUE 

CLOSE (10) 

RETURN 

991  OPENCUNIT: 15,  FILE= 'SILO ATT. ERR' ,  ACCESS= ' SEQUENTIAL ' , 

A  STATUS=' UNKNOWN') 

WRITEClS,*)  'ERROR  IN  OPENING  FILE  SILOATT.DAT' 

WRITEdS,*)  'ERROR  CODE  =  ',  lERROR 
CLOSEClS) 
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STOP 


1000  FORMAT(Al) 

1010  F0RMAT(A8) 

1020  F0RMAT(A80) 

LSD 

C*'*****************************4******t*:***it**************************^********* 

SUBROUTIIE  FORMULCBU,  LF,  LFMT.  MADJ.  KAXLF,  MAXMT,  SAXTC,  SAXTT, 
ft  MT,  lADJ,  lUNBU,  BUKFAIR.  BUNGOOD,  BUHPOOR, 

ft  lUMTAR,  PCT,  T.TC,  TCTT,  TMAX,  THIS,  TT. 

ft  Wl.  V2) 

Qt***************************************^************************************* 

c 

C  SAME:  FORMUL 

C 

C  REV  DATE  OF  I 

C  SO.  CHARGE  DESIGNER  CCR  DESCRIPTIOR 

C  -  J - 

C  000  9FEB93  R.  PACE  ORIGINAL  RELEASE 

C  I 

C  DESCRIPTION:  WRITES  THE  DATA  INTO  AN  MPS-IFORHAT  FILE 
C  i 

C  CALLED  BY:  SILOATT  | 

C  i 

C  CALLS:  NONE  I 

C  i 

C  CALLING  SEQUENCE:  { 

C  I 

CIS  BU  -  ARRAY  STORING  THE  INDICES  JF  ijHE  BACKUP  LAUNCH  FIELDS 
C  LF  -  NUMBER  OF  UUNCH  FIELDS  I 

C  LFMT  -  THE  NUMBER  OF  LAUNCH  FIELD/MISSILE  TYPE  COMBINATIONS  THAT  ARE 
C  POSSIBLE  (I.E.,  THE  NUMBER  OF  MISSILES  OF  THAT  TYPE  IS  HOT  ZERO) 

C  MADJ  -  ARRAY  STORING  THE  ADJUSTED  NUMBER  OF  MISSILES  OF  EACH  TYPE  AT 
C  EACH  LAUNCH  FIELD 

C  MAXLF  -  MAX  NUMBER  Or  LAUNCH  FIELDS 

C  MAXMT  -  MAX  NUMBER  OF  MISSILE  TYPES 

C  MAXTC  -  MAX  NUMBER  OF  TARGET  COMPLEXES 

C  MAXTT  -  MAX  HUMBER  OF  TARGET  TYPES 

C  NT  -  THE  NUMBER  OF  MISSILE  TYPES 

C  HADJ  -  ARRAY  STORING  THE  ADJUSTED  HUMBER  OF  TARGETS  OF  EACH  TYPE  AT 
C  EACH  TARGET  COMPLEX 

C  BUMBU  -  THE  NUMBER  OF  BACKUP  LAUNCH  FIELDS 

C  NUNFAIR  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  DO  NOT  CONTAIN  ANY  GOOD 

C  TARGETS  BUT  DO  CONTAIN  AT  LEAST  ONE  FAIR  TARGET 

C  NUMGOOD  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  CONTAIN  AT  LEAST  ONE  GOOD 
C  TARGET 

C  NUMPOOR  -  THE  HUMBER  OF  TARGET  COMPLEXES  THAT  DO  HOT  CONTAIN  ANY  GOOD  OR 
C  FAIR  TARGETS 
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C  HUMTAR  -  ARRAY  STORING  THE  TOTAL  HUMBER  OF  TARGETS  AT  EACH  TARGET  COMPLEX 

C  PCT  -  PERCENT  OF  THE  MINIMUM  FLIGHT  TIME  THAT  THE  DURATION  IS  ALLOWED 

C  TO  CE 

C  T  -  ARRAY  STORING  FLIGHT  TIMES  OF  MISSILES  (BY  TYPE)  FROM  EACH  LAUNCH 

C  FIELD  TO  EACH  TARGET  COMPLEX 

C  TC  -  NUMBER  OF  TARGET  COMPLEXES 

C  TCTT  -  THE  NUMBER  OF  TARGET  COHPLEX/TARGET  TYPE  COMBINATIONS  THAT  ARE 

C  POSSIBLE 

C  TMAX  -  THE  LARGEST  POSSIBLE  FLIGHT  TIME 

C  TMIH  -  THE  SMALLEST  POSSIBLE  FLIGHT  TIME 

C  TT  -  NUMBER  OF  TARGET  TYPES 

C  W1  -  THE  RELATIVE  WEIGHT  OF  MINIMIZING  THE  EARLIEST  FLIGHT  TIME 

C  W2  -  THE  RELATIVE  WEIGHT  OF  MINIMIZING  THE  DURATION  OF  THE  ATTACK 

C 

C  FILES:  HONE 
C 

C  LOCAL  VARIABLES: 

C 

C  COUNT  -  COUNTER 

C  COUNT 1  -  COUNTER 

C  I,  J,  K,  L  -  COUNTERS 

C  TERROR  -  STORES  THE  VALUE  OF  lOSTAT 

C  ITERLIM  -  UPPER  LIMIT  ON  THE  NUMBER  OF  LP  ITERATIONS  TO  BE  PERFORMED 

C  BY  ZOOM 

C  MAXNODES  -  UPPER  LIMIT  ON  THE  NUMBER  OF  NODES  ALLOWED  IN  THE  BRANCH  AND 
C  BOUND  TREE  IN  ZOOM 

C  MAXSAVE  -  UPPER  LIMIT  ON  THE  HUMBER  OF  WARM  BASES  SAVED  DURING  THE 
C  BRANCH  BOUND  SEARCH  IN  ZOOM 

C  NUMIHT  -  HUMBER  OF  0/1  VARIABLES  IN  THE  MIP 

C  NUMROWS  -  HUMBER  OF  ROWS  IN  THE  HIP 

C  HUMVARS  -  NUMBER  OF  VARIABLES  IN  THE  MIP  (EXCLUDING  SLACKS) 

C  P,  Q  -  COUNTERS 

C  R  -  R7  -  COUNTERS 

C 

C  ERROR  MESSAGES r  HONE 
C 

c 

C  LOCAL  VARIABLE  DECLARATIONS; 

INTEGER  BU(MAXLF),  COUNT,  COUNTl,  I,  lERROR,  ITERLIM,  J,  K,  L,  LF, 
ft  LFMT,  HADJ(MAXLF.HAXNT).  MAXNODES,  MAXSAVE,  MT, 

ft  NADJ(HAXTC,MAXTT),NUMBU,  HUMFAIR,  NUMGOOD,  NUMINT, 

ft  HUMPOOR,  NUMROWS,  NUMTAR(MAXTC) ,  HUMVARS,  P,  Q,  R,  R1 , 

ft  R2,  R3,  R4,  RS,  R6,  R7,  TC,  TCTT,  TT 

REAL  PCT,  T(MAXLF,MAXTC,MAXMT),  TMAX,  TMIN,  Wl,  W2 


0PEN(UHIT=40,  FILE=’SILOATT.MPS’ ,  ACCESS= 'SEQUENTIAL ’ , 
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*  STATUS='UHKHOWH'.  rOSTAT=IERRQR,  ERR=991) 

C  THE  FIRST  LINE  CONTAINS  THE  NAME  OF  THE  MODEL.  THIS  IS 
C  FOLLOWED  BY  THE  'ROWS*  SECTION: 

10  WRITE(40.1000),  'NAME',  'SILO  ATTACK  PLAN(  MIN)',  'ROWS', 

A  'N',  'OBJ' 

q  =  0 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  30  J  =  1,  TC 

C  STEP  THROUGH  EACH  TARGET  TYPE: 

DO  20  L  =  1,  TT 

C  IF  THERE  ARE  TARGETS  OF  TYPE  L  AT  TARGET  COMPLEX  J,  THEM  WRITE  A  "ROWS"  LINE 
C  FOR  AN  EQUALITY  CONSTRAINT  FOR  THAT  TARGET  COMPLEX/TARGET  TYPE  COMBINATION : 
IF  (HADJ(J,L)  .GT.  0)  THEN 

q  =  Q  +  1 

WRITE(40,1010),  'E',  'C',  Q 
END  IF 

20  CONTINUE 
30  CONTINUE 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  60  I  =  1,  LF 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  40  P  =  1,  NUMBU 

C  IF  THE  CURRENT  LAUNCH  FIELD  IS  A  BACKUP  ONE,  GO  TO  THE  NEXT  FIELD: 

IF  (I  .EQ.  BU(P))  GO  TO  60 
40  CONTINUE 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  SO  K  =  1,  MT 

C  IF  THERE  ARE  MISSILES  OF  TYPE  K  AT  LAUNCH  FIELD  I,  THEN  WRITE  A  "ROWS"  LINE 
C  FOR  A  LESS  THAN  OR  EQUAL  TO  CONSTRAINT  FOR  THAT  LAUNCH  FIELD/MISSILE  TYPE 
C  COMBINATION: 

IF  (MADJ(I,K)  .GT.  0)  THEN 

Q  =  Q  +  1 

WRITE(40,1010),  'L',  'C',  Q 
END  IF 

SO  CONTINUE 

60  CONTINUE 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  100  I  =  1,  LF 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  70  P  =  1,  NUMBU 
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C  IF  THE  CURRENT  LAUNCH  FIELD  IS  A  BACKUP  LAUNCH  FIELD,  GO  THE  NEXT  ONE: 

IF  (I  .EQ.  BU(P))  GO  TO  100 
70  CONTINUE 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  90  J  =  1,  TC 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  80  K  =  1.  MT 

C  IF  THERE  ARE  MISSILES  OF  TYPE  K  AT  UUNCH  FIELD  I.  THEM  WRITE  A  "ROWS"  LINE 
C  FOR  A  LESS  THAN  OR  EQUAL  TO  CONSTRAINT  FOR  THAT  LAUNCH  FIELD/TARGET  COMPLEX/ 
C  MISSILE  TYPE  COMBINATION: 

IF  (HADJd.K)  .GT.  0)  THEN 
Q  =  g  +  1 

HRITE(40,1010),  *L',  »C',  Q 
END  IF 

80  CONTINUE 

90  CONTINUE 

100  CONTINUE 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  140  I  s  1,  LF 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  110  P  =  1,  MUMBU 

C  IF  THE  CURRENT  LAUNCH  FIELD  IS  A  BACKUP  LAUNCH  FIELD,  GO  TO  THE  NEXT  ONE: 

IF  (I  .EQ.  BU(P))  GO  TO  140 
110  CONTINUE 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  130  J  =  1 ,  TC 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  120  K  =  1,  MT 

C  IF  THERE  ARE  MISSILES  OF  TYPE  K  AT  LAUNCH  FIELD  I,  THEN  WRITE  A  "ROWS"  LINE 
C  FOR  A  GREATER  THAN  OR  EQUAL  TO  CONSTRAINT  FOR  THAT  LAUNCH  FIELD/TARGET 
C  COMPLEX/MISSILE  TYPE  COMBINATION: 

IF  (MADJ(I,K)  .GT.  0)  THEN 
Q  =  Q  +  1 

WRITE(40,1010) ,  'G',  'C',  Q 
END  IF 

120  CONTINUE 

130  CONTINUE 
140  CONTINUE 
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C  IF  THE  USER  DOES  HOT  WANT  TO  USE  S(I,J)  VARIABLES,  COHMEHT  OUT  THIS 
C  SECTIOH  AND  INCLUDE  THE  NEXT: 

C  STEP  THROUGH  EACH  LAUNCH  FIELD/TARGET  COMPLEX  COMBINATION  (EXCLUDING 
C  BACKUPS) ,  AND  WRITE  A  "ROWS"  LINE  FOR  A  LESS  THAN  OR  EQUAL  TO  CONSTRAINT 
C  FOR  EACH  ONE: 

DO  160  P  =  1.  (LF  -  NUMBU)*TC 

q  =  q  +  1 

HRITE(40,1010),  'L‘,  *C' ,  Q 
160  CONTINUE 

C**********^^^*^i^**tf^tr^*^t^t*t^^**^4******************************t*************** 

C  DO  166  I  =  1,  LF 
C  DO  160  P  =  1,  NUMBU 

C  IF  (I  .Eq.  BU(P))  GO  TO  166 

C  160  CONTINUE 

C  DO  163  J  =  1,  TC 

C  DO  162  K  =  1.  NT 

C  IF  (HADJd.K)  .EQ.  0)  GO  TO  162 

c  q  =  q  +  1 

C  HRITE(40,1010),  *L*.  'C',  q 

C  162  CONTINUE 

C  163  CONTINUE 

C****************************************************************************** 

C  IF  THE  USER  DOES  NOT  WANT  TO  USE  S(I,J)  VARIABLES,  COMMENT  OUT  THIS 
C  SECTIOH: 

C  STEP  THROUGH  EACH  LAUNCH  FIELD/TARGET  COMPLEX  COMBINATION  (EXCLUDING 
C  BACKUPS) ,  AND  WRITE  A  "ROWS"  LINE  FOR  A  GREATER  THAN  OR  EQUAL  TO  CONSTRAINT 
C  FOR  EACH  ONE: 

DO  160  P  =  1,  (LF  -  NUMBU) ♦TC 

q  =  q  +  1 

WRITE(40,1010),  'G',  'C' ,  Q  i 

160  CONTINUE  1 

C*****************«i*******************************************f'***************m 

Qt****************************************************************************i 

C  IF  THE  USER  DOES  HOT  WANT  TO  USE  S(I,J)  VARIABLES,  COHMEHT  OUT  THIS  \ 

C  SECTION:  \ 

C  STEP  THROUGH  EACH  TARGET  COMPLEX,  AND  WRITE  A  “ROWS"  LINE  FOR  AH  EQUALITY 
C  CONSTRAINT  FOR  EACH  ONE: 

DO  170  J  =  1,  TC 

q  =  q  +  1 

C  IF  THE  USER  HANTS  TO  ALLOW  A  TARGET  COMPLEX  TO  BE  ATTACKED  BY  TWO  LAUNCH 
C  FIELDS,  COMMENT  OUT  THIS  STATEMENT  AND  INCLUDE  THE  NEXT  ONE: 

WRITE (40, 10 10),  'E',  'C' ,  Q 
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C  WRITE(40,1010)  'L*.  *C' .  Q 

l^^C^^^^^i^ll^t***********************^*’******************************************** 

17C  COKTIHUE 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  210  I  =  1,  LF 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  180  P  =  1,  lUMBU 

C  IF  THE  CURRENT  LAUNCH  FIELD  IS  A  BACKUP  ONE,  GO  TO  THE  NEXT  LAUNCH  FIELD: 

IF  (I  .EQ.  BU(P))  GO  TO  210 
180  CONTINUE 

C  STEF  THROUGH  E^CH  TARGET  COMPLEX: 

DO  200  J  =  1,  TC 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  190  K  =  1.  MT 

C  IF  THERE  ARE  MISSILES  OF  TYPE  K  AT  LAUNCH  FIELD  I,  THEN  WRITE  A  "ROHS"  LINE 
C  FOR  A  GREATER  THAN  OR  EQUAL  TO  CONSTRAINT  FOR  EACH  LAUNCH  FIELD/TARGET 
C  COMPLEX/MISSILE  TYPE  COMBINATION: 

IF  (MADJd.K)  .GT.  0)  THEN 

q  =  q  +  1 

WRITE(40,1010),  *G',  'C',  Q 
END  IF 

190  CONTINUE 

200  CONTINUE 
210  CONTINUE 

C  STEP  THROUGH  EACH  TARGET  COMPLEX  AND  WRITE  A  "ROWS"  LINT  FOR  A  GREATER 
C  THAN  OR  EQUAL  TO  CONSTRAINT  FOR  EACH  ONE: 

DO  220  J  =  1,  TC 
q  =  Q  +  1 

WRITE(40,1010),  'G',  ’C'.  Q 
220  CONTINUE 

Q  =  q  +  1 

WRITE(40,1010),  'E',  *C’,  Q 

C  STEP  THROUGH  EACH  TARGET  COMPLEX  AND  WRITE  A  "ROWS"  LINE  FOR  A  LESS  THAN 
C  OR  EQUAL  TO  CONSTRAINT  FOR  EACH  ONE: 

DO  22S  J  =  1,  TC 
Q  =  Q  +  1 

WRITE(40,1010) ,  'L',  'C’,  Q 
225  CONTINUE 

C  STEP  THROUGH  EACH  TARGET  COMPLEX  AMD  WRITE  A  "ROWS"  LINE  FOR  A  GREATER  THAN 
C  OR  EQUAL  TO  CONSTRAINT  FOR  EACH  ONE: 

DO  230  J  =  1,  TC 

q  =  q  +  1 
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WRITE(40.1010),  'G',  'C',  Q 
230  CONTINUE 
q  =  Q  +  1 

WRITE(40,1010).  'E',  'C,  Q 
q  =  p  +  1 

WRITE(40,1010),  'L*.  'C\  q 

C  STEP  THROUGH  EACH  COMBINATION  OF  TARGET  COMPLEXES  WITH  GOOD  TARGETS  AND 
C  TARGET  COMPLEXES  WITH  FAIR  TARGETS  AID  WRITE  A  "ROWS”  LINE  FOR  A  LESS  THAI  OR 
C  EqUAL  TO  CONSTRAINT  FOR  EACH  ONE: 

DO  240  P  =  1,  NUMGOOD^IUMFAIR 
q  =  q  +  1 

WRITE(40,1010).  'L«,  'C',  q 
240  CONTINUE 

C  STEP  THROUGH  EACH  COMBINATION  OF  TARGET  COMPLEXES  WITH  FAIR  TARGETS  AND 
C  TARGET  COMPLEXES  WITH  POOR  TARGETS  AID  WRITE  A  "ROWS"  LINE  FOR  A  LESS  THAN  OR 
C  EqUAL  TO  CONSTRAINT  FOR  EACH  ORE: 

DO  250  P  =  1.  HUMFAIR*NUMPOOR 

q  =  q  +  1 

WRITE(40,1010),  *L',  'C*.  q 
260  CONTINUE 

WBITE(40,1020),  'E*.  'C\  q  +  1,  'E*.  »C',  q  +  2 

C**************************^*»******t************«*******m************t*^*4i**** 

C  IF  THE  USER  WANTS  TO  ALLOW  A  TARGET  COMPLEX  TO  BE  ATTACKED  BY  TWO  LAUNCH 
C  FIELDS,  ADD  IK  THE  FOLLOWING  STATEMENT: 

C  WRITE(40,1010)  *L*,  'C',  q  +  3 

C**************************************************»******^***^****^^n^^i^^^^*^*tf 

C****************************************************************************** 

C  IF  THE  SOLUTION  IS  UNBALANCED,  THE  USER  NAT  INCLUDE  THE  NEXT  SECTION: 

c  q  =-q  +  1 
C  DO  256  I  =  1,  LF 
C  DO  262  P  =  1,  NUMBU 

C  IF  (I  .Eq.  BU(P))  GO  TO  266 

c  q  =  q  +  1 

C  ,  WRITE(40,1010)  'L',  'C',  q 

C  262  CONTINUE 
C  265  CONTINUE 


c***** ♦******♦♦*»♦♦*♦**♦*♦*♦*♦♦♦♦•♦♦♦♦♦♦♦*♦♦*♦♦♦*♦♦•♦♦♦***♦•♦♦*♦♦♦*♦*♦*♦♦♦♦♦*** 

C  IF  THE  USER  DOES  NOT  WANT  EACH  LAUNCH  FIELD  TO  HAVE  A  FCTBU  BACKUP,  THEN  ADD 
C  THE  NEXT  SECTION  (ALSO  DECLARE  TOTMIS)  : 

C  TOTMIS  =0 

C  DO  267  I  =  1,  LF 

C  DO  266  K  =  1,  MT 

C  TOTMIS  =  TOTMIS  +  NADJ(I,K) 

C  256  CONTINUE 
C  267  CONTINUE 
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C  Q  =  Q  +  1 

C  WRITE(40,1010)  'L'.  »C',  Q 


C  THE  lEXT  SECTIOH  IS  THE  'COLUMNS*  SECTION; 

HRITE(40.1110),  'COLUMNS’ 

C  VARIABLE  X(I,J.K.L): 

COUNT  =0 
R  =  TCTT 

R1  =  LFMT  +  TCTT 

C****************************************************************************** 
C  IF  THE  SOLUTION  IS  UNBALANCED,  ADD  THE  FOLLOVIHG  ASSIGNMENT  STATEMENT: 

C  j  R2  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  NUMBU)*TC  +  3*TC 
C  ^  *  +  NUMGOOD*HUHFAIR  +  NUMFAIR*NUMPOOR  +  5 

C**************t********-************************t****************************** 

Qtl0*****i‘****t*t**’Hf*’¥’H‘*********************<l‘*****************^*************** 

C  IF  THE  USER  DOES  HOT  WANT  EACH  LAUNCH  FIELD  TO  HAVE  A  PCTBU  BACKUP,  ADD  THE 
C  FOLLOWING  ASSIGNMENT  STATEMENT: 

C  :  R2  =  LFMT  ♦  TCTT  +  3*LFMT*TC  +  2*(LF  -  HUMBU)*TC  +  3*TC 

C  *  +  HUMGOOD*NUMFAIR  +  NUMFAIR*NUMPOOR  +  S 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

I  DO  300  I  =  1,  LF 
COUNTl  =  0 

I  q  =  o 

I 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  260  P  =  1,  NUMBU 

C  IF  THE  CURRENT  LAUNCH  FIELD  IS  A  BACKUP  ONE,  GO  TO  THE  NEXT  LAUNCH  FIELD: 

IF  (I  .EQ.  BU(P))  GO  TO  300 
260  CONTINUE 

R  =  R  +  1 


C****************************************************************************** 

C  IF  THE  SOLUTION  IS  UNBALANCED,  ADD  THE  FOLLOWING  STATEMENT: 

C  R2  =  R2  +  1 

C****************************************************************************** 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  290  J  =  1,  TC 
Q  =  q  +  1 
COUNT  =  0 

IF  (COUNTl  .GT.  1)  R  =  R  -  1 
COUNTl  =  0 
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C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  280  K  =  1,  MT 


C  IF  THERE  ARE  HO  MISSILES  OF  TYPE  K  AT  LAUNCH  FIELD  I,  GO  TO  THE  HEXT  MISSILE 
C  TYPE: 

IF  (MADJd.K)  .EQ.  0)  GO  TO  280 
IF  (COUNT  .GT.  1)  Q  =  Q  -  1 
COUNT  =0 

COUNT!  =  COUNT!  +  ! 

IF  (COUNT!  .GT.  !)  R  =  R  +  ! 

R!  =  R!  +  ! 


C  STEP  THROUGH  EACH  TARGET  TYPE: 

DO  270  L  =  1,  TT 


C  IF  THERE  ARE  NO  TARGETS  OF  TYPE  L  AT  COMPLEX  J,  GO  TO  THE  NEXT  TARGET  TYPE: 
IF  (NADJ(J,L)  .EQ.  0)  GO  TO  270 
COUNT  =  COUNT  +  ! 

IF  (COUNT  .GT.  !)  P  =  Q  ♦  1 


C  OTHERWISE,  X(I,J,K,L)  WILL  BE  IN  CONSTRAINT  Q  WITH  A  COEFFICIENT  OF  !.0,  AND 
C  HILL  BE  IN  CONSTRAINTS  R,  Rl,  AND  R1  +  LFMT*TC  WITH  COEFFICIENTS  OF  !.0. 

C  -!.0,  AND  -!.0  RESPECTIVELY: 

*X'. 

'X», 

'X’. 

’X» 


HRITE(40.1030). 

HRITE(40,!030), 

HRITE(40,!030), 

HRITE(40,1030), 


I.  J, 
I, 

I. 


K. 


L. 

J.  K,  L, 
J.  K,  L, 
K. 


I,  3, 

LFMT*TC.  -!.0 


L. 


'C»,  Q.  !.0 
»C»,  R,  !.0 
»C»,  R!,  -!.0 
*C'  R!  + 


C***********************************************»*»’*******tf***<ti************<^*** 

C  IF  TEE  SOLUTION  IS  UNBALANCED,  DECLARE  AND  ASSIGN  A  VALUE  TO  JSTAR  AND  ADD 
C  THE  FOLLOWING  STATEMENT: 

C  IF  (J  .EQ.  JSTAR)  HRITE(40, !030)  'X’,  i,  JSTAR, 

C  k  K,  L,  'C’,  R2,  !.0 

C********************************************************************rH„H,*****^ 

C  IF  THE  USER  DOES  NOT  WANT  TO  HAVE  A  PCTBU  AT  EACH  LAUNCH  FIELD,  ADD  THE 
C  FOLLOWING  STATEMENT: 

C  WRITE(40,!030)  'X»,  1,  J,  K,  L,  'C',  R2  +  !,  !.0 

C*******************************m*t************************tm***t**t*r^*^****t** 

270  CONTINUE 

280  CONTINUE 

290  CONTINUE 
300  CONTINUE 

HRITE(40,!040),  ’INT’ 


C4<**«**^****4i*******4i4i4i*i*4i*4i*^*4i4i«4i*4t***********i»4>*4i:*4>*<»**^*4>*4i***4t*4>«**i^**«*4i4t 

C  IF  THE  USER  DOES  NOT  WANT  TO  USE  S(I.j)  VARIABLES.  THEN  COMMENT  OUT  THIS 
C  SECTION  AND  INCLUDE  THE  NEXT.  YOU  MUST  ALSO  DECLARE  THE  VARIABLES  'M'  AND 
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C  ROW(I,J,K).  YOU  MUST  ALSO  CHAHGE  SOME  OTHER  SECTIONS  AS  ANNOTATED  BELOW: 

C  VARIABLE  Y(1,J,K): 

R  =  LFHT  +  TCTT 

R1  =  LFMT  +  TCTT  +  2*LFMT*TC 

R2  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  IUHBU)*TC  +  TC 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  340  I  =  1,  LF 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  310  P  =  1,  NUMBU 

C  IF  THE  CURRENT  LAUNCH  FIELD  IS  A  BACKUP  ONE,  Gl  TO  THE  NEXT  LAUNCH  FIELD: 

IF  (I  .EQ.  BU(P))  GO  TO  340 
310  CONTINUE 

R2  =  R2  -  LFMTfTC 
R3  =  0 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  330  J  =  1,  TC 
R1  =  R1  +  1 
B2  =  R2  -  R3 
R3  =  0 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  320  K  =  1,  MT 

C  IF  THERE  ARE  NO  MISSILES  OF  TYPE  K  AT  LAUNCH  FIELD  I,  GO  TO  THE  NEXT  MISSILE 
C  TYPE: 

IF  (MADJd.K)  .EQ.  0)  GO  TO  320 
R  =  R  +  1 
R2  =  R2  +  1 
R3  =  R3  +  1 

C  OTHERWISE.  Y(I,J,K)  WILL  BE  IN  CONSTRAINT  R  WITH  A  COEFFICIENT  OF  1.0, 

C  CONSTRAINT  R  +  LFMT*TC  WITH  A  COEFFICIENT  OF  THE  HUMBER  OF  TARGETS  AT  COMPLEX 
C  J,  CONSTRAINT  R1  WITH  A  COEFFICIENT  OF  -1.0,  AND  CONSTRAINT  Rl  + 

C  (LF  -  NUMBU) *TC  WITH  A  COEFFICIENT  OF  -1.0: 

WRITE (40, 1050),  'Y',  I,  J,  K,  'C',  R,  1.0 

WRITE(40,1060),  'Y',  I,  J,  K,  'C'.  R  +  LFMT*TC, 

k  REAL(HUMTAR(J)) 

WRITE(40,1050),  ’Y',  I,  J,  K,  'C'.  Rl,  -1.0 

WRITE(40.10E'>),  ’Y>,  I,  J,  K,  ’C*.  Rl  +  (LF  -  NUMBU) 

k  ♦TC,  -1.0 

WRITE(40.1050),  'Y',  I.  J,  K.  'C'.  R2.  -T(I,J,K) 

320  CONTINUE 

R2  =  R2  +  LFMT 
330  CONTINUE 
340  CONTINUE 
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C  R  =  TCTT  +  LFMT  +  2*LFMT*TC  +  2*(LF  -  »UMBU)*TC  +  TC 

C  VARIABLE  S<I,J): 

R  =  LFMT  +  TCTT  +  2*LFMT<-TC 

C*************’^*******************************************************fii***t***t 

C  IF  THE  USER  WAHTS  TO  ALLOW  A  TARGET  COMPLEX  TO  BE  ATTACKED  BY  TWO  LAUHCH 
C  FIELDS,  ADO  THE  FOLLOWIHG  ASSIGIMEMT  STATEHEKT: 

C  R2  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  MUMBU)*TC  +  3*TC 
C  *  -t-  WUMGOODOHUMFAIR  +  VUMFAIR^^IUMPOOR  5 

C****************************************************************************** 


C  STEP  THROUGH  EACH  LAUSCH  FIELD: 

DO  370  I  =  1,  LF 

R1  =  LFMT  +  TCTT  +  2*LFMT*TC  +  2*(LF  -  IUMBU)*TC 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  360  P  =  1,  NUMBU 

C  IF  THE  CURRENT  LAUNCH  FIELD  IS  A  BACKUP  ONE,  GO  TO  THE  NEXT  LAUNCH  FIELD: 

IF  (I  .Eq.  BU(P))  GO  TO  370 
360  CONTINUE 

C  STEP  THROUGH  EACH  TARGET  COMPLEX.  S(I,J)  WILL  BE  IN  CONSTRAINT  R  WITH  A 
C  COEFFICIENT  OF  1.0,  IN  CONSTRAINT  R  +  (LF  -  NUMBU) *TC  WITH  A  COEFFICIENT  OF 
C  THE  NUMBER  OF  MISSILE  TYPES,  AND  IN  CONSTRAINT  R1  WITH  A  COEFFICIENT  OF  1.0: 
DO  360  J  =  1,  TC 
R  =  R  +  1 
R1  =  R1  +  1 

WRITE(40,1060),  'S',  I,  J,  'C',  R,  1.0 
WRITE(40,1060),  'S',  I,  J,  'C',  R  +  (LF  -  HUMBU)*TC, 
ft  REAL(MT} 

WRITE(40,1060),  'S',  I,  J,  'C',  Rl,  1.0 

C****************************************************************************** 

C  IF  THE  USER  WANTS  TO  ALLOW  A  TARGET  COMPLEX  TO  BE  ATTACKED  BY  TWO  LAUNCH 
C  FIELDS,  ADD  THE  FOLLOWING  STATEMENT: 

C  WRITE(40.1060)  'S',  I,  J,  'C',  R2  +1,  1.0 


360  CONTINUE 
370  CONTINUE 


C****************************************************************************** 

C  M  =  LFMT  (OR  SOME  OTHER  LARGE  ENOUGH  NUMBER) 

C  R  =  LFMT  +  TCTT 

C  Rl  =  LFMT  +  TCTT  +  2*LFMT*TC 

C  R2  =  LFMT  +  TCTT  +  4*LFMT*TC 

C  DO  340  I  =  1,  LF 
C  DO  310  P  =  1,  NUMBU 

C  IF  (I  .EQ.  BU(P))  GO  TO  340 
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310  C08TIHUE 

DO  330  J  =  1.  TC 
R1  =  R1  +  1 
DO  320  K  =  1.  MT 

IF  (MADJCI.K)  .EQ.  0)  GO  TO  320 
ROW(I,J,K)  =  R1 
320  CONTIHUE 

330  COHTIHUE 

340  CONTINUE 

DO  376  I  =  1.  LF 

DO  350  P  =  1,  HUMBU 

IF  (I  .EQ.  BU(P))  GO  TO  376 
360  COHTIHUE 

R2  =  R2  -  LFHT*TC 
R3  =  0 

DO  370  J  =  1,  TC 
R2  =  R2  -  R3 
R3  =  0 

DO  360  K  =  1.  MT 

IP  (MADJd.K)  .EQ.  0)  GO  TO  360 
R  =  R  +  1 
B2  =  R2  +  1 
R3  =  R3  +  1 

WRITE(40,1050)  'Y',  I,  J,  K,  'C' 


R,  1.0 


C  WRITE(40,10S0)  'Y',  I,  J,  K,  'C»,  R  +  LFMT*TC, 

C  t  REAL(NUHTAR(J)) 

C  WRITE(40,1050)  'Y*.  I.  J,  K,  'C' ,  BOW(I,J.K),  REAL(M) 

C  DO  363  II  =  1,  LF 

C  DO  351  P  =  1,  HUMBU 

C  IF  (II  .EQ.  BU(P))  GO  TO  363 

C  351  COHTIHUE 

C  IF  (I  .EQ.  II)  GO  TO  363 

C  DO  352  KK  =  1,  MT 

C  IF  (MADJdl.KK)  .EQ.  0)  GO  TO  352 

C  WRITE(40jl060)  'Y',  I,  J,  K,  'C», 

C  *  1  ROWdl.J.KK),  1.0 

C  352  COHTIHUE 

C  353  CONTINUE  j 

C  WRITE(40,1060)  'Y',,  I,  J.  K,  »C',  R2,  -TCl.J.K) 

C  360  CONTINUE  \ 

C  370  COHTIHUE  \ 

C  375  CONTINUE  \ 

(;*^^,*i,^im*************************i‘********************’**m********************** 

C  VARIABLE  DELTA (J) : 

R  =  LFMT  +  TCTT  +  3+LFMT+TC  +  2*(LF  -  HUMBU)*TC  +  TC 
R1  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  HUMBU)*TC  +  2*TC 


C***********1‘****ii***************************************:t********rl,i,:HHi^,*t***** 

C  IF  THE  USER  DOES  NOT  WANT  TO  USE  3(1,1)  VARIABLES.  COMMENT  OUT  THE  ABOVE 
C  ASSIGNMENT  STATEMENTS  FOR  R  AND  R1  AND  IHCLUDE  THE  BELOW  ONES: 
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C  B  -  LF»r  +  TCTT  +  4'*LFMT*TC 
C  R1  =  R  +  TC 

C*'t<4>*<4<*i»**4>**4>*********^***********«*4'*«***i****************#****«*«**«**4>*«**** 

C  STEP  THROUGH  EACH  TARGET  COMPLEX.  DELTA(J)  VILL  BE  II  COHSTRAIIT  R  WITH  A 
C  COEFFICIEHT  OF  TKAX  AID  IN  COMSTRAIHT  R1  +  1  WITH  A  COEFFICIEIT  OF  1.0: 

DO  380  J  =  1,  TC 
R  =  R  +  1 

WRITE(40.1065),  'DELTA',  J,  'C',  R,  TMAX 
WRITE(40.1066),  'DELTA',  J,  'C',  R1  +  1,  1.0 
380  CONTIHUE 

WRITE (40, 1040),  'IHT' 

C  VARIABLE  Z(J): 

R1  =  LFMT  +  TCTT  +  2*LFMT*TC  +  2*(LF  -  IUI!BO)*TC  +  TC 

R2  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  IUMBU)*TC  +  TC 

R3  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  IUMBU)*TC  +  4*TC  +  3 

R4  =  0 

R6  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  HUMBU)*TC  +  4*TC  +  3 
*  +  HUMGOOD*BUMFAIR 

C  IF  THE  USER  DOES  HOT  WANT  TO  USE  S(I,J)  VARIABLES,  COMMENT  OUT  THE  ABOVE 
C  ASSIGNMENT  STATEMENTS  FOR  Rl,  R2,  R3,  R4,  AND  R5  AND  INCLUDE  THE  BELOW  ONES: 

C  Rl  =  LFMT  +  TCTT  +  3*LFMT*TC 
C  R2  =  Rl  +  LFMT+TC 

C  R3  =  R2  +  3*TC  +  3 

C  R4  =  0 

C  RS  =  R3  +  NUMGOOD*NUMFAIR 

C******t*********t**’l‘^*<^*******m***********‘*************t*****************^**** 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  440  J  =  1,  TC 

C  Z(J)  WILL  BE  IN  CONSTRAINTS  Rl  +  1  TO  Rl  +  LFMT  WITH  COEFFICIENTS  OF  1.0: 

DO  390  R  =  1,  LFMT 

WRITE(40,1070),  'Z' ,  J,  'C',  Rl  +  R,  1.0 
390  CONTINUE 

Rl  =  Rl  +  LFMT 

R2  =  R2  +  1 


C  Z(J)  WILL  BE  IB  CONSTRAINT  R2  WITH  A  COEFFICIEHT  OF  -1.0: 

WRITE (40, 1070),  'Z',  J,  'C',  R2,  -1.0 

C  Z(J)  WILL  BE  IN  CONSTRAINT  R2  +  TC  +  1  WITH  A  COEFFICIENT  OF  -1.0: 
WRITE(40,1070),  'Z',  J,  'C',  R2  +  TC  +  1,  -1.0 
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C  2(J)  WILL  BE  II  COHSTRAINT  R2  +  2*TC  +  1  WITH  A  COEFFICIEHT  OF  -1.0: 

WRITE (40, 1070),  '2',  J,  'C',  R2  +  2*TC  +  1,  -1.0 

C  IF  THERE  ARE  GOOD  TARGETS  AT  COMPLEX  J,  THE!  2(J)  WILL  CE  IH  COHSTRAIHTS 
C  R3  +  1  TO  R3  +  SUMFAIR  WITH  COEFFICIEITS  OF  1.0: 

IF  (IADJ(J,1)  .GT.  0)  THE! 

DO  400  R  =  .1,  lUHFAIR 
R3  =  R3  +  1 

WRITE(40,1070),  *2',  J,  'C*,  R3,  1.0 
400  COITIHUE 

C  ELSE,  IF  THERE  ARE  FAIR  TARGETS  AT  COMPLEX  J, 

ELSEIF  (HADJ(J,2)  .GT.  0)  THE! 

R4  =  R4  1 

R6  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  HlfflBU)*TC  +  4*TC  +3 
R7  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  HUMBU)*TC  +  4*TC  +3 

C*****************«*«*i»«*«««4<****««******  <******i)<#'»*«**<'4>4'^****#**«***4<*«****4i 

C  IF  THE  USER  DOES  lOT  WANT  TO  USE  S(I,J)  VARIABLES,  COMMENT  OUT  THE  ABOVE 
C  ASSIGNMENTS  FOR  R6  AND  R7  AND  INCLUDE  THE  BELOW  ONES: 

C  R6  =  LFMT  +  T(^rr  +  4*LFMT*TC  +  3*TC  +3 

C  R7  =  R6 

C****************************************************************************** 

C  THEN  STEP  THROUGH  EACH  COMPLEX  WITH  GOOD  TARGETS: 

DO  410  R  =  1,  NUMGOOD 

C  2(J)  WILL  BE  IH  CONSTRAINT  R6  +  R4  WITH  A  COEFFICIENT  OF  -1.0: 

WRITE(40,1070),  '2',  J,  'C’,  R6  +  R4,  -1.0 
R6  =  R6  +  NUMFAIR 
410  CONTINUE 

C  STEP  THROUGH  EACH  COMPLEX  WITH  POOR  TARGETS: 

DO  420  R  =  1,  NUMPOOR 

C  2(J)  WILL  BE  IN  CONSTRAINT  R7  +  NUMGOOD+NUMFAIR  +  R4  WITH  A  COEFFICIENT 
C  OF  1.0: 

WRITE(40,1070),  '2',  J,  *C',  R7  +  NUMGOOD*HUMFAIR  + 

A  R4.  1.0 

B7  =  R7  +  NUMFAIR 
420  CONTINUE 

C  ELSE,  COMPLEX  J  HAS  ONLY  POOR  TARGETS: 

ELSE 

C  STEP  THROUGH  EACH  COMPLEX  WITH  FAIR  TARGETS: 

DO  430  R  =  1 ,  NUMFAIR 
R5  =  R5  +  1 
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C  Z(J)  WILL  BE  II  COHSTRAIKT  RS  WITH  A  COEFFICIEIT  OF  -1.0; 

WRITE(40,1070),  ‘Z*.  J,  'C',  RS,  -1.0 
430  COITIHUE 

EID  IF 
440  COHTIIUE 

C  VARIABLE  THIRST: 

R  =  LFMT  ♦  TCTT  ♦  3*LFNT*TC  +  2*(LF  -  »UMBU)*TC  ♦  TC 


C#******«*<»****************4'«*«**«*********************«********************«* 

C  IF  THE  USER  DOEF  lOT  WAIT  TO  USE  THE  S(I.J)  VARIABLES.  COHMEHT  OUT  THE  ABOVE 
C  ASSIGKHENT  STATEHEIT  FOR  R  AID  IICLUOE  THIS  OIE: 

C  R  »  LFHT  ♦  TCTT  +  4*LFMT*TC 


C  STEP  THROUGH  EACH  TARGET  COMPLEX.  TFIRST  WILL  BE  IN  CONSTRAINTS  R  ♦  1  TO 
C  R  ♦  TC  WITH  COEFFICIEITS  OF  1.0: 

DO  460  J  »  1,  TC 
R  =  R  ♦  1 

WRITE(40,1080),  'TFIRST*,  'C',  R,  1.0 
4S0  CONTIIUE 

R  »  LFMT  ♦  TCTT  +  3*LFHT*TC  ♦  2*(LF  -  lUMBUj^TC  ♦  2*TC  ♦  1 

C  STEP  THROUGH  EACH  TARGET  COMPLEX.  TFIRST  WILL  PE  IN  CONSTRAINTS  R  +  1  TO 
C  R  ♦  TC  WITH  COEFFICIENTS  OF  1.0: 

DO  4S6  J  =  1,  TC 
R  =  R  +  1 

WRITE(40,1080),  'TFIRST',  'C',  R,  1.0 
455  CONTINUE 

R  =  LFHT  ♦  TCTT  +  3*LFMT*TC  +  2*(LF  -  NUMBU)*TC  ♦  4«TC  +  2 

C***********************************************************************^  «•*** 

C  IF  THE  USER  DOES  HOT  WANT  TO  USB  THE  S(I,J)  VARIABLES,  COMMENT  OUT  THE  ABOVE 
C  ASSIGNMENT  STATEMENT  FOR  R  AND  INCLUDE  THIS  ONE; 

C  R  =  LFMT  +  TCTT  +  4*LFMT*TC  ♦  3*TC  +2 

C***********************************^*****************************************. 

WRITE(40,1080),  'TFIRST',  'C,  R,  1.0 
WRITE{40,1080),  'TFIRST',  'C',  R  +  1,  -PCT 
WRITE(40,1080),  'TFIRST',  'C',  R  +  2  +  NUHGOOD*NUMFAIR  ♦ 

R  NUMFAIR*NUMPOOR,  1.0 

C  VARIABLE  TLAST: 

R  =  LFMT  +  TCTT  ♦  3*LFMT*TC  ♦  2*(LF  -  NUMBU)*TC  ♦  3*TC  +  1 
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C*********’*******'*** '********«  ***************4'************  **********  *********** 

C  IF  THE  USER  DOES  lOT  WANT  TO  USE  THE  S(I.J)  VARIABLES,  COHMEMT  OUT  THE  ABOVE 
C  ASSIGNMENT  STATEMENT  FOR  R  AND  INCLUDE  THIS  ONE: 

C  R  »  LFMT  +  TCTT  +  4*LFMT*TC  ♦  2*TC  *1 


C********************************************************* ********************* 


C  STEP  THROUGH  EACH  TARGET  COMPLEX.  TLAST  WILL  BE  IN  CONSTRAINTS  R  ♦  1  TO 
C  R  ♦  TC  WITH  COEFFICIENTS  OF  1.0; 

DO  460  J  »  1,  TC 
R  =  R  ♦  1 

WRITE(40,1080).  'UAST',  'C*.  R,  1.0 
460  CONTINUE 

WRITE(40,1080),  'TLAST',  'C',  R  ♦  1,  -1.0 
C  VARIABLE  DUR: 

WRITE(40.1080),  'DUR',  'C',  R  +  1,  1.0 
WRITE(40.1C80).  'DUR',  'C',  R  *  2,  1.0 
R1  a  R  +  4  ♦  NUMGOOD*NUMFAIR  ♦  H0MFA1R*SUMP00R 
WRITE(40,1080) ,  'DUR',  'C',  Rl,  1.0 

C  VARIABLE  01: 

WRITE(40.1090) ,  '01',  'OBJ',  W1 
WRITE(40,1096),  '01',  'C',  Rl  -  1,  -1.0 


C  VARIABLE  02; 

WRITE(40,1090),  '02',  'OBJ',  W2 
WRITE(40,1095),  '02',  'C',  Rl,  -1.0 


C  THE  LAST  SECTION  LISTS  THE  RIGHT  HAND  SIDE  VALUES  FOR  EACH  CONSTRAINT: 
WRITE(40,1040) ,  'RHS' 

Q  a  0 


C  STEP  THROUGH  EACH  TARGET  COMPLEX: 
DO  480  J  =  1,  TC 

C 'STEP  THROUGH  EACH  TARGET  TYPE: 

DO  470  L  a  1,  TT 


C  ;F  THERE  ARE  NO  TARGETS  OF  TYPE  L  AT  COMPLEX  J,  GO  TO  THE  NEXT  TARGET  TYPE: 
IF  (HADJ(J,L)  .EQ.  0)  GO  TO  470 
q  =  Q  +  1 

WRITE (40, 11 00),  'RHS',  'C' ,  Q,  REAL(NADJ(J,L)) 

CONTINUE 
48|0  CONTINUE 
q  =  TCTT 


C  STEP  THROUGH  EACH  LAUNCH  FIELD: 
DO  510  I  =  1,  LF 
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C  STEP  THROUGH  EACH  bACKUP  LAOHCH  FIELD: 

DO  400  P  »  I,  lUHBU 

C  IF  THE  CURREIT  UUICH  FIELD  IS  A  BACKUP  OIE,  GO  TO  THE  lEXT  LAUICB  FIELD: 

IF  (I  .EQ.  BU(P))  GO  TO  510 
490  com  HUE 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  600  K  ■  1,  KT 

C  IF  THERE  ARE  10  MISSILES  OF  TYPE  K  AT  LAUICH  FIELD  I,  GO  TO  THE  lEXT  MISSILE 
C  TYPE: 

IF  (MADJ(I.K)  .EQ.  0)  GO  TO  600 
Q  «  Q  +  1 

HRITE{40,1100).  'RHS',  ‘C*.  Q,  REAL(MADJ(I,K)) 

BOO  com  RUE 

610  COMTINUE 


c  IF  THE  USER  DOES  lOT  WAIT  TO  USE  S(I,J)  VARIABLES,  THEM  IICLUDE  THE 
C  FOLLOWING  SECTIOR: 

C  Q  «=  q  +  2*LFMT*TC 
C  DO  618  I  «  1,  LF 
C  DO  611  P  -  1.  HUMBU 

C  IF  (1  .EQ.  BU(P))  GO  TO  618 

C  611  COITIHUE 

C  DO  617  J  »  1,  TC 

C  DO  616  K  »  1,  MT 

C  IF  (MADJd.K)  .EQ.  0)  GO  TO  616 

c  q  >  q  «  1 

C  VRITE(40,1100).  >RHS',  'C',  Q,  M 

C  616  COITIHUE 

C  617  COITIHUE 

C  618  COITIHUE 


C'*******«**********************************************************4*«********* 


C************«****«************')>********************4>************«************v 

C  IF  THE  USER  DOES  lOT  WAIT  TO  USE  THE  S(I,J)  VARIABLES,  COMMEIT  OUT  THE 
C  FOLLOHIIG  SECTION; 

q  =  TCTT  ♦  LFMT  +  2*LFMT*TC  ♦  2*(LF  -  HUMBU) *TC 


c  q  =  q  + 

C  STEP  THROUGH  EACH  TARGET  COMPLEX.  THE  RIGHT  HAND  SIDE  OF  COISTRAIITS  Q  ♦  1 
C  TO  q  ♦  TC  WILL  BE  1.0: 

DO  620  J  =  1,  TC 

q  =  q  +  1 

VRITE(40,1100),  'RHS*,  'C',  Q.  1.0 
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C*********  *************'*•******************************************************* 

C  IF  THE  USER  WARTS  TO  ALLOW  A  TARGET  COMPLEX  TO  BE  ATTACKED  BT  TWO  LAUICK 
C  FIELDS,  COHMEHT  OUT  THIS  STATEMEHT  AKD  IHCLUDE  THE  lEXT: 

WRITE(40,1100).  'RHS*,  'C'.  Q,  1.0 
C  WRITE(40,1100)  'RHS*,  *C*,  Q,  2.0 

C***********'***^*************************************************************** 

520  CORTIRUE 


q  »  Tcrr  +  lfmt  +  3*lfmt*tc  +  2*(lf  -  rumbu)*tc  ♦  2*tc  ♦  i 


C*********'*************************«*****«*****4>**********4'*****4>***4i******4>«** 

C  IF  THE  USER  DOES  ROT  WANT  TO  USE  TEE  S(I.J)  VARIABLES,  COMMENT  OUT  THE 
C  ABOVE  ASSIGNMENT  STATEMENT  FOR  Q  AND  INCLUDE  THIS  ONE; 

C  q  =  TCTT  +  LFMT  ♦  4*LFMT*TC  +  TC  ♦  1 

WRITE (40, 1100),  'RHS',  *C',  q.  TC  -  1.0 

q  a  TCTT  +  LFMT  +  3*LFMT*TC  ♦  2*(LF  -  NUMBU)*TC  +  4*TC  +  3 
A  ♦  NUMGOOD*NUMFAIR  ♦  BUMFAIR*NUMPOCR 

C****************************************************************************** 

C  IF  THE  USER  DOES  NOT  WANT  TO  USE  THE  S(I,J)  VARIABLES,  COMMERT  OUT  THE 
C  ABOVE  ASSIGNMENT  STATEMENT  FOR  q  AND  INCLUDE  THIS  ORE: 

C  q  =  q  +  2*TC  +  2  +  NUMGOOD^NUMFAIR  +  NUMFAIR*HUMPOOR 

C****************************************************************************** 

WRITE(40,1100),  'RHS',  'C,  q  +  1,  TMIN 

C  IF  THE  USER  WARTS  TO  ALLOW  A  TARGET  COMPLEX  TO  BE  ATTACKED  BY  TWO  LAUNCH 
C  FIELDS,  INCLUDE  THESE  STATEMENTS: 

C  q  =  LFMT  +  TCTT  +  3*LFMT*TC  +  2*(LF  -  RUMBU)*TC  +  4*TC 
C  ft  NUMGOOD«NUMFAIR  *  RUMFAIR*NUMPOOR  *  6 

C  WRITE(40,1100)  'RHS*,  *C',  q  +  1,  TC  ♦  1 

C'****************************************************************************** 


C**a***a«***********«**4>**«*******a**«****«********a*«****4>aaa*aa*4>*a  ****«****« 

C  IF  THE  SOLUTION  IS  UNBALANCED,  ADD  THE  FOLLOWING  SECTION: 

C  q  -  LFMT  +  TCTT  +  3*LFMT«TC  •»■  2*(LF  -  NUMBU)*TC  +  4*TC 
C  It  *  NUMGOOD*RUNFAIR  *  HUMFAIR*NUMPOOR  *  6 
C  DO  640  I  =  1,  LF 

C  DO  S30  P  =  1,  NUMBU 

C  IF  (I  .EQ.  BU(P))  GO  TO  640 
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C  630  COXTIHUE 

C  Q  =  Q  +  1 

C  WRITE(4C,H00)  *RHS»,  »C'.  Q,  lUMTAR(JSTiR) 

C  640  COXTIIUE 

C*********************i»********«******tf*********  ■<*************•***•«**«******# 


C*************4'************4'*************************************************** 

C  IF  THE  USER  DOES  HOT  UAHT  EACH  LAU1CH  FIELD  TO  HAVE  A  PCrSU  BACKUP,  ADD  THE 
C  FOLLOWIHG  STATEMEHTS: 

C  q  «  LFMT  +  TCTT  +  3*LFirr*TC  +  2*(LF  -  HUMBU)*TC  ♦  4*TC 
C  *  +  HUMGOOD^tUMFAIR  +  ■UlffAIR*HUMPOOR  +6 

C  WRITE(4C,1100)  *RHS*,  'C*.  Q  ♦  1.  TOTMIS 

VR1TE(40,1120),  'EHDATA* 

CLOSE (40) 

C  THIS  SECTIOH  CREATES  TREi  SPECS  FILE  FOR  ZOOM. 


0PE*(UIIT=60,  FILE=*SILOATT.SPC»,  ACCESS=*SEQUEITIAL* , 
k  STATUS= 'UHKHOWm'' .  IOSTAT=IERROR,  ERR=992) 


HRITE(B0,1130), 


'BEoitK’ 


•MIIIHIZE' 


I 


■UMROWS  *  TCTT  ♦  LFmIt  ♦  3*LFMT*TC  +  2*(LF  -  HUM3U)*TC  +  4*TC 
k  *  lUNGOOD^ilUHFAIR  *  IUNFAIR*IUHPOOR  t  5 

lUMVARS  =  LFMT*TCTT  *  LFMT*TC  (LF  -  ■UHBU)*TC  +  2*TC  +  6 
■UMIHT  =  LF«T*TC  +  (LF  -  ■U1IBU)*TC  ♦  TC 

c***************'**********!'**************************************************** 

C  IF  THE  USER  DOES  HOT  WAHT  TO  USE  THE  S(I,J)  VARIABLES,  COMMEHT  OUT  THE 
C  ABOVE  ASSIGIMENT  STATENEbTS  FOR  HUNROUS,  lUNVARS,  AID  HUNIHT  ARD  INCLUDE 
C  THESE: 


C  HUMROWS  =  TCTT  +  LFMT  +  4*LFMT*TC  +  3*TC 
C  t  *  RUMGOOD^HUMFAIR  *  HUMFAIR^HUMPOOR  -f  5 

C  HUMVARS  t  LFMT*TCTT  +  LFMT^TC  +  2*TC  +  6 
C  HUMIHT  =  LFMT*TC  +  TC 


ITERLIM  =  1000000 
HAZNODES  -  100000 
NAZSAVE  =  20 

VRITE(60,1140),  'ROWS',  HUMROWS,  'COLUMNS',  HUMVARS,  'BRANCH', 
ft  'YES',  'BOUNDS',  'HONE',  'GAP',  TMAX,  'INTEGER', 

ft  HUMIHT,  'LIMIT',  ITERLIM,  'MAX  NODES',  MAXHODES, 

ft  'MAX  SAVE',  HAXSAVE,  'OBJECTIVE',  'OBJ', 

ft  'PRINT  CONTINUOUS  O',  'PRINT  LPl  O', 

ft  'PRINT  LP2  O',  'PRINT  HEURISTIC  O', 

ft  'PRINT  BRANCH  O',  'PRINT  TOUR  0 ' ,  'QUIT', 

) 
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ft  '10',  'EID* 

CL0SE(50) 

RETURI 

991  0PE«(UIIT=16.  FILE=*SILOATT.ERR»,  ACCESS='SEQUEITIAL' , 
ft  STATUS=»UIKIOWE') 

WRITE(16,*)  'ERROR  I»  OPEIIHG  FILE  SILOATT.MPS' 
VRITE(16.*)  'ERROR  CODE  =  '.  lERROR 
CLOSEdS) 

STOP 


992  0PEI(U1IT=1S,  FILE='S:L0ATT.ERR',  ACCESS =' SEQUEITIAL ' , 
ft  STATUS='U»KI0W1(') 

WRITEdS,*)  'ERROR  IH  OPERIHG  FILE  SILOATT.SPC' 
WRITEdB,*)  'ERROR  CODE  =  ',  lERRDR 
CLOSEdS) 

STOP  - 


1000  F0RMAT(A4, 
1010  F0RMAT(2X. 
1020  F0RMAT(2X, 
1030  F0RM>T(4X, 
1040  F0RMAT(A3) 
1060  F0RMAT(4X, 
1060  F0RMAT(4X, 
106S  F0RMAT(4X, 
1070  F0RMAT(4X, 
1080  F0RHAT(4X, 
1090  F0RHAT(4X, 
1096  F0RMAT(4X, 
1100  F0RMAT(4X, 
1110  F0RMAT(A7) 
1120  FORMAT (A6) 
1130  F0RMAT(A5, 
1140  F0RMAT(A4, 


IX,  A22,  /,  A4,  /,  2X,  Al,  IX,  A3) 

Al,  IX,  Al,  13) 

Al,  IX,  Al,  13,  /,  2X,  Al,  IX,  Al,  13) 

Al,  II,  II,  II,  II,  6X,  Al,  13,  6X,  F8.3) 


ft 

ft 

ft 

ft 


A6, 

AS, 

A9, 

A12 


/.  A8) 

IX,  14,  /,  A7,  IX,  14,  /,  A6,  IX,  A3,  /,' 

IX,  A4,  /,  A3,  IX,  F6.2,  /,  A7,  IX,  13,  /, 

IX,  17,  /,  A9,  JX,  16,  /,  A8,  IX,  13,  /, 
lx,  A3,  /,  A18,  /,  All,  /,  All,  /,  A17,  /,  A14,  /, 
,  /,  A4  ,  IX,  A2,  /,  A3) 


END 


C*****«*******>^*****************f***********#*'***'^****4i****4>*4'*'***<»«***«****4>* 

C******«*4>*****«*****4i**************«********«******4i«i^**4i*4>*«**«:****4i*4>****4>» 

SUBROUTIKE  PRMTSOL(BU,  CPLX,  FIELD,  LF,  MADJ,  HAXLF,  MAXHT, 
ft  HAXTC,  MAXTT,  MAXVAR,  MT,  MTYPE,  MADJ, 

ft  NUMBU,  T,  TC,  TT,  TTYPE,  VARVAL) 

C***************** 

c 

C  MAME:  PRMTSOL 
C 


91 


CCR  DESCRIPTIOI 


C  REV  DATE  OF 
C  10.  CBAIGE  DESIGIER 

C  -  -  - 

C  000  9FEB93  R.  PACE  ORIGIRAL  RELEASE 

C 

C  DESCRIPTIOI:  PRUTS  THE  SOLOTIUI  IITO  TWO  OUTPUT  FILES 

C 

C  CALLED  BY:  SILOATT 
C 

C  CALLS:  lOIE 
C 

C  CALLIIG  SEQUEICE: 

C 

CII  BU  -  ARRAY  STORIIG.TBE  BACKUP  LAUICH  FIELO:> 

C  CPU  -  ARRAY  ST  TUG  THE  lA-MES  OF  THE  TARGET  COMPLEXES 

C  FIELD  -  ARRAY  S'^GMIG  ThE  lAMES  OF  THE  UUHCH  FIELDS 

c  LF  -  HUMPER  OF  LAUICH  FIELDS 

C  MADJ  -  ARRAY  STORIIG  THE  ADJUSTED  lUMBER  OF  MISSILES  BY  TYPE  AT  EACH 

C  UUICH  FIELD 

C  NAXLF  -  MAX  lUNBER  OF  LAUHCH  FIELDS 

C  MAXMT  -  MAX  lUMBER  OF  MISSILE  TYPES 

C  MAXTC  -  MAX  lUMBER  OF  TARGET  COMPLEXES 

C  MAXTT  -  MAX  lUMBER  OF  TARGET  TYPES 

C  MAXVAR  -  MAX  lUMBER  OF  VARIABLES 

C  MT  -  HUMBER  OF  MISSILE  TYPES 

C  MTYPE  -  ARRAY  STORIIG  THE  lAMES  OF  TEE  MISSILE  TYPES 

C  MADJ  -  ARRAY  STORIIG  THE  ADJUSTED  lUMBER  OF  TARGETS  BY  TYPE  AT  EACH 

C  TARGET  COMPLEX 

C  HUMBU  -  THE  lUMBER  OF  BACKUP  LAUICH  FIELDS 

C  T  ARRAY  STORIIG  FLIGHT  TIMES  OF  MISSILES  (BY  TYPE)  FROM  EACH  LAUHCH 

C  FIELD  TO  EACH  TARGET  COMPLEX 

C  TC  -  HUMBER  OF  TARGET  COMPLEXES 

C  TT  -  HUMBER  OF  TARGET  TYPES 

C  TTYPE  -  ARRAY  STORIIG  THE  lAMES  OF  THE  TARGET  TYPES 

C  VARVAL  -  ARRAY  STORIIG  THE  SOLUTIOI  VECTOR 
C 

C  FILES: 

C 

C  SILOATT. SOL 
C  SIHBAC.SOL 
C  SILOATT. ERR 
C 

C  LOCAL  VARIABLES: 

C 

C  DELTA  -  ARRAY  STORIHG  THE  SOLUTIOI  VALUES  OF  THE  DELTA(J)  VARIABLES 
C  DUR  -  STORES  THE  SOLUTIOH  VALUE  L.  ‘^HE  VARIABLE  DUR 

C  I,  II,  J,  K,  L,  P,  Q  -  COUHTERS 

C  TERROR  -  STORES  THE  VALUE  OF  lOSTAT 

C  S  -  ARRAY  STORIIG  THE  SOLUTIOI  VALUES  OF  THE  S(I.J)  VARIABLES 

C  TFIRST  -  STORES  THE  SOLUTIOI  VALUE  OF  THE  VARIABLE  TFIRST 

C  TLAST  -  STORES  THE  SOLUTIOI  VALUE  OF  THE  VARIABLE  TLAST 


C  TOTMIS  -  ARRAY  STORIHG  THE  TOTAL  HUMBER  OF  MISSILES  OF  TYPE  K  THAT  ARE 

C  ALLOCATED  FROM  LAUHCH  FIELD  I  TO  TARGET  COMPLEX  J 

C  X  -  ARRAY  STORIHG  THE  SOLUTIOH  VALUES  OF  THE  X(I,J,K.L)  VARIABLES 

C  Y  -  ARRAY  STORIHG  THE  SOLUTIOH  VALUES  OF  THE  Y(1,J,K)  VARIABLES 

C  Z  -  ARRAY  STORIHG  THE  SOLUTIOH  VALUES  OF  THE  2(J)  VARIABLES 

C 

C  ERROR  MESSAGES:  HOHE 
C 

Cm*******^**t-*************************************************t***************** 

C 

C  LOCAL  VARIABLE  DECLARATIOKS: 

IHTEGER  BU(MAXLF),  DELTA(IO),  I,  lERROR,  II,  J,  X,  L,  LF, 
k  MADJ(MAXLF,MAXMT),  MT.  KADJ(MAXTC,MAXTT) ,  lUMBU,  P,  Q, 

k  S(IO.IO),  TC.  TT,  Y(10,10.5) 

REAL  DUR.  T(MAXLF.MAXTC,MAXMT),  TFIRST,  TLAST, 

k  TOTMISdO.lO.B),  X(i0,10,S,5),  Z(10) 

DOUBLE  PRECISION  VARVAL(MAXVAR) 

CHARACTER'S  CPLX(MAXTC),  FIELD (MAXLF) 

CHARACTERS  MTYPE(MAXMT) ,  TTYPE(MAXTT) 

C 

C  THIS  SECTION  ASSIGNS  THE  VALUES  OF  THE  ZOOM  VARIABLES  IHTO  THE 
C  APPROPRIATE  VARIABLES  FOR  THE  HIP. 

C 

q  =  0 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  60  I  =  1,  LF 

C  IF  IT  IS  OHE  OF  THE  BACKUP  ONES,  GO  TO  THE  HEXT  LAUHCH  FIELD: 

DO  10  P  =  1 ,  NUHBU 

IF  (I  .EQ.  BU(P))  GO  TO  50 
10  CONTINUE 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  40  J  =  1,  TC 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  30  K  =  1,  MT 

C  IF  THERE  ARE  NO  MISSILES  OF  TYPE  K  AT  LAUHCH  FIELD  I,  GO  TO  THE  HEXT 
C  MISSILE  TYPE: 

IF  (MADJ(I,K)  .EQ.  0)  GO  TO  30 

C  STEP  THROUGH  EACH  TARGET  TYPE: 

DO  20  L  =  1,  TT 
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C  IF  THERE  ARE  10  TARGETS  OF  TYPE  L  AT  COMPLEX  J,  GO  TO  THE  HEXT  TARGET  TYPE: 

IF  (IADJ(J,L)  .EQ.  0)  GO  TO  20 

C  OTHERWISE,  SET  X(I,J,K,L)  EQUAL  TO  THE  lEXT  VALUE  II  THE  SOLUTIOI  VECTOR: 

Q  =  q  +  1 

X(I,J,K,L)  =  VAKVaL(Q) 

20  comiuE 

30  corriiuE 

40  COHTIHUE 

SO  COHTIIUE 

C  STEP  THROUGH  EACH  LAUICH  FIELD: 

DO  90  I  =  1.  LF 

fc  IF  IT  IS  ONE  OF  THE  BACKUP  DIES,  GO  TO  THE  lEXT  LAUNCH  FIELD: 

DO  60  P  =  1,  NUMBU 

IF  (I  .EQ.  BU(P))  GO  TO  90 
60  CONTINUE 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DC  80  J  =  1,  TC 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  70  K  =  1,  MT 

C  IF  THERE  ARE  HO  MISSILES  OF  TYPE  K  AT  LAUNCH  FIELD  I.  GO  TO  THE  NEXT 

C  MISSILE  TYPE: 

IF  (MADJd.K)  .EQ.  0)  GO  TO  70 

C  OTHERWISE,  SET  Y(I,J,K)  EQUAL  TO  THE  NEXT  VALUE  II  THE  SOLUTIOI  VECTOR: 

q  =  q  +  1 

Y(I,J,K)  =  VARVAL(q) 

70  COITINUE 

80  CONTINUE 
90  CONTIIUE 


C*******************^***************'*>******’t'****'************'*>^*************4i*** 

C  IF  THE  USER  DOES  NOT  WAIT  TO  USE  THE  S(I.J)  VARIABLES,  THEN  pO  NOT  DECLARE 
C  S(I,J)  (SEE  ABOVE)  AID  COMMENT  OUT  THE  FOLLOWING  SECTION: 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  120  I  =  1,  LF 


C  IF  IT  IS  ONE  OF  THE  BACKUP  ONES,  GO  TO  THE  NEXT  UUNCH  FIELD: 
DO  100  P  =  1,  lUMBU 

IF  (I  .EQ.  BU(P))  GO  TO  120 
100  CONTINUE 


C  STEP  THROUGH  EACH  TARGET  COMPLEX: 
DO  110  J  =  1,  TC 
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C  SET  S(I,J)  EQUAL  TO  THE  HEXT  VALUE  IH  THE  SOLUTIC:  ’TCTOR: 
q  =  Q  +  1 

S(I,J)  =  VARVAL(Q) 

110  COITIIUE 
120  CORTINUE 


C  STEP  THROUGH  TARGET  COMPLEX: 

DO  130  J  =  1,  TC 

C  SET  DELTA (J)  EQUAL  TO  THE  MEXT  VALUE  II  THE  SOLUTIOI  VECTOR: 

Q  =  Q  +  1 

DEI.T^CJ)  =  VARVAL(q) 

130  COITIIUE 

C  STEP  THROUGH  EACH  TARGET  COMPLEX: 

DO  140  J  =  1,  TC 

C  SET  2(J)  EQUAL  TO  THE  lEXT  VALUE  II  THE  SOLUTIOI  VECTOR: 

Q  =  Q  +  1 
Z(J)  =  VARVAL(q) 

140  COITIIUE 

C  SET  TFIRST,  TLAST,  ASO  DUR  EQUAL  TO  THE  lEXT  THREE  VALUES  II  THE  SOLUTIOI 
C  VECTOR: 

TFi;ST  s  VARVAL(Qt-l) 

TLAST  =  VARVAL(Q+2) 

DUR  =  VARVAL(Q+3) 

C  THIS  SECTIOI  PRODUCES  THE  OUTPUT  FOR  THE  J-8  SOLITION  FILE: 

0PEN(UHIT=16,  F1LE=’ SILO ATT. SOL',  ACCESS= 'SEQUENTIAL' , 

A  STATUS='UHKNOWH'.  IOSTAT=IERROR,  ERR=991) 

WRITE(16,1000) 

C  STEP  THROUGH  EACH  LAUNCH  FIELD: 

DO  170  I  =  1,  LF 

C  STEP  THROUGH  EACH  BACKUP  LAUNCH  FIELD: 

DO  160  P  =  1,  lUMBU 

C  IF  THE  CURRENT  LAUNCH  FIELD  IS  FOR  BACKUP  ONLY,  THEI  "SLIDE"  THE  lANE  OF 
C  EACH  SUBSEQUENT  LAUNCH  FIELD  FORVARD  IN  THE  LIST  (THEREBY  ERASING  THE  NAME 
C  OF  THE  BACKUP  FIELD: 

IF  (I  .EQ.  BU(P))  THEN 
DO  ISO  II  =  I,  LF 

FIELD(II)  =  FIELD(II+1) 

GO  TO  170 

ISO  CONTINUE 

END  IF 
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160  COHTIMUE 
170  COHTHUE 

LF  =  LF  -  lUMBU 

WRITE(16,1010)  (FIELDCI),  1=1, LF) 

C  STEP  THROUGH  EACH  TARGET  CCMPLEX: 

DO  210  J  =  1,  TC 

C  STEP  THROUGH  EACH  TARGET  TYPE: 

DO  200  L  =  1,  TT 

C  IF  THERE  ARE  EG  -fARGETS  OF  TYPE  L  AT  COMPLEX  J,  CO  TO  THE  lEXT  TARGET  TYPE 
IF  (IADJ(J,L)  .EQ.  0)  GO  TO  200 

C  OTHERWISE,  WRITE  THE  lAME  OF  THE  CURREHT  TARGET  COMPLEX  TO  'SILOATT.SOL' , 

C  AHD  THEH  STEP  THROUGH  EACH  LAUHC3  FIELD: 

WRITE(16.1020)  CPLX(J) 

DO  190  I  =  1,  LF 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  i80  K  =  1,  MT 

C  IF  HO  MISSILES  OF  TYPE  K  HAVE  BEEH  ALLOCATED  FROM  UUHCH  FIELD  I  TO  TYPE 
C  L  TARGETS  AT  TARGET  COMPLEX  J,  THEH  GO  TO  THE  lEXT  MISSILE  TYPE: 

IF  (X(X,J,K,L)  .EQ.  0)  GO  TO  180 

C  THIS  IF  STATEMENT  BLOCK  WRITES  THE  SOLUTIOH  DATA  UHDER  THE  COLUMN  FOR  THE 
C  APPROPRIATE  LAUNCH  FIELD: 

IF  (I  .EQ.  1)  THEN 

IF  (Z(J)  .GE.  TFIRST)  THEH 
Wr.ITE(16,1030)  TTYPE(L),  MTYPE(K), 

*  INT(I(I,J,K,L)),  Z(J) 

ELSE 

HRITE(16,1030)  TTYPE{L) ,  MTYPE(K), 

*  irr(X(I,J,K,L)),  TFIRST 
END  IF 

ELSE  IF  (I  .EQ.  2)  THEN 
IF  (Z(J)  .GE.  TFIRST)  THEN 
HRITE(16,1040)  TTYPE(L),  MTYPE(K), 
k  IHT(X(I,J,K.L)),  Z(J) 

ELSE 

WRITE(16,1040)  TTYPE(L),  HTYPE(K), 
k  IIT(X{I,J,K,L)),  TFIRST 

END  IF 

ELSE  IF  (I  .EQ.  3)  THEN 
IF  (Z(J)  .GE.  TFIRST)  THEN 
WRITE(16,10S0)  TTYPE(L),  MTYPE(K), 
k  INT{X(I,J,K,L)),  Z(J) 

ELSE 

WRITE(16,1050;  -TYPECD ,  HTYPE(K) , 

*  INT(X(I,J,K,L)),  TFIRST 
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EID  IF 
ELSE 

IF  (Z(J)  .GE.  TFIRST)  THEI 
WRITE(16,1060)  TTYPE(L),  MTyPE(K) . 

*  IlTCXd.J.K.D),  Z(J) 

ELSE 

WRITE(16,1060)  mPE{L),  MTYPEd) , 
ft  IIT(X(I.J,K,L)),  TFIRST 

EID  IF 
EID  IF 

180  COITIIUE 

190  COITIIUE 
200  COITIIUE 
210  COITIIUE 
C 

C  THIS  S'lCTIOI  PRODUCES  THE  OUTPUT  FOR  THE  SIIBAC  SOLUTIOI  FILE: 

C 

0PEI(UIIf=17.  FILE=* SIIBAC. SOL*,  ACCESS= 'SEQUEITIAL* , 
ft  STATUS=*UIKIOWH’.  IOSTAT=IERROR,  ERR=992) 

C  STEP  THROUGH  EACH  LAUICR  FIELD: 

D0240I,=1,LF 

C  STEP  THROUGH  EilCH  TARGET  COMPLEX: 

DO  2S0  J  =  1,  TC 

C  STEP  THROUGH  EACH  MISSILE  TYPE  AID  SET  TOTMISd, J,K)  EQUAL  TO  0: 

1^0  220  X  =  1,  NT 
!  TOTMISd, J.K)  =  0 
220  (ioiTIIUE 

230  COITilUE 
240  COITIIUE 

I 

I 

C  STEP  THROUGH  I EACH  TARGET  COMPLEX: 

DO  280  J  i=  1,  TC 

C  SUM  UP  THE  TOTAL  lUMBER  OF  MISSILES  OF  TYPE  K  ATTACKIHG  COMPLEX  J  FROM 
C  LAUICH  FIELD  I: 

C  STEP  THROUGH  EACH  LAUICH  FIELD: 

DO  270  I  =  1,  LF 

C  STEP  THROUGH  EACH  MISSILE  TYPE: 

DO  260  K  =  1,  MT 

C  STEP  THROUGH  EACH  TARGET  TYPE: 

DO  260  L  =  1,  TT 

C  IICREASE  THE  CURREIT  VALUE  OF  TOTMISd, J,K)  BY  THE  VALUE  OF  Xd,J,K,L) 
TOTMISd, J,K)  =  TOTMISd, J,K)  +  X(I,J,K.L) 

250  COITIIUE 
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C  IF  THERE  ARE  HISSILES  OF  TYPE  K  ALLOCATED  FROM  LAUHCH  FIELD  I  TO  TARGET 
C  COMPLEX  J.  WRITE  THE  SOLUTIOM  DATA  TO  'SIIBAC  SOL’: 

IF  (TOTMIS(I.J.K)  .CT.  0.0)  THEM 
IF  (T0TM1S(I,J,K)  .'JE.  10)  TBEI 

WRITEdT.lOrO)  CPLX(J)(1:3),  MTYPE(K) , 
ft  FIELD(I)(1:3),  IirT(10*T0TMIS(l,J.K)) 

ELSE  IF  (TOTMISd.J.K)  .CE.  1)  TREH 

WRITEdT.lOYl)  CPLX(J)d:3),  MTTPE(K)  , 
ft  FIELD(l)d;3),  lMTdO*TOTHIS{I , J.K)) 

ELSE 

WRITEd7,1072)  CPLX(J)d:3),  MTYPE(K), 
ft  FIELD(l)d;3).  IHTdO*TOTMIS(I,J.K)) 

EHD  IF 
EID  IF 
COHTIIUE 
27b  COHTIIUE 
280  COHTIIUE 

CLOSECIO) 

CLOSEdT) 

HETURI 

801  0PEH(0IIT«1S,  FILE»»SILOATT.ERR»,  STATUS* ’UHKHOWI’) 

WRITEdS,*)  'ERROR  II  OPEIIMC  FILE  SILOATT.SOL* 

WRITEdS,*)  'ERROR  CODE  «  lERROR 
CLOSECIE) 

002  0PEI(UIIT*1B,  FILE- 'SILOATT. ERR',  STATUS- ’UIKIOWI’) 

WRITEClS,*)  'ERROR  II  OPEIIHG  FILE  SIIBAC. SOL’ 

VRITEdS,*)  'ERROR  CODE  >  '.  lERROR 
CLOSECIB) 

1000  F0RMAT(33X,  'LAUICI  FIELDS',  /,  'TOT  FLOS’,  20X, 
ft  '(lUMBER/TYPE/IMPACT  TIME)',  /,  IIX) 

1010  FORMATdlX,  10(A8,  9X)) 

1020  F0RNAT(A8) 

1030  FORMATdX,  A4,  61,  A4,  '/’,  13,  F8.3) 

1040  FORMATdX,  A4,  23X,  A4,  '/’,  13,  ’/’,  F#.3) 

1060  FORMATdX,  A4,  40X,  A4,  '/',  13,  F8.3) 

1060  FORMATdX,  A4,  67X,  A4,  '/’,  13,  '/’,  F6.3) 

1070  F0RMAT(A3,  A4,  A3,  13) 

1071  F0RMAT(A3,  A4,  A3,  'O',  12) 

1072  F0RMAT(A3,  A4,  A3,  '00',  11} 

EID 


Appendix  C.  Output  Files 

C.  1  Solution  Summary  for  J-8  Review 

The  following  is  the  output  file  produced  by  the  FORTRAN  program  which  is 
in  the  format  that  J-8  wants  so  that  they  can  review  the  solution: 


LAUNCH  7iELDS 


TGT  ELDS 

(TYPE/NUMBER/ rMP ACT  TIME) 

NEW  BERN 

DURHAM 

RALEIGH 

WILMINGT 

MACON 

GOOD 

GOOD/  5/35.150 

GOOD 

ATL/TYB 

FAIR/  14/35.150 

GOOD 

GOOD/ 

3/35.980 

GOOD 

FAIR/ 

15/35.980 

ATL/TYB 

FAIR 

SAVANNAH 

GOOD 

COLUMBUS 

FAIR 

BRUNSWIG 

GOOD 

ATHENS 

GOOD/ 

18/35.980 

FAIR/  21/34.027 

FAIR/  12/37.430 

GOOD/  24/37.430 

POOR 

FAIR/  2/37.430 

C.8  Solution  Summary  for  SINBAC 


'  The  following  is  the  output  file  produced  by  the  FORTRAN  program  which  is 
in  the  format  that  SINBAC  needs  in  order  to  analyze  the  solution: 


MACG00DWILO5O 
HACFAIRWIL140 
ATLG00DNEW210 
ATLFAIRNEWj  50 
SAVFAIRDUR210 
C0LFAIRRAL12O 
BRUGO0DWIL24O 
ATHFAIRRAL020 


Appendix  D.  GAMS  Input  File 


The  following  is  the  GAMS  input  file  that  was  used  to  solve  the  unclassified 
sample  problem: 


$OFFSYMXREF  OFFSYMLIST 
SETS 

I  launch  fields  /  N-BERN,  DURHAM,  RALEIGH,  WILM  / 

J  target  fields  /  MACON,  ATL-TYB,  SAVAN,  COLUMB,  BRUNS,  ATHENS  / 

K  missile  types  /  GOOD,  FAIR  / 

K2(K)  missile  types  /  GOOD,  FAIR  / 

L  target  types  /  GOOD,  FAIR,  POOR  / 

G(J)  fields  with  good  tgts  /  MACON,  ATL-TYB,  SAVAN,  BRUNS  / 

F(J)  fields  with  fair  tgts  /  COLUMB  / 

P<J)  fields  with  poor  tgts  /  ATHENS  /; 


TABLE  T(I,K,J)  flight  times 


MACON 

ATL-TYB 

SAVAN 

COLUMB 

BRUNS 

ATHENS 

N-BERN. GOOD 

34.16 

35.98 

33.89 

35.33 

36.49 

38.22 

N-BERN. FAIR 

32.61 

33.24 

32.66 

33.10 

33.61 

34.82 

DURHAM. GOOD 

0 

0 

0 

0 

0 

0 

DURHAM. FAIR 

32.11 

32.405 

32.45 

32.51 

32.96 

34.71 

RALEIGH. GOOD 

0 

0 

0 

0 

0 

0 

RALEIGH. FAIR 

32.95 

33.475 

33.09 

33.42 

33.92 

35.34 

WILM. GOOD 

35.15 

34.84 

35.60 

36.26 

37.43 

40.78 

WILM. FAIR 

32.30 

32.565 

32.67 

32.70 

33.14 

34.94; 

TABLE  M(I,K)  number  of  missiles  of  type  K  at  field  I 


GOOD 

FAIR 

N-BERN 

21 

15 

DURHAM 

0 

22 

RALEIGH 

0 

15 

WILM 

29 

17; 

TABLE  N(J,L)  number  of  targets  of  type  L  at  field  J 

GOOD  FAIR  POOR 
19  0  0 


MACON 


ATL-TYB  18 

18  0 

SAVAN  21 

0  0 

COLUHB  0 

12  0 

BRUNS  24 

0  0 

ATHENS  0 

0  2  ; 

SCALARS  TMAX 

maximum  flight  time  possible 

/40.78/ 

THIN 

minimum  flight  time  possible 

/32.11/ 

PCT 

user-specified  percentage 

/l.O/ 

TF 

number  of  target  fields 

/6/ 

.  Q 

most  launch  fids  that  pan  hit 

each  tgt  fid 

MISSNUN 

the  number  of  missile  types 

/2/ 

W1 

relative  weight  of  obj  1  /lO/ 

W2 

relative  weight  of  obj  2  /!/  ; 

VARIABLES 

X(I,J,K,L) 

num  of  msls  of  type  k  from  i  to  tgts  of  type  1 

Y(I,J,K) 

indicator  Vciriable  for  if  msls 

of  type  k  from 

S(I,J) 

indicator  variable  for  if  msls 

allocated  from 

/I/ 


Z(J) 

DEL(J) 

TFIRST 

TLAST 

DUR 

OPT 

SIPLUS 

S2PLUS 

SIHINUS 

S2MINUS 


largest  flight  time  of  missiles  allocated  to  j 
indicator  for  disjunctive  constraints 
flight  time  of  the  missile  that  impacts  first 
flight  time  of  the  missile  that  impacts  last 
duration  of  the  attack  (tlast  -  tfirst) 
the  optimal  solution 

is  over  its  goal 
is  over  its  goal 
is  under  its  goal 
is  under  its  goal  ; 


i  to  j 


the  amount  that  obj  1 
the  amoTint  that  obj  2 
the  amount  that  obj  1 
the  amount  that  obj  2 


FREE  VARIABLE  OPT 


POSITIVE  VARIABLES  X,  Z,  TFIRST,  TLAST,  DUR,  SIPLUS,  S2PLUS,  SIHINUS, 
S2MINUS  ; 


BINARY  VARIABLES  S,  Y,  DEL  ; 


EQUATIONS 

FLTTIME 

OBJl 

0BJ2 

DEMAND (J,L) 

SUPPLY (I, K) 

IND1(I,J,K) 

IND2(I,J,K) 

IND3(I,J) 

IND4(I,J) 

MAXFT(I,J,K) 

DISJUNl(J) 

DISJUN2 


objective  function  | 
goal  for  tfirst  \ 

goal  for  dur  \ 

attack  each  target  \ 

cannot  shoot  more  than  aVailable 
set  yijk  to  1  if  missiles  of  type  k  are 
allocated  from  i  to  j 
set  sij  to  1  if  missiles 
and  to  0  otherwise 
identify  the  max  fit  time  into  j 
disjunctive  constraints  to  find  the  minimum 
flight  time 


allocated  from  i  to  j  , 
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TFIR(J)  ensure  tfirst  is  set  to  the  smallest  z(j) 

TLAS(J)  find  the  maximum  flight  time 

DURl  duration  equals  tlast  minus  tfirst 

DUR2  duration  must  be  .le.  a  user-specified  percent  of  tfirst 

HAXHIT(J)  do  not  attack  each  tgt  field  with  more  than  1  launch  fid 

PRTY1(G,F)  attack  good  tgts  before  fair  ones 

PRTY2(F,P)  attack  fair  tgts  before  poor  ones; 

FLTTIHE  . .  OPT  »E=  W1*S1PLUS  +  H2*S2PLUS  ; 

OBJl  . .  TFIRST  ♦  SIMINUS  -  SIPLUS  =E=  TMIN  ; 

0BJ2  . .  DUR  +  S2MINUS  -  S2PLUS  =E=  0  ; 

DEMAND(J,L)$(N(J,L))  ..  SUH((I,K)$(M(I ,K)) ,  X(I,J,K,L))  =E=  N(J,L)  ; 
SUPPLY(I.K)$(M(I,K))  ..  SUM((J,I.)$(N(J,L)),  X(I,J.K,L))  =L=  H(I,K)  ; 
IND1(I,J,K)$(T(I,K,J))  ..  Y(I.J,K)  =L=  SUM(L$(N(J,L)) ,  X(I,J,K,L))  ; 
IND2(I,J,K)$(T(I,K,J))  ..  Y(I,J,K)  =G=  (SUM(L$(N(J,L)) ,X(I. J,K,L)))/ 
(SUM(L,N(J,L)))  ; 

IND3(I,J)  ..  S(I,J)  =L=  SUH(K$(M(I,K)),  Y(I,J,K))  ; 

IND4(I,J)  S(I,J)  =G=  (SUM(K$(H(I,K)),  Y(I , J,K)))/MISSNUH  ; 
MAXFT(I,J,K)$(T(1,K,J))  ..  Z(J)  =G=  Y(I . J,K)*T(I,K, J)  ; 

DISJUNl(J)  ..  TFIRST  =G=  Z(J)  +  DEL(J)*v-TMAX)  ; 

DISJUN2  ..  SUM(J,  DEL(J))  =E=  TF  -  1  ; 

TFIR(J)  . .  TFIRST  =L=  Z(J)  ; 

TLAS(J)  ..  TLAST  =G«  Z(J)  ; 

DURl  . .  DUR  =E*  TLAST  -  TFIRST  ; 

DUR2  . .  DUR  eL=  PCT*TFIRST  ; 

MAXHIT(J)  ..  SUM(I,  S(I,J))  =E=  Q  ; 

PRTYKG.F)  ..  Z(G)  =L*  Z(F)  ; 

PRTY2(F,P)  ..  Z(F)  =L=  Z(P)  ; 

MODEL  ALLOCATE  /ALL/; 

OPTION  ITERLIH  =  1000000  ; 

OPTION  RESLIH  =  50000  ; 

OPTION  WORK  =  100000  ; 

OPTION  OPTCR  =  0.01  ; 

OPTION  LIHROW  =  0  ; 

OPTION  LIMCOL  =  0  ; 

SOLVE  ALLOCATE  USING  MIP  MINIMIZING  OPT  ; 
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Appendix  E.  Complete  Solutions  and  Allocation  Summaries 
Case  1:  See  chapter  IV,  section  4.3,  p.  41. 


Case  2:  Pci  =  0.05,  Wi  =  10,  IVj  =  1 
SOLUTION; 

Objective  Function  Value  =  37.159 


'T/irit 

= 

35.648 

Ttast 

= 

37.43 

Dur 

= 

1.782 

= 

3.538 

O2 

= 

1.782 

- 

= 

3 

-^1212 

= 

18 

X1221 

15 

•^2121 

= 

19 

^2623 

= 

2 

X3422 

= 

12 

-^4311 

21 

^"4511 

= 

8 

■^4521 

= 

16 

Yul 

= 

1 

^122 

= 

1 

^2,2 

1 

^262 

= 

1 

^42 

= 

1 

5^31 

= 

1 

X 

^451 

= 

1 

^52 

= 

1 

5i2 

= 

1 

S21 

= 

1 

^26 

= 

1 

•S'34 

= 

1 

S43 

= 

1 

S45 

= 

1 

= 

37.43 

Z2 

= 

35.98 

Z3 

35.648 

^4 

= 

37.43 

Zs 

= 

37.43 

Zs 

=z 

37.43 

= 

1 

62 

sr 

1 

64 

1 

^5 

= 

1 

1 

all  other  variables  =  0. 
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Case  3:  Pet  =  0.1,  =  2,  1^2  =  1 

SOLUTION: 

Objective  Function  Value  =  7.237 


Pfirtt 

34.027 

= 

37.43 

Dur 

= 

3.403 

0, 

= 

1.917 

O2 

= 

3.403 

^1211 

= 

3 

X1212 

= 

18 

^1221 

= 

15 

-^2321 

= 

21 

^3422 

= 

12 

XiUl 

= 

19 

^4sn 

= 

10 

-^4521 

= 

14 

■^4623 

= 

2 

Vi21 

=: 

1 

il22 

= 

1 

L232 

= 

1 

^3.2 

=i 

1 

= 

1 

L451 

= 

1 

^452 

= 

1 

^462 

= 

5i2 

=: 

1 

S23 

= 

1 

‘S’34 

= 

1 

S4I 

1 

S45 

= 

1 

■S’46 

= 

1 

Zi 

= 

35.15 

Z2 

= 

35.98 

^3 

= 

34.027 

Z4 

= 

37.43 

Zs 

= 

37.43 

^6 

= 

37.43 

Si 

= 

1 

S2 

= 

1 

^4 

— 

1 

Ss 

= 

1 

Se 

= 

1 

all  other  variables  = 

0. 

- 
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ALLOCATION  SUMMARY: 


Table  9.  Case  3  Allocation  Summary 


LAUNCH 

FIELD 

MISSILE 

TYPE 

TARGET  COMPLEX 

TARGET  TYPE 

Macon 

Atl/Tyb 

Savannah 

Columbus 

Bcaseswick 

Athens 

good 

good 

fair 

good 

fair 

good 

poor 

New  Bern 

good 

fair 

3 

15 

18 

Durham 

fair 

21 

Raleigh 

fair 

12 

V/ilmington 

good 

fair 

19 

10 

14 

2 

Charlotte 

Lumbcrton 
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Case  4:  Pet  =  0.1,  Wi  =  1,  VVj  =  2 
SOLUTION: 

Objective  Function  Value  =  5.32 


Tjirti 

= 

37.43 

I'last 

= 

37.43 

Dur 

= 

0 

Oi 

= 

5.32 

O2 

= 

0 

-^1211 

= 

3 

^1212 

= 

18 

-^1221 

= 

15 

■^2121 

= 

19 

■^2623 

= 

2 

-^3422 

= 

12 

-^'*311 

= 

21 

-^4511 

— 

8 

^4521 

= 

16 

y\i\ 

= 

1 

5^122 

=: 

1 

T212 

— 

1 

1^62 

=z 

1 

5^42 

= 

1 

T431 

= 

1 

^451 

= 

1 

^452 

= 

1 

S\3 

= 

1 

521 

1 

526 

= 

1 

Sm 

1 

543 

1 

•^45 

= 

1 

Zi  =  37.43  Zj  =  37.43  Z3  =  37.43 

Z4  =  37.43  Z5  =  37.43  Zg  =  37.43 


ALLOCATION  SUMMARY: 


Table  10.  Case  4  Allocation  Summary 


LAUNCH 

FIELD 

MISSILE 

TYPE 

TARGET  COMPLEX 

TARGET  TYPE 

Macon 

Atl/Tyb 

Savannah 

Columbus 

Bcaseswick 

Athens 

good 

good 

good 

fair 

good 

poor 

New  Bern 

good 

fair 

3 

IS 

18 

Durham 

■fair 

19 

2 

Raleigh 

fair 

12 

Wilmington 

good 

fair 

■ 

21 

8 

16 

Charlotte 

■ 

Lumberton 

Case  5;  Pci  =  0.05,  VKi  =  2,  =  1 

SOLUTION: 

Objective  Function  Value  =  8.858 


TfiTMt 

= 

35.648 

Ttast 

= 

37.43 

Dur 

= 

1.782 

0, 

= 

3.538 

O2 

= 

1.782 

■^1211 

= 

3 

-^1212 

= 

18 

^1221 

15 

-^2121 

= 

19 

-^2623 

= 

2 

^3.22 

= 

12 

=: 

5 

-^.321 

= 

16 

-^4511 

= 

24 

== 

1 

^122 

= 

1 

5^212 

=: 

1 

^262 

= 

1 

y342 

= 

1 

5^31 

= 

1 

^432 

= 

1 

^51 

= 

1 

5i2 

=; 

1 

521 

= 

1 

526 

1 

53. 

= 

1 

5.3 

= 

1 

5.5 

1 

= 

35.648 

^2 

= 

35.98 

^3 

= 

35.648 

37.43 

^5 

= 

37.43 

37.43 

Sj 

= 

1 

Sz 

— 

1 

<5. 

1 

6, 

= 

1 

6e 

= 

1 

all  other  variables  =  0. 
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ALLOCATION  SUMMARY: 


Table  11.  Case  5  Allocation  Summary 


LAUNCH 

FIELD 

MISSILE 

TYPE 

TARGET  COMPLEX 

TARGET  TYPE 

Macon 

Atl/Tyb 

Savannah 

Columbus 

Bcaseswick 

Athens 

good 

good 

fair 

good 

fair 

Kss&ai 

poor 

New  Bern 

good 

fair 

3 

15 

18 

Durham 

fair 

19 

2 

Raleigh 

fair 

s 

12 

Wilmington 

good 

fair 

■ 

■ 

5 

16 

24 

Charlotte 

■ 

■ 

Lumberton 

Case  6;  Pet  =  0.05,  W'l  =  1,  W^2  =  2 
SOLUTION: 

Objective  Function  Value  =  5.32 


Tfirit 

= 

37.43 

Tiaat 

= 

37.43 

Dur 

= 

0 

Ox 

=■ 

5.32 

O2 

= 

0 

■^1211 

= 

3 

^1212 

= 

18 

^\22\ 

15 

-^2321 

=z 

21 

^3422 

= 

12 

-^4111 

= 

5 

^<121 

=r 

14 

^45X1 

= 

24 

>^4623 

= 

2 

5^121 

1 

^122 

=r 

1 

I232 

1 

5^42 

= 

1 

= 

1 

^4,2 

= 

1 

K,51 

= 

1 

^62 

= 

1 

5i2 

=' 

1 

■S’23 

= 

1 

«S34 

== 

1 

5^1 

.= 

1 

545 

= 

1 

S4G 

= 

1 

Zx 

= 

37.43 

^2 

= 

37.43 

Z3 

=: 

37.43 

= 

37.43 

= 

37.43 

Ze 

37.43 

Sx 

= 

1 

52 

= 

1 

^3 

— 

1 

S4 

= 

1 

Ss 

=: 

1 

all  other  variables  =  0. 
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Case  7:  Pet  =-•  0.1,  IV,  =  1,  W',  =  1 
SOLUTION: 


Objective  Function  Value 

37.43  r,„t  =  37.43 

=  5.: 

Dur 

Oi 

= 

5.32 

0, 

= 

0 

'^1211 

= 

3 

-^1212 

18 

■^122: 

^2121 

1- 

19 

A'3412 

= 

12 

■^362; 

'^4.11 1 

= 

21 

■^4511 

= 

8 

-^4521 

^121 

1 

T122 

= 

1 

VjlJ 

V342 

= 

1 

^2 

= 

1 

5^31 

K|51 

= 

1 

ns2 

= 

1 

= 

1 

Six 

= 

1 

534 

•Ssf. 

= 

1 

543 

= 

1 

Sif, 

^1 

=r 

37.43 

Zi 

= 

37.43 

Z3 

= 

37.43 

Z5 

= 

37.43 

Ze 

= 

1 

^2 

= 

1 

^3 

^4 

1 

^5 

= 

1 

all  other  variables  =  0. 


ALLOCATION  SUMMARY: 


Table  14.  Case  8  Allocation  Summaxy 


LAUNCH 

FIELD 

MISSILE 

TYPE 

TARGET  COMPLEX 

TARGET  TYPE 

Macon 

Atl/Tyb 

Savannah 

Columbus 

Bcaseswick 

Athens 

good 

good 

fair 

good 

fair 

good 

poor 

New  Bern 

good 

fair 

3 

15 

18 

Durham 

fair 

21 

Raleigh 

fair 

12 

2 

Wilmington 

good 

fair 

5 

14 

■ 

■ 

24 

Charlotte 

1 

Lumberton 

! 

■■■ 

■■ 

■■■■ 

■■■■I 

HHHHI 

■■■ 

Case  9:  Pet  =  0.15,  Wi  =  10,  =  1 

SOLUTION: 

Objective  Function  Value  =  9.26 


TjiTSt 

== 

32.548 

Tla,t 

= 

37.43 

Dur 

= 

4.882 

Oi 

= 

0.438 

O2 

= 

4.882 

-^1211 

== 

3 

■^1212 

= 

18 

^1221 

= 

15 

-^2321 

21 

Jl3422 

= 

12 

-^4111 

= 

19 

-^4511 

= 

10 

-^4521 

= 

14 

-^4623 

= 

2 

1 

^122 

= 

1 

I232 

= 

1 

^342 

= 

1 

= 

1 

I45I 

= 

1 

^452 

= 

1 

V462 

1 

S\2 

1 

523 

1 

534 

= 

1 

Sa\ 

= 

1 

S45 

= 

1 

546 

= 

1 

Z, 

= 

35.15 

Z2 

= 

35.98 

^3 

32.548 

^4 

= 

37.43 

z^ 

= 

^17.43 

^6 

= 

37.43 

<5i 

=: 

1 

S2 

1 

^4 

1 

^5 

= 

1 

Se 

=  ■ 

'1 

all  other  variables  =  0. 


117 


ALLOCATION  SUMMARY: 


Table  15.  Case  9  Allocation  Summary 


LAUNCH 

FIELD 

MISSILE 

TYPE 

TARGET  COMPLEX 

TARGET  TYPE 

Macon 

Atl/Tyb 

Savannah 

Columbus 

BcMeswick 

Athens 

good 

good 

fair 

good 

fair 

good 

poor 

New  Bern 

good 

fair 

3 

15 

18 

Durham 

fair 

*■ 

21 

Raleigh 

lair 

12 

Wilmington 

good 

fair 

19 

10 

14 

2 

Charlotte 

Lumberton 


Appendix  F.  User’s  Guide 

F.  1  Introduction 

This  appendix  contains  instructions  for  using  the  FORTRAN  program  to  solve 
a  missile  allocation  problem. 

F.2  How  to  Use  SILOATT 

1.  Make  sure  that  the  files  SILOATT.EXE  and  SILOATT.DAT  are  loaded  in  the 
current  directory  on  the  VAX/VMS  computer  system. 

2.  Edit  the  file  SILOATT.DAT: 

(a)  Input  the  data  for  number  of  missiles  by  type  at  each  launch  field  into 
the  first  table  following  these  instructions: 

i.  Each  row  of  this  table  must  have  the  following  format:  launch  field 
name  (columns  1-8),  one  blank  space,  missile  type  (columns  10-13), 
four  blank  spaces,  and  number  of  missiles.  The  name  of  the  launch^ 
field  must  start  in  column  1  and  the  missile  type  must  start  in  column 
10.  If  the  name  of  a  launch  field  (or  missile  type)  is  less  than  eight 
(four)  characters,  put  in  enough  blanks  to  make  the  entire  field  eight 
(four)  characters.  If  the  name  of  a  launch  field  (or  missile  type)  is 
more  than  eight  (four)  characters,  input  no  more  than  eight  (four) 
characters.  The  number  of  missiles  may  begin  in  column  18  or  be¬ 
yond. 

Example: 

New  Bern  good  12 

New  Bern  fair  22 

Durham  fair  33 

Lumberto  good  20 

ii.  The  first  launch  field  listed  will  be  considered  by  the  program  as 
launch  field  1  from  then  on,  the  second  one  listed  will  be  considered 
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launch  field  2,  and  so  on.  In  the  example  above,  New  Bern  is  launch 
field  1,  Durham  is  launch  field  2,  and  so  on. 


iii.  Similarly,  the  first  missile  type  listed  will  be  considered  missile  type 
1,  the  second  one  listed  as  missile  type  2,  and  so  on.  In  the  example 
above,  good  is  missile  type  1  and  fair  is  missile  type  2. 

iv.  The  numbers  of  missiles  entered  must  be  integer  values  only. 

(b)  Input  the  data  for  number  of  targets  by  type  at  each  target  complex  into 
the  second  table  following  these  instructions: 


i.  Each  row  of  this  table  must  have  the  following  format:  target  complex 
name  (columns  1-8),  one  blank  space,  target  type  (columns  10-13), 
four  blank  spaces,  and  number  of  targets.  The  name  of  the  target 
complex  must  begin  in  column  1  and  the  target  type  must  begin  in 
column  13.  If  the  name  of  a  target  complex  (or  target  type)  is  less 
than  eight  (four)  characters,  put  in  enough  blanks  to  make  the  entire 
field  eight  (four)  characters.  If  the  name  of  a  target  complex  (or 
target  type)  is  more  than  eight  (four)  characters,  input  no  more  than 
eight  (four)  characters.  The  number  of  targets  may  begin  in  column 
18  or  beyond. 


Example: 

Macon  good  120 

Macon  fair  220 

Savannah  fair  330 

Athens  good  20 

ii.  The  first  target  complex  listed  will  be  consid^ered  by  the  program  ns 
target  complex  1  from  then  on,  the  second  one  listed  will  be  consid¬ 
ered  target  complex  2,  and  so  on.  In  the  example  above,  Macon  is 
target  complex  1,  Savannah  is  target  complex  2,  and  so  on. 
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iii.  Similarly,  the  first  target  type  listed  will  be  considered  target  type 
1,  the  second  one  listed  as  target  type  2,  and  so  on.  In  the  example 
above,  good  is  target  type  1  and  fair  is  target  type  2. 

iv.  The  numbers  of  targets  entered  must  be  integer  values  only. 

(c)  Input  the  data  for  the  flight  times  of  each  missile  type  from  each  launch 
field  to  each  target  complex  into  the  third  table  by  following  these  in¬ 
structions: 

i.  The  first  row  of  this  table  must  contain  the  names  of  the  target  com¬ 
plexes.  Input  the  names  starting  in  column  18  (leave  the  rows  with 
the  headings  in  please,  including  the  headings  “LAUNCH  FIELD” 
and  “MISS  TYPE”).  Each  target  complex  name  should  be  eight  char¬ 
acters  long  followed  by  one  blank  space.  These  names  must  appear 
exactly  £is  in  the  second  table  and  be  in  the  same  order  (but  list  each 
target  complex  only  once  here). 

ii.  Leave  one  blank  line  and  then  type  in  the  rest  of  the  rows  in  the 

following  format:  launch  field  name  (columns  1-8),  one  blank  space, 

missile  type  (columns  10-13),  four  blank  spaces,  and  the  flight  time 

of  the  appropriate  type  of  missile  from  the  current  launch  field  to 

each  target  complex.  The  names  of  the  launch  fields  must  begin 

in  column  1  and  the  missile  types  must  begin  in  column  10.  The 

flight  times  should  be  no  larger  than  999.99  and  should  be  rounded 

off  to  the  nearest  one-hundredth  of  a  minute.  Each  flight  time  entry 

should  stcirt  directly  under  the  ^irst  character  of  the  appropriate  target 

complex  and  should  take  up  six  spaces  followed  by  three  blank  spaces. 

Example: 

LAUNCH  MISS 

FIELD  TYPE  MACON  SAVANNAH  ATHENS 

New  Bern  good  34.55  35,11  39.99 

New  Bern  fair  35.66  36.92  40.82 
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Durham  fair  47,33  43.22  45,99 

Lumberto  good  37.82  32.94  IB. 39 

iii.  Note  that  each  launch  field  and  each  miss'  .  ,ype  must  appear  exactly 

the  same  as  it  does  in  the  first  table  and  be  in  the  exact  same  order. 

(d)  Input  the  value  of  the  weighting  factors  for  minimizing  the  earliest  flight 

time  (M^i)  and  duration  These  must  be  decimal  form  (real  num¬ 

bers). 

(e)  Input  the  percentage  of  the  earliest  flight  time  that  the  duration  is  allowed 
to  be  (Pet).  This  must  be  in  decimal  form  (i.e.,  0.1  instead  of  10%). 

(f)  Input  the  percent  backup  (Petbu)  that  is  required  to  be  in  each  launch 
field.  This  must  also  be  in  decimal  form. 

(g)  Indicate  whether  any  of  the  launch  fields  should  be  reserved  for  backup 
only  by  typing  ‘Y’  or  ‘y'  for  yes  and  ‘N’  or  ‘n’  for  no.  If  there  are  such 
launch  fields,  indicate  how  many  and  which  ones  (list  by  number,  not  by 
name). 

(h)  Indicate  whether  any  of  the  tairget  complexes  should  be  combined  by 
typing  ‘Y’  or  ‘y’  for  yes  and  ‘N’  or  ‘n’  for  no.  If  any  should  be  combined, 
list  which  ones  (smallest  number  first)  and  the  name  of  the  combined 
complex  (maximum  of  eight  characters). 

3.  At  the  command  prompt,  type 

def  for029  siloatt.spc 
def  for024  siloatt.mps 

which  tells  the  computer  that  the  files  siloatt .  spe  and  siloatt  .mps  should 
be  equated  to  FORTRAN  units  29  and  24  respectively.  This  needs  to  be  done 
because  ZOOM  accesses  these  units  directly  without  using  OPEN  and  CLOSE 
statements. 
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4.  At  the  command  prompt,  type 


run  siloatt 

which  causes  the  program  to  execute. 

5.  The  program  creates  the  following  output  files: 

SILOATT. SOL  -  Contains  the  allocation  summary. 

SILOATT. ERR  -  Contains  appropriate  error  messages.  Only  created  if 

-  there  was  a  problem  encountered  during  execution. 

SINBAC.SOL  -  Contains  the  allocation  in  the  SINBAC  format. 

FOR025.DAT  -  Created  by  ZOOM.  Contains  the  entire  solution  to  the 

-  mixed-integer  program. 

6.  If  a  feasible  solution  is  found,  the  program  creates  the  above  mentioned  .SOL 
files.  If  this  solution  is  unsatisfactory,  the  user  may  want  to  modify  the  solution 
directly  (just  make  sure  the  new  solution  is  still  feasible).  Otherwise,  the 
user  could  either  modify  the  input  data  and  re-run  the  program  or  modify  a 
constraint  by  changing  something  in  the  FORTRAN  code.  If  the  FORTRAN 
code  is  changed,  be  sure  io  re-compile  the  entire  program  (including  all  the 
ZOOM  files)  before  trying  to  re-run  it.  Also,  make  sure  to  record  the  changes 
made  so  that  they  can  be  un-made  later! 

7.  If  no  feasible  solution  exists,  an  error  message  stating  so  is  written  in  file 
SILOATT. ERR.  At  this  point,  the  user  could  again  choose  either  to  modify 
the  input  ^iata  file  or  the  FORTRAN  code  and  re-run  the  program. 
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